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Abstract 



NEC FASCES 
NEC OPES 
SOLA ARTIS 
SCEPTRA PERENNANT 
Not power, not wealth, 
only the reign 
of art and science 
shall persist 
Tycho Brahe 1546 - 1601 

This thesis deals with different aspects of a special kind of high redshift galaxy, 
namely Lya emitters. Lya emitters are galaxies found through their Lya emission, 
at redshifts larger than z >2 where the emission line has been redshifted into the 
optical or near-infrared regime. The thesis has two main parts; a lower redshift, 
observational part [z ~ 3) and a more technical/theoretical very high redshift part 
{z ~ 9). 

In the first, lower redshift part I present the analysis of a narrow-band image 
taken in the GOODS-S field, focused on a redshift for Lya oi z = 3.15. The image, 
covering a central ~ 7' x 7' part of the GOODS-S field revealed 25 Lya emitting 
candidates, of which one turned out to be a so-called Lya "blob". Lya blobs are 
large nebulae of gas emitting a large amount of light in the Lya line. They can be 
very bright {Liy Q, ~ lO^'^ erg s ^) and have projected diameters as large as 150 kpc. 
Three possible mechanisms have been proposed to explain this phenomenon; i) star 
formation and "super-winds", ii) AGN activity and Hi) cold accretion. The blob in 
GOODS-S turned out to be the first Lya blob best explained by cold accretion, and 
was published in a Letter to the Editor in A&A in June 2006 (Nilsson et al. 2006a). 
The remaining 24 Lya emitters were analysed together, in order to gain insight 
in the nature of galaxies selected through narrow-band imaging for redshifted Lya 
emission. An extensive SED fitting was performed and the results also included 
the discovery of an apparent filamentary structure and a comparison to Lyman- 
Break Galaxies at similar redshifts. The results from this analysis is accepted for 
publication in A&A in Nilsson et al. (2007). 

In the second part, I discuss future, very high redshift narrow-band surveys for 
Lya emitters. In February 2005, the idea of buying a narrow-band filter for the, 
then half way through its construction. Visible and Infrared Survey Telescope for 
Astronomy (VISTA) to search for Lya emitters at very high redshift started. In 
Chapter H] I describe the science case for the filters and the process of designing and 
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procurement of the filters. In this chapter is also described how the initial idea for 
a narrow-band survey was developed into what has become ELVIS - Emission Line 
galaxies with VISTA Survey, part of the Ultra- VISTA survey which will - hopefully 
- start observing in the Spring of 2008. The need for better predictions of what 
may be observed with this type of survey was realised early on, and in Chapter O I 
present a paper in which we try to make more accurate predictions of the outcome 
of very high redshift narrow-band surveys for Lya emitters using two theoretical 
models. 

Finally, in a project unrelated to the other two parts of the thesis, I present 
a search for a "Fundamental plane" of Lya emitters in the colour space produced 
by large-scale multi-wavelength surveys such as GOODS or COSMOS. The goal of 
the project was to decide the most efficient narrow-band colour selection method for 
detecting Lya emitters at several different redshifts. We also wanted to optimise the 
selection method to exclude interloper galaxies such as lower redshift [OII]-emitters. 
In this thesis I present the results of this study. 
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Chapter 1 
Introduction 



1.1 The high redshift Universe 

1.1.1 Brief history of our Universe and the redshift 

Our Universe was created in the Big Bang. In the first few fractions of a second the 
Universe was extremely hot, but it expanded and cooled rapidly. After this initial 
stage of "primeval soup" , where particles and radiation were coupled and high energy 
physics governed everything, a period of cooling began. At around SOO'OOO yrs after 
the Big Bang, the Universe had cooled and de-pressurised enough so that nuclei 
had formed and the electrons had been captured to make neutral atoms. This is 
called recombination, see Fig. ll.li At this time, the photons created in the Big Bang 
and immediately thereafter were released and could begin their journey through the 
Universe. This radiation is called the Cosmic Microwave Background (CMB) and 
is an almost perfect blackbody emission with a peak temperature of 2.73 K. The 
irregularities seen in the CMB indicate that there were slight inhomogeneities in 
the matter density at recombination. These inhomogeneities later turned into the 
structures we see today. 

After recombination, these peaks in the matter-density distribution continued 
to contract due to gravity and after some ~ 100 Myrs, the first stars were born. 
These were massive stars that burned their fuel quickly, exploded and enriched 
their surroundings. Soon, more stars were forming and galaxies started to take 
shape. Also, the first quasars may have lit up at this time. These galaxies and 
quasars produced a large amount of UV radiation and this radiation field ensured 
that the Universe would once again become ionised and the electrons separated 
from the atoms. This is called the re-ionisation and is believed to have occurred at 
z ^ 6 — 10, between 0.5 — 1.0 Gyrs after the Big Bang. The search for when exactly 
this happened is a hot topic today, as our telescopes push further and further back 
in time. The current record in redshift for an observed galaxy today is ^ = 6.96 (lye 
et al. 2006), and observations in the near future should be able to push this record 
to redshifts around 2; ~ 9, see also Chapter H] and O 

From redshift ~ 6 and on, the Universe has continued to expand, cool, gather 
mass in clusters and galaxies and produce stars. The redshift at which the star 
formation rate per volume peaked appears to be around redshift z ~ 1—2, see 
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Fig. I1.2[ Our Sun was formed some 4 Gyrs ago, corresponding to redshift z ~ 0.4. 

So, what is redshift? Redshift occurs as photons travel through the Universe, 
which simultaneously expands. An analogy is the change in pitch when a trains 
whistle or an ambulances sirens pass your ear. The sound waves are compressed 
in front of the train, thus the pitch goes up, and elongated behind the train, thus 
the decrease in pitch again. In space, the size of the Universe changes with time 
according to the scale factor a(t). Thus, photons travelling towards us will stretch 
its wavelength due to the expanding Universe and the light will be shifted towards 
the red, i.e. redshift. Redshift is defined as: 



-. ^observed '^(^2=o) /"I 1^ 

■^emitted (^\J' emitted ) 



1.1.2 Methods of finding high redshift galaxies 

In our constant search for knowledge we want to understand how the stars and 
galaxies were formed, and how they evolved. We can make models of star formation, 
but ultimately we need observations of the young Universe. Thus, we search for high 
redshift galaxies, whose light was emitted a long time ago, corresponding to temitted 
in Eq. II. 1[ That means that the light we observe from a galaxy at redshift 2; ~ 6 
was emitted from a very young galaxy, when the Universe was only ~ 1 Gyrs old. 
Looking further away in redshift means looking further into our past. And thus 
finding a large, representative sample of galaxies at different redshifts will give us 
insight into galaxy formation and evolution. 

There are many methods of finding star forming galaxies in the high redshift 
Universe, and each method explores different classes of galaxies. One of the major 
questions in observational cosmology is thus what are the relations between these 
classes. This thesis explores properties of high redshift galaxies selected by their 
Lya emission (see below). In this introduction I discuss several different methods 
used to detect high redshift galaxies. 

Lyman Break Galaxies selection 

One of the most common methods in the last decade, now comprising an impressive 
catalogue of thousands of galaxies, has been the selection based on photometric 
redshifts gained from observations of the "Lyman Break" in galaxies. The Lyman 
Break, located at 912 A in the spectrum of a galaxy, represents the cut-off energy 
where the single electron of a hydrogen atom is ionised. Almost all the UV light 
of a galaxy is absorbed at wavelengths shorter than this, as hydrogen is extremely 
abundant in galaxies. Thus, a "break" in the spectrum of the galaxy. The Lyman 
Break Galaxy (LBG) selection method was pioneered by Steidel et al. (1996; 1999; 
2000; 2003), but has also been used by many other groups (e.g. Madau et al. 1996; 
Pettini et al. 2001; Bunker et al. 2004; Stanway et al. 2004; Ouchi et al. 2004a,b; 
Wadadekar et al. 2006). An illustration of the method is found in Fig. II. 3[ As can 
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Figure 1.1: This plot illustrates the history of our Universe and 

marks a few note- worthy redshifts, see text. From www. astronomy. ohio- 
state.edu/^pogge/TeachRes/Artwork/Cosmology/index. html. 



be seen, a galaxy will appear to "drop out" in the bluest filter, hence the method 
is sometimes also referred to as the drop-out technique. Spectroscopic follow-up is 
necessary to confirm the high redshift nature of the galaxy. 

The Lyman Break technique can be used for a wide range of redshifts, when 
different filters are used for the selection. Table 11.11 gives the filters and redshifts 
commonly used, or proposed to be used. The technique is typically spectroscopically 
complete to an R band magnitude of i? ~ 25.5 for U-band drop-outs but the rate 
of confirmation falls for higher redshift LBG candidates (Giavalisco 2002). A few 
studies of the properties of LBGs such as masses, dust content, ages etc. have 
been made so far. These studies are well summarised in Giavalisco (2002). LBGs 



4 



Chapter 1 



1 234567 




3 0.4 06 0.8 



Figure 1.2: Star formation rate density history of the Universe. Data points are star 
formation rate densities from observed galaxies. Lines are best fit extrapolations. 
From Hopkins & Beacom (2006). 




Figure 1.3: Illustration of the LBG selection method. The top panel shows the 
typical shape of a galaxy spectrum, redshifted to ^ ~ 3. The Lyman Break can 
be seen at approximately 400 nm. The bottom panels show how such a galaxy 
would be seen as observed through the three broad-band filters U, G and R that are 
located on the blue and the red side of the break respectively. The galaxy will not 
be observed in the filter blueward of the break, but is clearly seen in the red filters. 
Image credit: Johan Fynbo. 
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Table 1.1: Redshifts and filters used in the LBG method. 



Blue filter 


Red filters 


Redshift range 


U 


B, G, R, V 


2.5 - 3.5 


B 


G, R, V 


3.5 - 4.5 


V 


i, z 


4.5 - 6.0 


i 


z, J 


6.0 - 7.5 


z 


J, H, K, 


7.5 - 11.5 


J 


H, K, 


11.5 - 15.5 



appear to have ages ranging from a few to several hundred Myrs or even up to 1 
Gyr. Stellar masses lie in the range 10^ < M^/Mq < 10^^. Papovich, Dickinson & 
Ferguson (2001) find very high extinction in their sample of LBGs, with Ay ^1 — 2, 
while Shapley et al. (2001) and Verma et al. (2007) find more modest values of 
Ay ~ 0.3 — 0.5 in their samples. Thus, it appears that LBGs are medium mass 
and medium dusty galaxies, with high star formation rates of several hundred solar 
masses per year (Shapley et al. 2001). Several groups have also detected clustering 
in LBGs, similar to that of Lya emitters (Steidel et al. 1998; Giavalisco et al. 1998; 
Ouchi et al. 2004b; see also sec. 11.3.51) . 

Damped Ly« Absorption selection 

When observing a quasar at high redshift, sometimes a gas cloud with a high col- 
umn density will be found along the same line of sight. This is not as uncom- 
mon as one might think. When a cloud with an HI column density larger than 
N{H1) > 2 X 10^° cm~^ is observed along the sight-line of a quasar it is called a 
Damped Lya Absorber (DLA) as its Lya absorption line will be damped. Thus, 
the DLA technique is based on finding high column density galaxies by searching 
for absorption lines in the spectra of quasars. An illustration of the method can be 
found in Fig. II. 4[ The search for DLA systems began already in the mid-80s but a 
large sample of DLA galaxies were not collected until ten years ago (for a review, 
see Wolfe, Gawiser & Prochaska 2005). 

A problem with this selection method is the apparent proximity of the quasar 
in the sky, making imaging follow-up difficult. However, DLA system can still give 
valuables insights into properties such as the neutral gas fraction in the Universe, 
chemical evolution and metallicity production across a large redshift range and gas 
kinematics in galaxies. 

Selection based on sub-millimeter emission 

The Spectral Energy Distribution (SED) of a dusty starburst galaxy will have its 
peak emission at wavelengths of ~ 0.1 mm, see Fig. 11.51 These galaxies will, however, 
not have significant emission in the optical or near-infrared parts of the spectrum. 
Thus, observing the sky in the sub-mm will reveal a class of galaxies otherwise left 



6 



Chapter 1 




5000 5500 6000 6500 7000 
Wavelength {A) 



Figure 1.4: Illustration of the DLA selection method. The figure shows the Keck/ESI 
spectrum of QSO PSS0209+0517 at redshift zqso = 4.17. Two DLA systems can 
be observed at redshifts zdla = 3.86 and zdla = 3.67. The labelling refers to 
absorption lines identified to be associated with the zdla = 3.86 system. From 
Wolfe, Gawiser & Prochaska (2005). 



unstudied. Due to the difficulties in building detectors for these wavelengths, it was 
not until the commissioning of the SCUBA detector (Holland et al. 1999) on the 
James Clerk Maxwell Telescope (JCMT) in 1997 that progress was made in this field 
of study. As the sample of sub-mm selected galaxies has grown to be ~ 100 (e.g. 
Blain et al. 2002; Ivison et al. 2005), the median redshift of this sample appears to 
be z ~ 2, although sub-mm emission should in theory be just as efficient to detect 
galaxies between z = 1 — 10. 

Two, related, problems with selecting galaxies based on sub-mm emission are 
that the beam sizes are still very large (of the order 10 arcsec) and that identifi- 
cation of optical/infrared counterparts are difficult. The process of identifying a 
counterpart is generally done by first searching for a radio counterpart, since radio 
observations have smaller beam sizes and radio emissions have been shown to be 
proportional to sub-mm emission (e.g. Carilli & Yun, 1999). Following identification 
of a radio counterpart, an optical or infrared counterpart may be found. However, 
even with radio identifications, optical/infrared identifications can be difficult to 
make. The sample of galaxies with confirmed counterparts display a wide variety of 
properties including AGN and starburst activity, as well as signs of recent mergers 
(Ivison et al. 2000; Blain et al. 2002). Future instruments such as SCUBA-II and 
ALMA will hopefully give us a better insight into the nature of sub-mm selected 
galaxies. 
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Figure 1.5: Composite SED of a starburst galaxy. The peak at sub-mm wavelengths 
is apparent. From Blain et al. (2002). 



Gamma Ray Burst selection 

Gamma Ray Bursts (GRBs) are violent explosions, more energetic than supernovae. 
The name derives from the fact that GRBs are found through a burst of gamma 
ray emission. GRBs are typically divided into two classes of events, long or short 
duration bursts; where the long bursts are believed to be related to relativisti- 
cally beamed extreme versions of (hundreds of times brighter than) supernovae (e.g. 
Hjorth et al. 2003), i.e. the death throws of a massive star, and the short bursts 
are believed to occur when two compact objects (i.e. neutron stars or black holes) 
collide (e.g. Hjorth et al. 2005). The method of finding high redshift galaxies using 
GRBs is to use the location of the GRB to search for its host galaxy. The gamma 
ray event of the GRB triggers a satellite based telescope, such as for instance the 
Swifl0 telescope, to slew to the part of the sky where the event occurred. In general, 
the accuracy of the positioning from the gamma ray telescope is poor, e.g. of the 
order of arcminutes, and it is followed up with an X-ray telescope, as X-rays are 
easier to pin-point and most GRB afterglows emit strongly also in X-rays. When 
the location of the GRB has been localised to a smaller error box, optical/near-IR 
telescopes take over and search for the optical/near-IR afterglow of the GRB. Thus, 
the GRB can finally be completely localised, to subarcsecond resolution, using an 
optical detection. After the GRB has faded, this position can be observed to greater 
depth to search for the host galaxy of the GRB. Host galaxies of GRBs are often 
blue, star-forming galaxies (e.g. Bloom et al. 1998; Sokolov et al. 2001; Gorosabel 
et al. 2005; Fruchter et al. 2006). 

Lya selection 

A very efficient method to detect high redshift galaxies is to observe the sky with 
a narrow-band filter focused on the Lya emission line at a particular redshift. A 

"'^http: / / www.swift.psu.edu/ 
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narrow-band filter is a filter that allows only a very small range of wavelengths to 
pass through it and blocks all other light. This will result in a narrow range of 
redshifts for Lya, typically Az ^ 0.05. In order to find the Lya emitters, the same 
field is then observed with one or two broad band filters with the same, or near the 
same, central wavelength as the narrow-band filter. The broad band observations 
then probe the continuum of the source and emission-line galaxies are found by 
comparing the narrow-band flux with the broad band flux. Objects with a high flux 
ratio of narrow-band vs. broad band measurement are selected. 

In the work that is presented in this thesis, the equivalent width of the line has 
been calculated in the selection process. Equivalent width (EW) is a measure of the 
strength of an emission- (or absorption-) line and is defined as 

EW = ^^ (1.2) 

J\,cont 

where Fune is the flux in the emission/absorption line and fx^cont is the flux density 
in the continuum at the central wavelength of the line. The EW is positive for 
emission lines and negative for absorption lines. Object with a fiat continuum, 
without emission or absorption lines, should have an EW of zero. The method used 
to find Lya emitters in our surveys is described in more detail in section I3.3.2[ In 
sec. II. 3[ results from past surveys for Lya emitters are reviewed. 

1.2 Lya in a historical perspective 

This thesis deals with high redshift galaxies found through their Lya emission lines. 
The Lyman series (Lya, Ly/5, Ly7. . . ) are the emission lines created when an 
electron falls from any higher energy level, back to the ground state of the simplest 
and most abundant atom in the Universe, the Hydrogen atom. The strongest line, 
with longest wavelength (1215.67 A), is the one from the second energy level to the 
first; Lya. The Lyman series are named after Theodore Lyman, a physicist who 
discovered them during the first two decades of the 20th century. 

In the mid 1920s, Menzel (1926) and Zanstra (1927) each separately discussed 
the ionisation of emission lines in planetary nebulae, including Lya, and in 1945 
Chandrasekhar discussed the radiative equilibrium of Lya. The first observations 
of highly redshifted Lya came with the discovery of quasars, or QSOs. It was in the 
1950s that star-like objects were found unexpectedly to have bright radio emission. 
In 1963 Schmidt published the first "large redshift object" (3C 273) and just two 
years later a quasar was identified to be at redshift 2.01 through its Lya emission 
(3C 9; Schmidt 1965). It was at the same time as these observations were published 
that theoretical predictions for Lya emission and absorption in the young Universe 
started. Gunn & Peterson (1965) observed that the continuum of the quasar on 
the blue side of the Lya emission line was not absorbed and put an upper limit on 
the absorption by neutral hydrogen clouds in the line of sight towards the quasar. 
Bahcall (1966) predicted the line strengths of the Lyman lines, and finally Partridge 
& Peebles (1967) published the first article with predictions on the observability 
of the Lya from young star forming galaxies, not QSOs, at high redshifts. Their 
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results have a remarkable level of agreement with the results of present day Lya 
surveys and sparked the field of narrow-band surveys for Lya emitters. 

It would, however, take a bit more than two decades before the first high red- 
shift Lya emitters were to be found. Immediately following the Partridge & Peebles 
(1967) paper, some groups attempted to discover these "primeval galaxies" (Par- 
tridge 1974; Davis & Wilkinson 1974; Meier 1976; Hogan & Rees 1979) but they were 
all unsuccessful. The narrow-band technique (see sec. I1.1.2p was brought into use 
in the late 1980s but the first surveys were unsuccessful (e.g. Pritchet & Hartwick 
1989, 1990; Rhee, Webb & Katgert 1989). The first detections came in the early 
1990s by Lowenthal et al. (1991), Wolfe et al. (1992), M0ller & Warren (1993) and 
Macchetto et al. (1993). From there on, narrow-band surveys became more and 
more frequent and successful. Examples are the LALA survey (Rhoads et al. 2003; 
Dawson et al. 2003), the Building the Bridge Survey (Fynbo et al. 2003) and the 
studies in the fields of radio galaxies (Venemans et al. 2007). To the publication 
date of this thesis, more than 550 Lya emitters have spectroscopic confirmation 
with redshifts ranging from 2; ~ 2 — 7, see Fig. II. 6[ 

1.3 What are Lya emitters? 
1.3.1 Objects that emit Lya 

A multitude of objects in the Universe emit Lja. To ionise the hydrogen atom, a 
photon with a wavelength shorter than 912 A is needed. This kind of UV radiation 
can come from several different sources. The most common sources are young, 
massive and short-lived stars of the spectral type O and B. The more massive the 
star, the shorter it lives and the bluer its spectrum is. O and B stars are the most 
massive ones known in the Universe, hence they emit almost all of their light in 
the UV range of the spectrum. They are also very luminous because they are so 
massive. Hence, Lya is often seen emitted from regions with intense star formation, 
because the short lifetime of the O and B stars make them burn out and die before 
they have time to move away from their birthplace. So, O and B stars and hence 
star forming galaxies are often intense sources of Lya. 

Another type of object that often emits Lya, or has a Lya halo around it, are 
quasars and active galactic nuclei (AGN, Schmidt 1965; McCarthy 1993; Villar- 
Martin et al. 2005). These objects are believed to consist of massive black holes 
accreting material from accretion disks surrounding them. The fact that AGNs have 
an accretion disk is in part crucial for the large Lya haloes seen around this type 
of object, as the outcoming UV light is highly collimated along the axis of the disk 
(Weidinger et al. 2005). That radio jets are observed in some AGN is also one of 
the primary arguments of why there should be an accretion disk (Urry & Padovani 
1995). The spectrum emitted by an AGN has three components, an X-ray power 
law spectrum, an infrared bump and a strong UV continuum emission component 
(Collin 2001). It is this UV component that ionises the hydrogen, which will then 
emit Lya in two highly ionised cones perpendicular to the disk (Haiman & Rees 
2001; Weidinger et al. 2005). Hence, AGN are also Lya emitters. 
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Figure 1.6: Distribution of spectroscopically confirmed Lya emitters, as of the date 
of submission of this thesis. The total number of confirmed Lya emitters is 565. 



A third possibihty for the production of Lya photons is the theory of "cold 
accretion" (e.g. Fardal et al. 2001; Dijkstra et al. 2006a,b; Nilsson et al. 2006a). 
The general idea is that if the Universe contains dark matter haloes that initially 
have no galaxy in them but plenty of neutral hydrogen gas, then the gas would 
slowly start to fall in, onto the halo because of the gravitational potential energy. 
This would cause the material to heat up because of the loss of potential, and 
this heat could theoretically be cooled off by emitting Lya photons. The idea was 
prompted after a similar effect had been seen in clusters of galaxies were the gas fell 
in, heated up and cooled via X-ray emission (Fabian 1994). On a galaxy size scale, 
this emission could be predominantly in Lya. Until recently this had never been 
observed. However, as can be read in Chapter [21 we have published the first paper 
presenting a probable observation of this phenomenon. 



1.3.2 Redshift distribution of Lya emitters 

Lya can be observed with ease in the optical regime, covering approximately the 
redshift range of 2 < 2; < 6.5. At even higher redshifts, observations are possible in 
the near-infrared. The distribution of spectroscopically confirmed Lya-emitters at 
the time of writing this thesis can be found in Fig. II. an . A much larger sample of 
candidate Lya-emitters exist, with photometric selection only. Spectroscopic follow- 
up generally has a success rate of 75 - 90% for 2; ~ 3 surveys (Fynbo et al. 2001; 
2003), and approximately 50 - 75% for 2; = 5 — 7 surveys (Taniguchi et al. 2005; 
Kashikawa et al. 2006; Shimasaku et al. 2006). 



have been trying to keep an up-to-date list of all spectroscopically confirmed Lya emitters. 
The true number may be slightly higher if I missed a publication. It will not be lower. 
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The multiply peaked distribution shown in Fig. 11.61 can be explained in the 
following way. The peak at redshift z ~ 3 is most likely due to the CCD being most 
sensitive at around 5000 A as well as to the abundance of [OIII] emission targeted 
narrow-band filters at most ground-based observatories. The peak at 2; ~ 4.5 is 
due to the LALA survey (Rhoads et al. 2000; Dawson et al. 2004). At higher 
redshifts, z > 5, CCD detector efficiency starts to drop off and the sky OH airglow 
lines gain in strength. Thus, narrow-band surveys at these redshifts are focused on 
the OH airglow "windows" where the sky background is low. Such windows exist at 
redshift z = 5.7, 6.5, 7.7, 8.8. . . Hence the two peaks in the distribution at redshifts 
z = 5.7 and 6.5. The current, undisputed, record in redshift to date is z = 6.96 from 
lye et al. (2006), although two tentative detections at z ~ 9 have been reported in 
Stark et al. (2007). Three surveys have been made ai z = 8.8 (Parkes, Collins & 
Joseph 1994; Willis & Courbin 2005; Cuby et al. 2007), with upper limits as results. 

1.3.3 Star formation rates in Lya emitters 

The star formation rate (SFR) of a galaxy denotes the amount of stellar mass 
produced in that galaxy per year. The current SFR of our galaxy, the Milky Way, is 
approximately 3 Mq yr^^ (Small 2002). When observing galaxies at high redshift, 
where stars or star forming regions are impossible to resolve, there are many SFR 
indicators in use. Almost any part of the electromagnetic spectrum can be used. 
Ranalli, Comastri & Setti (2003) argue that the X-ray flux of a galaxy can be 
used as a SFR indicator based on a relationship between the X-ray and radio/FIR 
fluxes. Radio and FIR fluxes have been shown to be tracers of star formation in 
Condon (1992) and Kennicutt (1998b) respectively. This is because radio emission 
is dominated by synchrotron emission of accelerated electrons in a magnetic plasma, 
typically from a supernovae remnant. These supernovae in turn are proportional to 
the production rate of massive, short-lived stars, which is proportional to the star 
formation rate. The FIR emission is dominated by thermal emission from the gas 
rich regions in which the stars are formed. These regions are heated by the UV light 
of the newly formed stars. The UV part of the spectrum (A ~ 1500 A) has long 
been argued to be a good tracer of the SFR of the galaxy due to the light from the 
bright, short-lived, UV intense massive stars (Kennicutt 1998a). Further, several 
emission-lines have been argued to be good tracers of star formation, especially the 
Ha emission line (Kennicutt 1983; Tresse et al. 2002) and the [Oil] emission line 
(Gallagher, Hunter & Bushouse 1989; Hashimoto et al. 1998). Assuming a certain 
recombination rate between Lya and Ha one can derive a simple relation between 
Lya luminosity and star formation rate (Kennicutt 1983; Brocklehurst 1971; Hu, 
Cowie & McMahon 1998). 

Star formation rates found in Lya emitters are typically ~ 1 — 10 Mq yr~^ (e.g. 
Cowie & Hu 1998; Hu, Cowie & McMahon 1998; Gronwall et al. 2007; Nilsson et 
al. 2007, see also references in Table [L2|) . Many authors flnd discrepancies between 
Lya SFR and the SFR derived from the UV continuum flux (e.g. Ajiki et al. 2003; 
Taniguchi et al. 2005; Gronwall et al. 2007), with larger UV derived SFRs than 
those derived from Lya. In Nilsson et al. (2007) we flnd SFRs from 0.5 to6 Mq/jt, 
but no discrepancy between continuum and emission line derived SFRs. Many 
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Table 1.2: Summary of Lya results for the cosmic star formation rate density found 
in the literature. Gronwall et al. 2007 correct for incompleteness in their sample, 
Ouchi et al. (2003) and Shimasaku et al. (2006) estimates have been calculated 
from the UV flux of the Lya emitters. Malhotra & Rhoads (2004) estimates have 
been derived from luminosity functions. 



studies of star formation rate densities have been published for Lya. Most of these 
can be found in Table 11.21 The values in this table are lower than what can be 
seen in Fig. 11.21 by almost a factor of ten. This discrepancy may be caused by dust 
extinction of the Lya/UV luminosity. Another possible explanation is that Lya 
emitters only make up 10% of the global star formation at high redshifts. Still, the 
disagreement between measurements is large, showing that the star formation rate 
density is a difficult quantity to measure. Problems with these calculations involve 
dust corrections, completeness corrections, uncertain survey volumes and not very 
well understood conversions between Lya luminosity and star formation rate. 

1.3.4 SED fitting results for Lya emitters 

Each galaxy has its own individual spectrum. Ideally, to extract all information 
about a galaxy, it is necessary to have the spectrum with high resolution. However, 
this is impractical for high redshift galaxies as each galaxy would take many hours 
to observe and only a small number of galaxies would be observed per year. The 



Reference 



Redshift psfr (Mq yr"^ Mpc"^) 



Palunas et al. 2004 
Madau et al. 1996 
Steidel et al. 1999 
Gronwall et al. 2007 
Kudritzki et al. 2000 
Nilsson et al. 2007 
Cowie & Hu 1998 
Hu et al. 1998 
van Breukelen et al. 2005 
Fujita et al. 2003 
Hu et al. 1998 
Ouchi et al. 2003 
Ajiki et al. 2003 
Rhoads et al. 2003 
Malhotra & Rhoads 2004 
Shimasaku et al. 2006 
Murayama et al. 2007 
Kodaira et al. 2003 
Malhotra & Rhoads 2004 
Taniguchi et al. 2005 



2.4 0.0024 

3.0 0.016 

3.0 0.011 

3.1 0.012 
3.13 0.0064 
3.15 0.013 
3.4 0.0047 

3.4 0.006 

3.5 0.0067 

3.7 0.00041 
4.5 0.01 

4.8 0.0063 
5.7 0.0012 
5.7 0.0005 
5.7 0.0018 
5.7 0.0023 
5.7 0.00072 
6.5 0.00052 
6.5 0.0036 
6.5 0.0013 
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alternative is to observe the galaxy in several broad-, intermediate- and/or narrow- 
band filters which in effect constitutes a spectrum with a very low resolution. The 
spectrum of the galaxy, as observed through a number of imaging filters is called the 
"spectral energy distribution" (SED). Once a well-sampled SED has been obtained 
for a galaxy it is possible to try to find the best fit between the SED and synthetic 
SEDs produced by combining populations of stars with varying ages, metallicities, 
dust properties, star formation histories, masses etc. Thus, we can gain information 
about the properties of the galaxy from the best fit synthetic SED. If the fit is made 
in a statistical way, information about the probabilities of the different properties 
can also be obtained. One of the main results in Nilsson et al. (2007) are the SED 
fits. Of course, to achieve a well constrained fit to the SED, it is necessary to have 
a good sampling, preferably from very short to very long wavelengths as different 
parts of the spectrum reveal information about different properties of the galaxy. 
For instance, the UV describes the population of massive, young stars while the 
infrared describes the old and dusty population. Thus, getting galaxy properties 
from SED fitting is constrained to galaxy samples in fields of the sky where large 
amounts of multi- wavelength data exist. 

Four papers describe SED fitting of Lya-emitters; Gawiser et al. (2006; z = 3.1, 
G06), Lai et al. (2007; z = 5.7, L07), Finkelstein et al. (2007; z = 4.5, F07) 
and Nilsson et al. (2007; z = 3.15, N07). The first and last paper examines 
the properties of Lya-emitters in GOODS-S, L07 uses data in GOODS-N and F07 
have obtained imaging of the LALA sample of Lya-emitters. All authors use the 
GALAXEV (Bruzual & Chariot 2003) synthetic spectra models with Salpeter IMF 
ranging from 0.1 - 100 Mq. G06, L07 and F07 use Calzetti et al. (1994; 1997; 2000) 
dust extinction laws, while N07 use Chariot & Fall (2000) dust models. All but N07 
find the best model by creating a grid of spectra with the different model parameters 
allowed and minimising the of each model, whereas N07 fit the best parameters 
with a Monte-Carlo Markov-Chain method (see more details in sec. 13.61) . N07 and 
G06 stack their sample and fit for the average parameters of the sample, F07 stack 
their sample in four stacks and L07 fit for each galaxy individually. The results from 
the four papers are summarised in Table 11.31 They are in good agreement with each 
other. Stellar masses of Lya emitters are not very well constrained, but lie in the 
range 10^~^ Mq. Dust contents are very low. Ages are also difficult to constrain, 
but lie in the range of a few to a few hundred Myrs. However, SED fitting is still 
a difficult and uncertain science and future work on this subject will undoubtedly 
yield new and exciting results. 

1.3.5 Large scale structure results for Lya emitters 

The luminous matter in the Universe is homogeneously distributed, on large scales, 
but on smaller scales the matter is highly clumped. Small perturbations in the initial 
mass distribution during the inflation of the Universe have collapsed into galaxy 
groups, clusters and filaments. Such large scale structure has been observed in many 
surveys (Totsuji & Kihara 1969; Peebles 1974; Bahcall & Soneira 1983; Hawkins et 
al. 2003; Eisenstein et al. 2005). The clustering of galaxies is generally measured 
by the angular correlation function, w{6), although many different suggestions have 
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Table 1.3: Summary of SED fitting results in the literature. Observed bands "Chl- 
4" refer to the Spitzcr telescope channels Chi (3.6 /um), Ch2 (4.5 /im), Ch3 (5.7//m) 
and Ch4 (8.0 fim). F07 give no result for the fitting of the dust component. In 
G06, the metallicity was set to solar metallicity, in L07 and F07 the metallicity was 
unconstrained. 



G06 L07 F07 N07 



Redshift 


3.1 


5.7 


4.5 


3.15 


Number of objects 


40 


3 


76 


23 


Observed bands 


UBVRIzJK 


BRViz 

Chl-4 


g'r'i'z' 


UBVizJHK 

Chl-4 


Stellar mass (M©) 


5 X 10^ 


1 X 109-10 


1 X 10^-9 


4.7+^^ X 10^ 


Dust 


Ay < 0.1 


E(B-V) - 0.0 - 0.5 




Ay = 0.26^°-^^ 


MetaUicity 


Z© (set) 


Unconstr. 


Unconstr. 


0.005 Zq 


Age (Myrs) 


~ 500 


5-100 


1-200 


850I420 



been made on how to calculate this value (Hamilton 1993; Landy & Szalay 1993; 
Landy, Szalay & Broadhurst 1998). The basic idea of the correlation function is 
to calculate the probability that two galaxies will be found within two infinitesimal 
solid angle elements SQi^2 separated by an angle 9. This probabihty is 



p = {i + w{e))T.^6ni6n2 (i-3) 

where S is the surface density of the galaxies. Thus, the greater the clustering at a 
certain angle, the greater w{6) is. If the sample is completely homogeneous, w{6) is 
zero. Most commonly, the correlation function is calculated by 



i?i?(0) 

In this equation, DD{Q) is the number of pairs of observed galaxies as a function 
of angle, RR{Q) is the similar number of randomly generated galaxy pairs after 
creating a large number of random fields with the same field geometry and sample 
size. DR{Q) are pairs counted when inserting an observed galaxy in each random 
frame. Correlation functions are often assumed to follow power-laws of the form 



w{e) - A^e~^ (1.5) 

with amplitude and slope p. This is valid for projected samples of galaxies in 
two dimensions. However, if the redshifts of the objects are known it is also possible 
to calculate a spatial correlation function. This will also have the form of a power 
law 
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Figure 1.7: Example of a correlation function. The correlation function is on the 
y-axis and the angular distance in arcsec on the x-axis. This function was best fit 
with an exponential form with Ay^ ~ 10 and j3 = 0.8. From Kovac et al. (2007). 



e = {r/ror (1.6) 

where tq is the so-called scale correlation length and 7 = (3 + 1. The correlation 
length gives the typical distance that a type of object is clustered on. An example 
of a correlation function can be seen in Fig. 11.71 

Clustering of Lya-emitters is a well-studied topic. Ouchi et al. (2003; 2005) and 
Shimasaku et al. (2003; 2005) present clustering results from the Subaru Deep Field 
(SDF) and Subaru-XMM Deep Field (SXDF) at redshifts z = 4.8 and 5.7. Ouchi 
et al. (2003) presents 87 candidates in SDF aX z = 4.8 and perform a correlation 
analysis on them. They find a very shallow /? of 0.1 and decide to fix it to 0.8, similar 
to that discovered for LBGs, because the Lya sample is too small to determine 
the slope. After fixing the slope, they find a correlation amplitude of = 29", 
corresponding to a correlation length of Tq = 3.5 ± 0.3h~^ Mpc. This correlation 
length becomes, after correcting for sample completeness, VQ^corr = Q.2±0.5h~^ Mpc. 
Shimasaku et al. (2003) and Ouchi et al. (2005) use a slightly different method of 
looking at overdensities in surface density within certain radii when studying samples 
oi z = 4.8 and z = 5.7 emitters in SDF and SXDF respectively. Both papers present 
overdensities with respect to field galaxies. 

Steidel et al. (1998; 2000), Hayashino et al. (2004) and Matsuda et al. (2005) 
describe the properties of a large cluster at 2; = 3.1 in the SSA22 field. Steidel et 
al. (2000) find 77 candidates which appear, by visual inspection, to be structured. 
However, the correlation analysis yielded no result. Hayashino et al. (2004) observed 
the same part of the sky, at the same redshift, but with an area ten times larger 
than Steidel et al. (2000). They find 283 candidates which again show no clear 
clustering in the correlation analysis. However, they study the surface density of 



16 



Chapter 1 



objects and find that the Lya emitter candidates are four times more clustered than 
what is predicted by the mass fluctuations of the CDM models. 

Several papers from Venemans et al. (2002; 2004; 2005; 2007) claim to have 
detected overdensities and protoclusters of Lya-emitters surrounding radio galaxies. 
In total their sample consists of ~ 300 candidates between 2 < z < 5.2 in eight fields, 
of which they have spectra for 139 confirmed Lya emitters. They claim that at least 
six out of their eight fields are overdense compared to field galaxy samples. Kovac 
et al. (2007) have studied the clustering properties of z = 4.5 Lya-emitters in 
the LALA fields. They study 151 Lya emitter candidates and find a correlation 
amplitude of = 6.73 ± 1.80 when the slope is fixed to (3 = 0.8. The amplitude 
corresponds to a scale length of tq = 3.20 ± 0A2h~^ Mpc, which after completeness 
correction becomes ro = 4.61 ± 0.60h~^ Mpc. 

Several filamentary structures have also been discovered, including one presented 
in sec. 13.5.21 M0ller & Fynbo (2001) discovered a 20 Mpc long filament with a cross 
section diameter of 1.6 Mpc at z = 3.04, consisting of 8 Lya-emitters. Matsuda et 
al. (2005) presented a filament in SSA22 consisting of 56 spectroscopically confirmed 
Lya-emitters at redshift z = 3.1. The length of the filament is 30 Mpc and its width 
is 10 Mpc. In Nilsson et al. (2007; see sec. I3.5.2p we find two projected, parallel 
filaments with a Aa significance. Follow-up spectroscopy is needed to confirm the 
filamentary structure. The distance between the filaments appear to be ~ 6 Mpc 
and the width of the filament on average ~ 1.5 Mpc. 

Finally, there are also studies which have shown Lya-emitters not to be clustered. 
In particular Shimasaku et al. (2005) and Murayama et al. (2007) find no clustering 
detections in their samples of 89 z = 5.7 emitters in SDF and 119 z = 5.7 emitters 
in COSMOS. Thus it is still unclear to what extent Lya emitters trace large scale 
structure, and if/how this depends on the fluxes of the emitters. 

1.3.6 Luminosity functions 

A luminosity function is a function describing the density of a certain type of galax- 
ies as a function of luminosity. It is always true that the brighter galaxies are more 
rare than the fainter ones. Comparing luminosity functions at different redshifts 
for the same type of object can show evolutionary changes in either density and/or 
luminosity of that type of object with time. Similarly, comparing luminosity func- 
tions of different types of objects at the same redshift will illuminate differences and 
similarities between the two types of objects. An example of a luminosity function 
can be seen in Fig. 11.81 Often, the luminosity function is showed as a cumulative 
luminosity function which means that the volume density at a certain luminosity is 
the density of all objects brighter than this luminosity. Thus, given a luminosity 
limit, it is easily possible to calculate the number of galaxies detected in a survey 
with a certain search volume. 

Luminosity functions have been shown to be well fit by a so-called "Schechter 
function" (Schechter 1976). The Schechter function has the following form: 



(j){L)dL = <i)*{L/L*Ye^^{-L/L*)dL/L* 



(1.7) 
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Figure 1.8: Example of a cumulative luminosity function. The x-axis shows the Lya 
luminosity of the Lya emitters and the y-axis the cumulative volume density. The 
red square points is the photometric sample of z = 6.5 candidates from Kashikawa 
et al. (2006) and the red triangles the spectroscopic sample from the same. The 
blue squares is the photometric sample of Shimasaku et al. (2006). The lines are 
the best fits made to the different samples. The reason that the faint end of the 
spectroscopic sample is lower than the photometric is spectroscopic incompleteness. 
From Kashikawa et al. (2006). 



where the shape of the function is determined by the three parameters a, 0* and L*. 
The a parameter determines the slope of the faint end of the luminosity function, 
L* represents the cut-off in the bright end of the function and 0* normalises the 
function. 

There are many difficulties in determining a luminosity function for Lya emit- 
ters. Firstly, a large sample of galaxies is needed in order to decrease the statistical 
error bars on the points. Ideally, the function should be made up of only spec- 
troscopically confirmed emitters, as the sample otherwise is contaminated by lower 
redshift interlopers. However, spectroscopic follow-up is offer confined to the bright 
end of the luminosity function, as the fainter an object is, the harder it is to detect 
spectroscopically. This causes the error bars to be large in the bright end due to 
low number statistics and in the faint end due to incompleteness issues. Secondly, 
when detecting Lya emitters by narrow-band imaging, several issues arise due to the 
shape of the narrow-band filter itself. For instance, the surveyed volume is smaller 
for fainter emitters, as these will only be detected at the centre of the filter and 
not at the wings. Also, some emitters may appear fainter or brighter than they 
are depending on how the photometric calibration has been done and where in the 
filter the emitter is located in redshift space. For an extensive discussion on this, 
see Gronwall et al. (2007). Because of these reasons, not many luminosity functions 
for Lya emitters have been published, and the results from the published ones differ 
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significantly. 

The results from previous work is most easily presented in a table, see Table 11.41 
In this table, the third column gives number of candidates on which the analysis is 
based, separated into spectroscopic sample and total photometric sample. Malhotra 
& Rhoads (2004) summarise several surveys, for more information on those surveys 
see that paper. The Gronwall et al. (2007) sample has spectroscopic follow-up, yet 
unpublished. The Schechter function parameter fits for the references with published 
fits can be found in Table 15.11 in Chapter [51 As can be seen in Table II. 4[ there are 
not many results on this topic yet, and many authors choose not to fit their data 
with Schechter functions, or to draw any firm conclusions. It is still unclear if any 
evolution is occurring from redshift z = 3 to z = 6.5. At lower redshifts, only one 
large survey has been presented so far (Gronwall et al. 2007). At higher redshift, 
several results have been presented but they disagree with each other to almost an 
order of magnitude in the density at the faint end of the luminosity function (e.g. 
Malhotra & Rhoads 2004; Shimasaku et al. 2006; Kashikawa et al. 2006, see also 
Fig. 15.21) . This is yet a subject where much improvement is to be expected in the 
coming decade. 



Table 1.4: Comparison of previous results for luminosity functions. For more information about this table, see text. 



Reference 


Redshift 


Sample size 
(bpec./rnot.) 


Area (arcmin^) 


Field 


Conclusions 


(Jucm et al. (2UUc5j 


4.86 


— /87 


543 




AT 1 4- n A 

i\lo evolution between z = 3.4 

J A O T~) * 1 J. 11 £ T 

and 4.8. Bright end slope oi Lya 
continuum luminosity function 

steeper than that of LEG LFs. 


Maier et al. (2003) 


4.8, 5.7 


2/- 


100 


CADIS 


Fabry-Perot detections. Small 
number of detections, discuss 
large scale structure influence 
on results. Find no change in LF 
between z = 3.5 and 5.7. 


tiu et al. (2004) 


5.7 


18/26 


918 


CO A oo 

SSA22 


AT 1 * J 

iNo conclusions made 


Malnotra & 


5.7, 6.5 








/- 1 • 1 • 1 1 • 1 1 

Compiles various published 


Rhoads (2004) 










surveys. Find no real evolution 
between z — 4.5 and 5.7. 

Conclude that the Universe is 
ionised at 2; ~ 6. Publish fits 
to a Schechter function. 


van Breukelen 


2.5 - 4.4 


14/- 


1.36 




IFU observations. Publish fit to 


et al. (2005) 










Schechter function. Find no 
evolution between z — 3.4 and 5.7. 


Tapken et al. 


5.7 


8/15 


46 


FORS 


Find more Lya emitters than 


al. (2006) 








Deep Field 


expected compared to 
Malhotra & Rhoads (2004) LF. 



Table 1.5: Table [L^ continued. 



Shimasaku et 
al. (2006) 


5.7 


28/89 


725 


SDF 


Publish fit for Schechter function. Find 
no change between z = 4.5 and 5.7. 


Kasliikawa et 
al. (2006) 


6.5 


17/58 


876 


SDF 


Similar selection criteria as Shimasaku 
et al. (2006). Publish Schechter function fit. Large 
spectroscopic follow-up rate. Find deficit in bright end 
of LF by a factor of 2 compared to 2; = 5.7 but cannot 
rule out cosmic variance. 


Gronwall et 
al. (2007) 


3.1 


— /162 


1008 


E-CDFS 


Publish Schechter function fit. 
Concludes that if Malhotra & Rhoads (2004) are 
correct then emitters at z = 3.1 are 2.5 times brighter 
or more numerous than at 2; = 5.7 but if Shimasaku et al. 
(2006) are correct, then the opposite result. 
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1.3.7 Modelling of Lyo; emitters 

Modelling Lya is a difficult business. Lya is a resonance line that can be affected 
by many occurrences, such as e.g. dust, in- and outflows, geometry, clumping etc. 
Many groups have made various models about the Lya emission line or Lya emitting 
galaxies. Some are summarised here. 

Laursen & Sommer-Larsen (2007) present a new three dimensional Lya radiative 
transfer code based on smoothed particle hydrodynamics simulations. The code is 
a Monte-Carlo code which propagates the Lya photons through the medium sur- 
rounding the galaxy. The code is shown to be successful in describing one of the 
observed properties of Lya emitters; that the Lya emission is often more extended 
than the continuum emission (M0ller & Warren 1998; Fynbo et al. 2001; 2003). 
A similar result is also achieved by Purlanetto et al. (2003; 2005), who also calcu- 
late the radiative transfer through a smoothed particle hydrodynamics simulation. 
Purlanetto et al. (2005) also present luminosity functions for Lya emitters at z — 3. 
Cantalupo et al. (2005) use the same approach to study the luminosity of fluores- 
cent Lya sources at high redshift. Verhamme, Schaerer & Maselli (2006) present 
another radiative transfer code with which they have studied the line profile of Lya. 
Yet another description of the same method can be found in Tasitsiomi (2006). The 
conclusion from this paper is that a 2; ~ 8 Lya galaxy would not be observable with 
current, ground-based observatories. 

Several papers by Dijkstra ct al. (2006a; 2006b; 2007a; 2007b; 2007c) also use 
a radiative transfer code for Lya photons, applied to a cosmological simulation. In 
the two papers from 2006, the authors describe the effects of infall of material onto 
a galaxy on the Lya emission. Dijkstra et al. (2007a) study the effect of the IGM 
on the emergent Lya and discuss the possibility to use these galaxies to determine 
re-ionisation at very high redshifts. The latter issue is also discussed in Dijkstra et 
al. (2007b), in relation to recent new results regarding the luminosity function of 
Lya emitters at very high redshift. Finally, Dijkstra et al. (2007c) indicate that 
up to 50% of all Lya emitters may harbour primeval stars (so-called Population III 
stars). 

After the discovery of two large Lya blobs by Steidel et al. (2000), several 
explanations for this phenomenon were presented in different publications. The 
explanations proposed, all based on trying to find a way of producing the energy 
needed to ionise the nebula, were i) star formation, where the energy comes from 
newly formed stars, ii) AGN activity, where the energy comes from the accretion 
of material onto a supermassive black hole, and in) cold accretion, where the en- 
ergy comes from loss of gravitational potential energy of cold gas accreting onto a 
dark matter halo. The first explanation, star formation, is based on the idea that 
a starburst galaxy will produce a large amount of young, massive stars emitting 
much UV, and thus Lya, radiation. These massive stars will rapidly burn out and 
become supernovae, and the explosions from these supernovae will create a sort of 
"superwind" which will blow out the ionised material surrounding the galaxy. Thus 
the galaxy will be enshrouded in a Lya emitting cloud. The modelling behind this 
scenario is described in Taniguchi et al. (2001), Ohyama et al. (2003), Mori et 
al. (2004) and Wilman et al. (2005). The second explanation, involving an AGN, 
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involves that the AGN can be obscured from the hne-of-sight, but still ionise the 
material surrounding it. This method is explained in Haiman & Rees (2001) and 
Weidinger et al. (2004; 2005). Finally, cold accretion, as described in e.g. Haiman, 
Spaans & Quataert (2000), Fardal et al. (2001), Dijkstra et al. (2006a,b) and Dekel 
& Birnboim (2006), is described as infalling gas onto a dark matter halo. This in- 
falling gas is heated due to the loss of gravitational energy. The heat can then be 
released through Lya emission. This topic is further discussed in Chapter [2] of this 
thesis. 

Two papers study the theoretical aspects of filaments discovered in Lya surveys; 
Weidinger et al. (2002) and Monaco et al. (2005). In Weidinger et al. (2002) it 
is proposed that filaments found with relative ease through narrow-band imaging 
can be used as an independent measure on cosmological parameters such as Qm and 
Q\. The idea is that, assuming that all filaments must be isotropically oriented, 
filaments will be stretched differently depending on the assumed cosmology of the 
Universe. Thus, a large sample where the individual alignments become washed out 
could give an independent measure of the cosmological parameters. In Monaco et 
al. (2005), the effect of peculiar motions on the redshifts in observed filaments is 
quantified. It is shown that the effect is small, but non-negligible and should be 
taken into account when drawing conclusions from filamentary structures. The idea 
of determining cosmological parameters with Lya emitters has been picked up in 
the HETDE}ql project, described in Hill et al. (2004). HETDEX proposes to use 
a new integral field unit (IFU) spectrograph called VIRUS (Visible IFU Replicable 
Ultra-cheap Spectrograph) to study several thousands of Lya emitters in the redshift 
range 1.8 < z < 3.8 every night on the Hobby-Eberly Telescope, thus being able 
to study large scale structure of this type of galaxy to extreme accuracy, and from 
that determine the cosmological parameters to a great level of certainty. 

Barton et al. (2004), Le Delliou et al. (2005; 2006), Thommes & Meisenheimer 
(2005), Stark, Loeb & Ellis (2007) and Kobayashi et al. (2007) make predictions for 
very high redshift Lya surveys. Barton et al. (2004) use cosmological hydrodynamic 
simulations of galaxy formation to study the detectability of z > 7 sources in the near 
future. Le Delliou et al. (2005; 2006) use a semi-analytic model based on the ACDM 
model and show that it fits well with lower redshift observations. Kobayashi et al. 
(2007) use a similar semi-analytical model but with a different parametrisation of 
the Lya escape fraction which they claim to be more physical. They also claim good 
fits with lower redshift results. Thommes & Meisenheimer (2005) and Stark, Loeb 
& Ellis (2007) instead use phenomenological models, assuming that Lya luminosity 
is proportional to the star formation rate of a galaxy, and that the star formation 
rate is proportional to the baryonic mass of the galaxy. The results from these 
papers are luminosity functions at several redshifts between z = 3 — 8. This work is 
followed up in Chapter [5l where we make predictions for high redshift narrow-band 
surveys. 

In preparation for the many surveys for very high redshift Lya emitters coming 
in the next decade, several groups discuss how to draw conclusions about the re- 
ionisation history of the Universe. Loeb & Rybicki (1999) calculate the brightness 
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distribution and tlie sliape of the Lya emission line in haloes surrounding high- 
redshift galaxies before re-ionisation was complete. They find that observations 
of Lya prior to re-ionisation will be a good probe of the neutral IGM. Haiman & 
Spaans (1999) look at dust effects on high redshift Lja. They discuss how changes 
in the Lya luminosity function may be used to determine the level of re-ionisation. 
Many authors have quantified the effect of the Gunn-Peterson absorption (Gunn & 
Peterson 1965) on the red wing of the Lya line as a function of redshift (Miralda- 
Escude 1998; Miralda-Escude & Rees 1998; Haiman & Loeb 1999; Haiman 2002; 
Santos 2004; Gnedin & Prada 2004; Haiman & Gen 2005; Tasitsiomi 2006). This 
effect comes from the vast amount of neutral hydrogen in the Universe at high 
redshift. Photons with wavelengths less than that of Lya at its place of origin 
will travel towards us, but is redshifted in the process. Depending on the original 
wavelength, it will at a certain redshift have the wavelength of Lya. At high redshift, 
the likelihood to encounter a neutral hydrogen cloud at this distance from the source 
is high, resulting in that almost all the photons blue-ward of Lya from high redshift 
galaxies are absorbed. At very high redshift, even photons in the red wing of the 
Lya line will be absorbed and thus the level of absorption, i.e. the asymmetry of 
the line, will reveal information of the level of re-ionised material surrounding the 
galaxy. 

McQuinn et al. (2007) show, using a radiative transfer simulation, that obser- 
vations of the clustering of very high redshift Lya emitters may reveal the level of 
re-ionisation independently of the luminosity function and the line profile. This will 
be a good method to use to cross-check with the results of other methods. 

In a very recent paper, Fernandez & Komatsu (2007) use the mass-to-light ratio 
of Lya galaxies to make luminosity functions. The general idea of their paper is 
to use a derived mass function, and a mass-to-light ratio which is a free parameter 
to calculate a luminosity function. From comparisons with observed high redshift 
luminosity functions, they then argue that Lya emitters are either starburst galaxies 
with low escape fractions, or normal galaxies with a higher escape fraction. The also 
see a hint of metallicity evolution between redshift z = 5.7 and 6.5, but no evidence 
that re-ionisation should have ended by z ~ 7. 

1.4 This thesis 

This thesis presents work I have done regarding Lya emitters. The different chapters 
describes various aspects of these galaxies. Both theoretical and observational, as 
well as low and high redshift work is presented. The main part of the thesis is 
divided into two parts. The first part concerns narrow-band imaging for medium 
redshift {z ~ 3.15) Lya emitters in the survey field of GOODS-S. The outcome of 
that project was two papers (Nilsson et al. 2006a; Nilsson et al. 2007), presented 
in Ghapter [2] and Ghapter [3] respectively. In Ghapter [2] we described the discovery 
of a Lya blob, which after careful analysis turned out to be the first of its kind 
best explained by cold accretion onto a dark matter halo. In Ghapter [3] all other 
results from the original narrow-band imaging are presented, including the imaging 
and spectroscopic results we have, the photometry of the candidates and spectra 
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of three confirmed candidates, a tentative discovery of a filamentary structure, a 
careful SED fitting procedure and a comparison with Lyman Break Galaxies. 

The second major project I have been heavily involved in from the start is a 
contribution to the Visible and Infrared Survey Telescope for Astronomy (VISTA). 
My original idea was to buy a set of narrow-band filters for VISTA, to detect a 
sample of z = 8.8 Lya emitters. VISTA is a new survey telescope in the near- 
infrared, with a camera with a field-of-view of 0.6 deg^. It will be dedicated to Public 
Surveys, enabling large, time- demanding projects to be carried out. I have been 
involved in all aspects of this project, from design and procurement of the filters, to 
writing a Public Survey proposal and later a Survey Management Plan. This work, 
described in Chapter HI prompted a new project, presented in Chapter [5l In this 
chapter we attempt to use two different theoretical models to calculate luminosity 
functions for Lja emitters at very high redshift {z > 7). We use these results to 
predict numbers of galaxies detected in a number of present and future surveys. 

In the final science chapter we explore the colour space that Lya emitters inhabit 
compared to that of "normal" galaxies or [OII]-emitters. We are writing a paper 
with the aim to try to develop a method to improve the narrow-band photometric 
selection technique using optical and near-infrared colours. The results from this 
project are found in Chapter [61 

In Chapter [7] I summarise the work presented in this thesis and in the final 
chapter. Chapter [HI I attempt to look forward at what may happen in this field of 
study in the near future, and what my contribution could be. 



Chapter 2 

A Lya blob in GOODS-S 



This paper has been pubhshed in Astronomy & Astrophysics Letters in June 2006 
(A&A, 452, L23). The authors are Nilsson, K.K., Fynbo, J.P.U., M0ller, P., 
Sommer-Larsen, J., & Ledoux, C. 

2.1 Abstract 

We report on the discovery of a z — 3.16 Lyman-a emitting blob in the Great 
Observatories Origins Deep Survey (GOODS) South field. The discovery was made 
with the VLT, through narrow-band imaging. The blob has a total Lya luminosity 
of ~ 10^^ erg s^^ and a diameter larger than 60 kpc. The available multi-wavelength 
data in the GOODS field consists of 13 bands from X-rays (Chandra) to infrared 
(Spitzer). Unlike other known Lya blobs, this blob shows no obvious continuum 
counter-parts in any of the broad-bands. In particular, no optical counter-parts are 
found in deep HST/ACS imaging. For previously published blobs, AGN (Active 
Galactic Nuclei) or "superwind" models have been found to provide the best match 
to the data. We here argue that the most probable origin of the extended Lya 
emission from this blob is cold accretion onto a dark matter halo. 

2.2 Introduction 

Narrow-band surveys for Lyman-a (Lya) emitting galaxies at high redshift have 
recently revealed a number of luminous (up to 5 ■ 10^'^ erg s~^), very extended (from 
a few times ten kpc to more than 150 kpc) Lya-emitting objects, so-called Lya 
"blobs" (Fynbo et al. 1999; Keel et al. 1999; Steidel et al. 2000; Francis et al. 
2001; Matsuda et al. 2004; Palunas et al. 2004; Dey et al. 2005; Villar-Martin 
et al. 2005). At least three mechanisms have been suggested as energy sources for 
Lya blobs. These are: i) hidden QSOs (Haiman & Rees 2001; Weidinger et al. 
2004, 2005), a) star formation and superwinds from (possibly obscured) starburst 
galaxies (Taniguchi et al. 2001; Ohyama et al. 2003; Mori et al. 2004; Wilman 
et al. 2005), and in) so-called cold accretion (Haiman, Spaans & Quataert 2000; 
Fardal et al. 2001; Keres et al. 2004; Mailer & Bullock 2004; Birnboim & Dekel 
2003; Sommer-Larsen 2005; Dijkstra et al. 2006(a,b); Dekel & Birnboim 2006). 

0^ 
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Cooling flows are phenomena observed in galaxy clusters for more than a decade 
(Fabian 1994). These are explained by gas which is cooling much faster than the 
Hubble time through X-ray emission in the centres of the clusters. However, cooling 
emission from a galaxy, or a group sized halo can be dominated by Lya emission 
(e.g. Haiman, Spaans & Quataert 2000; Dijkstra et al. 2006(a,b)). In this Letter 
we present the discovery of a Lya blob at redshift 2; ^ 3.16 located in the GOODS 
South field, which we argue is the first piece of evidence for cold gas accretion onto 
a dark matter halo. 

Throughout this paper, we assume a cosmology with Hq = 72 km s~^ Mpc~^, 
flra = 0.3 and fl\ = 0.7. All magnitudes are in the AB system. 



2.3 Observations and Data reduction 

A 400x400 arcsec^ section, centred on R.A. = 03'^32"^21.8^ Dec = -27°45'52" 
(J2000), of the GOODS South field was observed with FORSl on the VLT 8.2 m 
telescope Antu during two visitor mode nights on December 1-3, 2002. A total of 16 
dithered exposures were obtained over the two nights for a combined exposure time 
of 30 ksec, all with the narrow band filter OIII/3000+51 and using the standard 
resolution collimator (0.2x0.2 arcsec^ pixels). For this setup the central wavelength 
of the filter is 5055 A with a FWHM of 59 A, corresponding to the redshift range 
z = 3.126 - 3.174 for Lya. 

The observing conditions were unstable during the two nights with the seeing 
FWHM varying between 0.66" and 1.25" on the first night and 1.4" and 3.3" on the 
second night. The images were reduced (de-biased, and corrected for CCD pixel- 
to-pixel variations using twilight flats) using standard techniques. The individual 
reduced images were combined using a modified version of our code that optimizes 
the Signal-to-Noise (S/N) ratio for faint, sky-dominated sources (see M0ller & War- 
ren 1993, for details on this code). The modification of the code was necessitated by 
the highly variable seeing. The sky background was assumed to be constant. The 
FWHM of the PSF of the final combined narrow-band image is 0.8". 

For object detection, we used the software package SExtractor (Bertin & Arnouts 
1996). A full description of our selection of Lya emitters in the GOODS field will be 
given in a subsequent paper. In this Letter we discuss the nature of an extended, low 
surface brightness blob with a centroid (of the Lya emission) of R.A. = 03'^32™14.6'* 
and Dec = — 27°43'02.4" (J2000) detected in the combined narrow-band image. 

Follow-up MOS spectroscopy was obtained in service mode using FORSl/VLT 
UT2 over the time period December 2004 - February 2005. The total observing 
time was 6 hours. We used a 1.4" slitlet and grism 600V resulting in a wavelength 
range of 4650 A to 7100 A and a spectral resolution FWHM of approximately 700. 
The seeing varied between 0.77" and 1.2" during the spectroscopic observations. 

The GOODS archival data used here and their detection limits are listed in 
Table O 
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Table 2.1: Specifications of deep, multi-wavelengtli data available in the GOODS 
South field and the narrow-band image. The last column gives the 3a limit as 
detected in a 2" radius aperture and the narrow-band value gives the blob flux in 
this aperture. 



Filter/Channel 


Ac 


Filter 3a limit (2" aperture) 






FWHM (erg cm'^ s^^ Hz-^) 



X-rays {Chandra) 


4.15 keV 


3.85 keV 


9.90 


10" 


-34 


U {ESQ 2.2-m) 


3630 A 


760 A 


8.62 


10" 


-31 


B (HST) 


4297 A 


1038 A 


9.25 


10" 


-30 


Narrow ( VLT) 


5055 A 


60 A 


6.68 


10- 


-30 


V {HST) 


5907 A 


2342 A 


4.66 


10- 


-30 


i (HST) 


7764 A 


1528 A 


1.50 


10- 


-29 


z (HST) 


9445 A 


1230 A 


3.00 


10- 


-29 


J ( VLT) 


1.25 /im 


0.6 /im 


5.31 


10- 


-30 


H ( VLT) 


1.65 /im 


0.6 /im 


1.86 


10- 


-29 


Ks ( VLT) 


2.16 /im 


0.6 /im 


1.56 


10- 


-29 


Chi (Spitzer/IRAC) 


3.58 /im 


0.75 /im 


2.51 


10- 


-31 


Ch2 (Spztzer/IRAC) 


4.50 /im 


1.02 /im 


6.43 


10- 


-32 


Ch3 {Spitzer/IRAC) 


5.80 /im 


1.43 /im 


5.01 


10- 


-29 


Ch4 {Spitzer/IRAC) 


8.00 /im 


2.91 /im 


4.65 


10- 


-30 



2.4 Results 

The spectrum of the part of the Lya blob covered by the slitlet can be seen in the 
left-most panel of Fig. 12.11 The line has the asymmetric profile expected for a high 
redshift Lya emitter. We detect no other emission lines in the spectrum. The most 
likely interloper is [Oil] at redshift 0.36, but no emission is observed in the spectrum 
where e.g. H/9 or [OIII] are expected at this redshift, see Fig. l2.1[ This leads us to the 
conclusion that we are observing a Lya-emitting object aX z = 3.157. The observed 
FWHM velocity width of the emission line is 505 km s~^. The instrument FWHM 
of the set-up is 290 km s~^, hence the Lya intrinsic velocity width is marginally 
resolved. The intrinsic width is less than 500 km s~^. This is of the order or smaller 
than for other published blobs, with velocity widths of 500 — 2000 km s~^ (Keel et 
al. 1999; Steidel et al. 2000; Francis et al. 2001; Ohyama et al. 2003; Bower et al 
2004; Dey et al. 2005). 

A contour-plot of the blob superimposed on the HST/ ACS V-band image is 
shown in the middle panel of Fig. 12.11 and a plot of the surface brightness of the 
blob is seen in the right panel of the same figure. The full set of thumb-nail images of 
the blob in all 14 bands can be found in Fig. 12. 2[ No obvious continuum counterpart 
is detected in any band. The radial size is at least 30 kpc (60 kpc diameter) with 
fainter emission extending out to 40 kpc radius. This can be seen as the extension 
to the SW in the contour-plot in Fig. 12.11 The significance of the lowest contour 
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Figure 2.1: Left a) Flux calibrated spectrum of the blob emission line. The line 
has the characteristic blue side absorption, indicating high redshift. b) The part 
of the spectrum (binned with a binsize equal to half the resolution (1.1 A)) where 
H/3 and [OIII] should have been observed if the emission line was [Oil] at a redshift 
of 2; ~ 0.36. These lines are not observed and therefore we conclude the observed 
line is due to Lya at z = 3.16. Middle Contour-plot of narrow-band emission from 
the Lya blob overlaid the HST V-band image. The narrow-band image has been 
continuum subtracted by subtracting the re-binned, smoothed and scaled HST/V- 
band image. Contour levels are 2-10"'^, 4 -10"^ and 6- 10~^ erg s~^ cm~^ in restframe 
flux (corresponding to 1.2 ■ 10~^^, 2.5 ■ 10~^^ and 3.7 ■ 10~^^ in observed flux). The 
image is 18" x 18" (18" corresponds to a physical size of ~ 133 kpc). Numbers refer 
to those used in section 12. 4[ The dotted hues indicate the slitlet position for our 
follow-up spectroscopy. Right Plot of surface brightness as function of radius. The 
flux is the sky subtracted narrow-band flux. The PSF of the image is illustrated by 
the solid line, and the dotted line is the best fit model of Dijkstra et al. 2006. The 
deficit at ~ 45 kpc is due to the asymmetric appearance of the blob. 



levels is of the order of 2a per pixel. The total Lya luminosity, in a 30 kpc radius 
aperture, is LLya = (1.09 ± 0.07) ■ 10'^'^ erg s"^ This coincides, after correction for 
the smaller area sampled in the spectrum, to the Lya flux detected in the spectrum 
within errors. A conservative lower limit to the restframe equivalent width (EW) 
of the emission line can be calculated from upper limits on the broad-band fluxes 
in the HST B and V filters in the same aperture. This limit is EW > 220 A in the 
restframe. This is in the range of previously published Lya blobs, that have a Lya 
flux to B-band flux density range between 50 - 1500 A in the restframe (but typically 
these values are derived measuring the continuum flux in a smaller aperture than 
the emission line flux). 

There are seven objects detected in a wide range (> 8) of energy bands, within a 
10" radius surrounding the blob. 8 other objects are detected within the V-band and 
one further detected in the Spitzer/IRAC bands. The photometric redshift of these 
objects was calculated using the public Hyper^ code by Bolzonella et al. (2000). 
The resulting photometric redshifts for the eight objects with most data points (> 8) 
can be found in Table 12.21 The other eight objects detected in the V-band have 



"'^ http://webast.ast.obs-niip.fr/hypcrz/ 
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Figure 2.2: Thumbnail images of all available multi-wavelength data in the GOODS 
South field, centred on the Lya blob. All images are 18" x 18". 



Table 2.2: Photometric redshifts of objects surrounding the blob. Numbering refers 
to those given in Fig. 12.11 Errors given are 1 a. 



Obj # 


Dist. from blob 

(arcsec) 


Zphot 


x'/d.o.f 


Type 


Ay rest 


1 


4.6 


1 -I+U.4U 
-'-•-'--0.30 


1.3 


Burst 


0.20 


2 


4.6 


1 1+0.34 
-'-•-'--0.41 


8.6 


Burst 


1.20 


3 


6.8 


^•^-0.59 


4.8 


Spiral 


1.20 


4 
5 


8.4 
8.7 


"-'•"-0.63 
'-'•^-0.89 


8.1 
2.0 


Burst 
Burst 


0.40 
0.20 


6 


3.0 


^•^-1.54 


0.9 


Spiral 


1.20 


7 


6.3 


-'-•-'--0.81 


1.9 


Burst 


1.20 


8 


4.5 


q p:+1.27 
'-'•"J-3.48 


0.6 


Spiral 


0.00 



only a few detections across the spectrum and hence their photometric redshifts are 
unreliable. The redshift of object # 3 is similar to the blob redshift and indicates 
that this galaxy may be near to the blob. Object # 6 is an intriguing object, 
undetected in the deep optical and near-IR imaging but bright in the Spitzer/IRAC 
bands. Its photometric redshift is consistent with the redshift of the blob, but with 
a large uncertainty. The object is also detected in the Spitzer/MIPS 24 fim band. 
Based on the Spitzer magnitudes, and on the diagnostic colour-colour diagram of 
Ivison et al. (2004), object # 6 is best fit by a star-burst at high redshift {z ~ 5.5, 
consistent with the photometric redshift estimate), hence unrelated to the blob. 



2.5 Discussion 

We first consider that the Lya emission of the blob may be due to recombination of 
gas, photo-ionized by an AGN or a starburst galaxy. If the blob is a "passive" gas 
cloud illuminated and photo-ionized by a nearby AGN, then, following Cantalupo 
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et al. (2005), one can show that for an AGN with luminosity L,y = Lll('^/'^ll)~", 
and in order to result in a peak blob Lja surface brightness of T^Lya, the AGN has 
to be located at a distance of no more than 



where equality applies to the case where the blob gas is optically thick at the Lyman 
limit. No AGN has been detected in the deep Chandra image available within this 
distance. Worslcy ct al. (2005) argue that a significant fraction of the unresolved X- 
ray background in hard bands consists of highly obscured AGN. However, Worsley 
et al. (2005) also predict that the AGN responsible for this background are situated 
at z ~ 0.5 — 1.0. Furthermore, Silverman et al. (2005) present a luminosity function 
for AGN at higher redshift. To the detection hmit of the CDFS {Lx{z = 3.15) ^ 
1.9 • 10^3 s~^) and with our search volume (3x3 Mpc x = 0.05) we expect to 
detect only 0.06 AGN in our entire search volume. We also consider the possibility 
that galaxy # 3 can photo-ionise the blob. However, if we assume a power law for 
the spectrum and extrapolating from the HST/B and HST/V detections we find 
that the UV luminosity of galaxy # 3 is not sufficient to photo-ionise the blob, 
unless highly coUimated towards the blob. We have no reason to believe that this 
is the case. 

The second possibility is that the blob Lya emission is somehow related to 
starburst driven, superwind outflows. A starburst would be expected to be located 
within the blob to create such a Lya halo and no central continuum source has 
been detected. Even though a very massive starburst can be made invisible in 
the UV/optical range by dust obscuration, it should be visible in the IR, i.e. the 
Spitzer/IRAC bands. 

The third option is that the Lya emission is due to cold accretion of predomi- 
nantly neutral, filamentary gas onto a massive dark matter halo. For cold accretion, 
the bulk of the Lya emission is produced by collisional excitation, rather than re- 
combination. Recently, Dijkstra et al. 2006 (a,b) presented a theoretical model for 
Lya cooling fiows, along with predictions of the emission line profile and the shape 
of the surface brightness function. The S/N of our spectrum is not high enough to 
allow a comparison of emission line profiles. However, the surface brightness profile 
matches well the predictions for a centrally illuminated , collapsing cloud of Dijk- 
stra et al. 2006(a), see Fig. 1. Further tests are needed to determine how well their 
model fits. To test whether this blob can be filamentary gas accreting "cold" onto 
a companion galaxy, we also conducted the following experiment: we calculated the 
Lya surface brightness in a 100x100 kpc (projected) region for a proto-galaxy of 
"cooling" radiation only (so all contributions from regions with young stars were 
removed, as well as all emission, in general, from gas closer than 10 kpc to any 
star-forming region). The calculation was based on a cosmological simulation of the 
formation and evolution of an M31-like disk galaxy (Sommer-Larsen 2005; Portinari 
& Sommer-Larsen 2005). 

The results at z ~ 3 are presented in Sommer-Larsen (2005), and get to a 
surface brightness about an order of magnitude lower than the observed level. This 
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is interesting, and may point to a cold accretion origin of the blob Lya emission on 
a larger scale, such as filamentary gas accretion onto a galaxy-group sized halo. 
Another possibility is that the periods with high surface brightness are shorter 
than 2.5 Myr (the resolution of the simulation). Given that in a search volume 
of about 40000 co- moving Mpc^, only one such blob has been detected, it is actually 
comforting, that we could not reproduce the blob characteristics, by cold accretion 
onto this, randomly selected, M31-like galaxy. This has to be a rare phenomenon. 

A test for the cold accretion model would be to observe the Balmer lines. For 
coUisionally excited hydrogen, neglecting extinction effects, the flux in Ha should 
only be about 3.5 percent of the Lya flux, whereas for recombining, photo-ionized 
gas this ratio is ~ 11.5 % (Brocklehurst 1971). Hence, the relative Ha luminosity 
is expected to be significantly larger in the latter case. The situation is similar for 
H/5, and whereas the Ha line will be very difficult to detect from the ground, H/5 
should be observable. 

2.6 Conclusion 

We have here reported the results of an extensive multi-wavelength investigation 
of a redshift z = 3.16 Lya emitting blob discovered in the GOODS South field. 
The blob has a diameter larger than 60 kpc diameter and a total luminosity of 
Llyq ~ 10^^ erg s~^. Deep HST imaging show no obvious optical counterpart, 
and the lack of X-ray or IR emission suggest there are no AGN or dusty starburst 
components associated with at least the centroid of the blob. Two galaxies within 
a 10" radius have photometric rcdshifts consistent with the redshift of the blob, but 
follow-up spectroscopy is needed to establish if there is a connection. We have run 
simulations of Lya surface brightness arising from cold accretion and found that such 
extended Lya emission may be explained by accretion along a filament onto a galaxy 
group sized dark matter halo. Another possibility is that such emission in very 
short lived, i.e. significantly shorter than the 2.5 Myr resolution of our simulation. 
We argue that other previously suggested origins of Lya blobs (hidden AGN and 
"super- winds" ) can be ruled out in this case due to the lack of detected continuum 
counter-parts. Hence, though our cold accretion simulation cannot perfectly match 
our data, it is the only explanation that is plausible. Our results combined with the 
fact that previously studied blobs appear to be caused by superwinds and/or AGN 
in turn implies that the energy sources for blob Lya emission are diverse. 
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Chapter 3 

Lya emitters in the GOODS-S 
field 

This paper has been pubhshed in Astronomy & Astrophysics in August 2007 (A&A, 
471, 71). The authors are Nilsson, K.K., M0ller, P., MoUer, O., Fynbo, J.P.U., 
Michalowski, M.J., Watson, D., Ledoux, C, Rosati, P., Pedersen, K., & Grove, L.F. 

3.1 Abstract 

Conteoct Lya-emitters have proven to be excellent probes of faint, star-forming galax- 
ies in the high redshift universe. However, although the sample of known emitters 
is increasingly growing, their nature (e.g. stellar masses, ages, metallicities, star- 
formation rates) is still poorly constrained. 

Aims We aim to study the nature of Lya-emitters, to find the properties of a typical 
Lya-emitting galaxy and to compare these properties with the properties of other 
galaxies at similar redshift, in particular Ly man-break galaxies. 

Methods We have performed narrow-band imaging at the VLT, focused on Lya 
at redshift z ^ 3.15, in the GOODS-S field. We have identified a sample of 
Lya-emitting candidates, and we have studied their Spectral Energy Distributions 
(SEDs). 

Results We find that the emitters are best fit by an SED with low metallicity 
{Z/Zq = 0.005) , low dust extinction (Ay ~ 0.32) and medium stellar masses 
of approximately 10^ Mq. The age is not very well constrained. One object out of 
24 appears to be a high redshift Lya-emitting dusty starburst galaxy. We find fila- 
mentary structure as traced by the Lya-emitters at the 4o" level. The rest-frame UV 
SED of these galaxies is very similar to that of Lyman Break Galaxies (LBGs) and 
comply with the selection criteria for [/-band drop-outs, except they are intrinsically 
fainter than the current limit for LBGs. 

Conclusion Lyct-emitters are excellent probes of galaxies in the distant universe, 
and represent a class of star-forming, dust and AGN free, medium mass objects. 
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3.2 Introduction 

The possibility to use the Lya emission hne to study galaxies in early stages of 
their formation was outlined already by Partridge & Peebles (1967) nearly 40 years 
ago, but early surveys (see Pritchet 1994 for a review) failed to produce anything 
other than upper limits. The unexpected faintness of the objects caused it to take 
almost three decades before the narrow-band technique was successfully used to 
identify the first high redshift Lya emitting galaxies that were not dominated by 
Active Galactic Nuclei (e.g., Lowenthal et al. 1991; M0ller & Warren 1993; Hu & 
McMahon 1996; Petitjean et al. 1996; Francis et al. 1996; Cowic & Hu 1998). It is 
only recently, with the advent of 8 m class telescopes and sensitive detectors, that 
larger samples of Lya selected objects have been reported (e.g. Steidel et al. 2000; 
Malhotra & Rhoads 2002; Fynbo et al. 2003; Ouchi et al. 2003; Hayashino et al. 
2004; Venemans et al. 2005). 

Already during the early studies, two interesting suggestions were raised. First, 
it was found that there is a tendency for Lya selected objects to "line up" as strings 
in redshift space, and that they therefore may be excellent tracers of filaments at 
high redshifts (Warren & M0ller 1996; Ouchi et al. 2004a,b; Matsuda et al. 2005). 
Secondly, they were found to have very faint broad band magnitudes and therefore 
could be good tools in detecting faint, high redshift galaxies (Fynbo et al. 2001; 
Fujita et al. 2003; Venemans et al. 2005; Gawiser et al. 2006). 

The first of those suggestions was explored theoretically via modeling of struc- 
ture formation including assignment of Lya emission to the models (Furlanctto et 
al. 2003; Monaco et al. 2005), and has been confirmed observationally (M0ller & 
Fynbo 2001; Hayashino et al. 2004), who also proposed to use such structures for 
a new cosmological test to measure Q/^ by looking at the "length-to-radius" ratio 
of filaments observed from the side or end-on. This proposed test was subsequently 
explored in detail by Weidinger et al. (2002). 

The second suggestion has gained significant interest because there are now sev- 
eral additional, but independent, ways of identifying high redshift but optically faint 
galaxies, e.g. Damped Lya Absorbers (DLA) galaxies (Wolfe et al. 1986; M0ller et 
al. 2002; Wolfe, Gawiser & Prochaska 2005), Gamma- Ray Burst (GRB) host galax- 
ies (Fruchter et al. 2006), sub-mm galaxies (Chapman et al. 2004). Most galaxies 
from such searches are however too faint to be identified in current ground-based 
optical fiux limited samples. A significant project (the Building the Bridge Survey, 
BBS) aimed at addressing those issues is currently underway at the ESQ Very Large 
Telescope (VLT) (Fynbo et al. 2001; Fynbo et al. 2003). The very faintness of the 
objects, however, renders it difficult to make any detailed comparisons. While some 
DLA galaxies have been imaged with HST (Warren et al. 2001; M0ller et al. 2002) 
and likewise GRB hosts (Jaunsen et al. 2003; Fynbo et al. 2005), and sub-mm 
galaxies (Small et al. 2004; Pope et al. 2005; Schmitt et al. 2006), only a very small 
subset of the Lya selected galaxies have been imaged with HST (Pascarelle et al. 
1996; Venemans et al. 2005; Overzier et al. 2006) and furthermore, such images are 
mostly too shallow for a detailed study. 

The GOODS-S (Giavahsco et al. 2004) provides a unique opportunity to obtain a 
deep, high resolution, multi-band data set of a complete and unbiased sample of Lya- 
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emitters (or LEGOs for Lya Emitting Galaxy-building Objects; M0ller & Fynbo 
(2001)). We have therefore started a program to collect a complete, unbiased sample 
of LEGOs in the GOODS-S. This allows a detailed study of the global properties, e.g. 
photometry and morphology of LEGOs, as well as SED fits including photometry 
from the large, available multi-wavelength data-set. 

This paper is organised as follows; in section 2 we present the imaging obser- 
vations, data reductions and candidate selection process. In section 3 we present 
spectroscopic observations of three candidates as well as the results from these ob- 
servations. Sections 4, 5 and 6 contain the discussion of various aspects of the 
LEGO candidate sample; first the basic characteristics of the LEGO sample, then 
the SED fitting and finally a comparison to Lyman-Break Galaxies. The conclusion 
is presented in section 7. 

Throughout this paper, we assume a cosmology with Hq = 72 km s~^ Mpc~^ 
(Freedman et al. 2001), f^m = 0.3 and = 0.7. Magnitudes are given in the AB 
system. We survey a co- moving volume of ~ 3300 Mpc'^. 



3.3 Imaging 

3.3.1 Narrow band observations and data reduction 

A 400x400 arcsec^ section of the GOODS-S field centred on R.A. = 03''32'^21.9" and 
Dec = -27°45'50.7" (J2000) was observed with FORSl on the VLT 8.2 m telescope 
Antu during two visitor mode nights on December 1-3, 2002. The log of observations 
is given in Table 1. A total of 16 dithered exposures were obtained over the two 
nights for a combined exposure time of 30000 seconds, all with the narrow band 
filter OIII/3000+51 and using the standard resolution collimator (0.2x0.2 arcsec^ 
pixels). For this setup the central wavelength of the filter is 505.3 nm with a FWHM 
of 5.9 nm which corresponds to the redshift range z = 3.131 — 3.180 for Lya. The 
transmission curve of the filter is shown in Fig. 13. 1[ The four spectrophotometric 
standards FeigellO, LDS749B, LTT3864, and LTT3218 were observed on the same 
nights. 

The observing conditions were unstable during the two nights with the seeing 
FWHM, as measured on the images, varying between 0.66" and 1.25" on the first 
night and between 1.4" and 3.3" on the second night. The images were reduced (de- 
biased and corrected for CCD pixel-to-pixel variations) using standard techniques. 
The individual reduced images were combined using a modified version of our code 
that optimizes the Signal-to- Noise (S/N) ratio for faint, sky-dominated sources (see 
M0ller & Warren, 1993, for details on this code). The 5 a detection limit of the 
combined narrow-band image as measured in circular apertures with radius twice the 
full width half maximum of point sources, i.e. with radius 1.6", is mag(AB)= 26.1. 
The combined narrow-band image is shown in Fig. 13. 2[ 
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Figure 3.1: Transmission of selection filters. The VLT FORSl narrow-band filter is 
drawn with a solid line. Dashed lines show the HST B and V filters. 



Table 3.1: Log of imaging observations with FORSl. 



date 


total exp. 


seeing 


range 


01-02.12.2002 


5.54 hours 


0.66"- 


-1.25" 


02-03.12.2002 


2.78 hours 


1.43"- 


■3.30" 



3.3.2 Selection of LEGOs in the fields 

For the selection of LEGO candidates we used the narrow-band image as detection 
image and the HST/ACS B- (F435W) and l^-band (F606W) images as selection 
images. The HST data is part of the public data in GOODS-S (Giavalisco et al. 
2004). This selection set-up, with a broad-band filter on either side of the narrow- 
band filter, appears to be one of the most eflicient configurations for selection of 
emission-line objects (Hayes & Osthn 2006; Hayes priv. communication). The HST 
images were re-binned to the pixel size of the narrow-band image. Due to the 
smaller field-of-view of the HST/ACS images, the narrow-band image was cut into 
six sub-images to match the size of the HST images. 

Our selection method consists of three consecutive steps. First, using the soft- 
ware package SExtractor (Bertin & Arnouts 1996), we select all objects identified in 
the narrow-band image. The narrow-band image is scanned with a detection thresh- 
old equal to the background sky-noise and requiring a minimum area of 5 connected 
pixels above this threshold. Centred on each candidate object we then extract pho- 
tometry from the narrow, V, and i?-band images using identical circular apertures 
of 2" diameter. Note that through this process we make no attempt to identify 
broad band counterparts to the narrow band objects. The broad band photometry 
is extracted in apertures defined solely from the centroid positions in the narrow 
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Figure 3.2: The VLT/narrow-band image of the 400x400 arcsec^ field with the 
positions of selected LEGO candidates (see Sect. 13.3.21) shown with circles. Spec- 
troscopically confirmed candidates are marked with crosses and their redshifts are 
indicated. North is up and East is left. The right, uppermost spectroscopically 
confirmed candidate is the Lya blob (Nilsson et al. 2006a). 



band image. There is always a small but finite possibility that an unrelated for- 
or background object could fall inside the 2" aperture which would complicate the 
search for emission line objects. This complication is minimised by not re-centering 
the aperture on objects in the broad band images. 

The second step is to accept only candidates that are detected at S/N> 5 in 
the narrow band circular aperture and are found at least 20 pixels from the edge of 
each image. This leaves us with a catalogue of 2616 narrow band objects within the 
resulting 385 x 400 arcsec^ field. The term "narrow band objects" is used here only 
to underline that while they are 5a detections in the narrow band image, they may 
or may not have been detected in the broad band images. 

In the third step we select the subset of our catalogue comprising of potential 
emission line objects. Via interpolation of the fiux levels in the B and V bands 
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it is easy to calculate the continuum flux at the central wavelength of the narrow 
band filter, and from there to obtain the equivalent width (EW) of a potential line 
and its associated propagated statistical error. Note that we here calculate the 
"EW of the aperture" , which means that if there is only one object in the aperture, 
then we find the EW of that object. If there are additional unrelated neighbours 
inside the aperture, then the calculated EW will be smaller than the actual EW. 
Our listed EWs (see Table 13. 2p are in that sense conservative lower limits. For the 
present work we are interested only in those with positive values of EW (emission 
line objects), and only the subset of those where the significance of the line is 
high enough to provide a high probability that it is reliable, i.e. providing a high 
efficiency of spectroscopic follow-up work. For the current field we already have a 
few spectroscopic confirmations (see Sec. 13. 4p and we conservatively chose to cut 
at the EW significance level {2.9a) where all confirmed LEGOs are included. This 
corresponds to a formal probability of 99.6% for confirmation of each object and 
we find 106 such objects. Of these, three are associated with a Lya-blob that we 
found in this field (Nilsson et al., 2006a). Lya-blobs (e.g. Fynbo et al. 1999; Steidel 
et al. 2000; Matsuda et al. 2004) are large luminous nebulae, with sizes up to 
150 kpc, emitting solely Lya emission. The Lya luminosity can reach 10^^ erg s~^. 
79 of the 106 objects with EW excess are ruled out after visual inspection due to 
stellar artifacts, saturated sources or for lying on parts of the image where the HST 
sky background is poorly constrained. Thus, we are finally left with a sample of 
24 LEGO candidates. These candidates are presented in Table 13.21 The formal 
probability of 99.6% for confirmation of the emission line results in 0.10 spurious 
detections in our sample of 24 objects. A colour-colour plot of the narrow-band 
sources can be found in Fig. 13.121 

3.3.3 Continuum counterparts and final photometry 

Following the selection process outlined above, we examined the narrow-band and 
broad-band images to find continuum counterparts to the narrow-band objects. The 
process of finding counterparts to the narrow-band objects is complicated by the 
very different PSF characteristics of the narrow- and broad-band images. This is 
illustrated in Fig. 13.31 For some of the candidates, no obvious counterpart was 
detected, but rather several counterpart candidates were found with small offsets. 
To determine what continuum objects were associated with the narrow-band source, 
two co-authors separately inspected the images visually. If only one counterpart was 
in the vicinity of the centroid of the narrow-band source, this object was identified 
as sole counterpart. If several sources were detected, their magnitudes were mea- 
sured and the statistical probability that they would appear in the area surrounding 
the narrow-band centroid was evaluated from number counts of galaxy searches. We 
then separately determined which counterparts we considered credible counterparts, 
and the lists were compared. Only counterparts assigned by both authors were ac- 
cepted as counterparts. This yielded 2 LEGO candidates without counterparts, 14 
candidates with single counterparts and 8 LEGOs with two or more counterparts. 
Aperture photometry was then performed on all candidates in the narrow-band 
and their selected counterparts in the HST broad band images. We used aperture 
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Table 3.2: Data on first selection. Coordinates are in J2000. Equivalent widths 
are calculated from the "first selection", i.e. with 2" radius apertures, centred on 
the narrow-band source centroid, see section [3.3.21 S/N is the signal-to- noise of the 
equivalent width. 



LEGO 


EWobs 


0"EW 


O /AT 

S/N 


K.A. 


Dec 


GOODS-S 


(A) 


(A) 








1 


896 


57 


16 


03:32:14.83 


-27:44:17.5 


2 


59 


15 


3.8 


03:32:17.62 


-27:43:42.3 


3 


172 


24 


7.1 


03:32:18.56 


A ' A A A 

-27:42:48.4 


4 


901 


48 


19 


03:32:31.46 


-27:43:37.2 


5 


184 


30 


6.0 


03:32:30.02 


-27:48:37.1 


6 


517 


23 


22 


03:32:30.82 


-27:47:52.8 


7 


153 


32 


4.7 


03:32:13.40 


-27:47:43.9 


8 


85 


20 


4.3 


03:32:30.79 


-27:47:19.2 


9 


100 


24 


4.2 


03:32:12.49 


-27:42:45.8 


10 


52 


16 


3.3 


03:32:13.30 


-27:43:29.9 


11 


154 


33 


4.6 


03:32:27.03 


-27:47:28.3 




ioi 


66 


A P. 
4.0 


Uo.oz. io.yy 


-Z / .44.4 / .U 


13 


53 


18 


2.9 


03:32:14.58 


-27:45:52.5 


14 


37 


4 


8.3 


03:32:18.82 


-27:42:48.3 


15 


88 


11 


7.9 


03:32:31.56 


-27:46:26.9 


16 


34 


7.5 


4.4 


03:32:20.72 


-27:42:33.8 


17 


202 


46 


4.4 


03:32:28.93 


-27:45:31.5 


18 


269 


55 


4.9 


03:32:17.26 


-27:45:21.0 


19 


47 


12 


3.9 


03:32:15.80 


-27:44:10.3 


20 


58 


15 


3.9 


03:32:30.96 


-27:43:14.2 


21 


53 


16 


3.4 


03:32:28.73 


-27:47:54.1 


22 


58 


15 


4.0 


03:32:17.77 


-27:42:52.1 


23 


362 


57 


6.3 


03:32:15.37 


-27:48:18.5 


24 


136 


24 


5.7 


03:32:18.89 


-27:47:03.6 



radii of two times the FWHM of each image. Aperture corrections were calculated 
for each image for point sources (see Table 13.51) . For candidates where multiple 
components were assigned, small apertures were placed either i) centred on each 
counterpart separately if the counterparts are further apart than two times the ra- 
dius of the aperture, ii) centred on a coordinate half way between the counterparts 
if the distance between counterparts is less than one times the radius of the aperture 
or Hi) apertures slightly shifted from the central coordinate of the counterparts to 
ensure that different apertures do not overlap if the distance between counterparts 
is less than two, but more than one, times the aperture radius. The magnitudes, 
EWs and star formation rates (SFRs) for all candidates can be found in Table 13. 3[ 
see also Sec. 13.5.11 Multiple candidates are marked with a star. To investigate 
how correct our method is for measuring fluxes of multiple objects, we also mea- 
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Figure 3.3: In this figure, we illustrate the difficulty in determining which broad- 
band counterpart was associated with the narrow-band source. Images are 12" across 
and are centred on the narrow-band source. The circles mark the apertures used 
in the initial photometry, see section I3.3.2[ The two left panels show a simple case 
with only one possible counterpart. The two right panels show a more complex case. 
In this case, all objects within the circle were assumed to be counterparts. In both 
cases a presumably unrelated object is seen at the edge of the aperture. 



sured the photometry using larger apertures; for the six multiple counterparts with 
distances between the counterparts less than two times the radius of the aperture 
(LEGO GOODS-S_# 9, 10, 12, 14, 19 and 22), we applied new apertures with 
radii half of the distance between counterparts plus the original aperture radius. 
For five candidate counterparts, the difference in flux measured was within 1.5a 
of the previously measured value. Hence we conclude that our original measure- 
ments are correct. The remaining candidate, LEGO GOODS-S_# 14, has a very 
complicated morphology, see also Fig. 13.31 For this candidate, the flux increased to 
Mb = 25.49 ± 0.11, My = 24.47 ± 0.04, Mj = 24.36 ± 0.13 and My = 24.49 ± 0.13. 
This reduces the observed EW to 91 ± 8 A. 

In summary we first performed simple circular aperture photometry in order 
to select candidates based on our conservative emission-line definition (Table [372]) . 
For the selected candidates, we then searched for continuum counterparts in the 
high resolution HST images and carried out detailed final photometry where such 
counterparts were found. This final photometry provides the relevant magnitudes 
that describe the objects and is reproduced in Table 13. 3[ 



Table 3.3: Final photometry of candidates in narrow- and broad-bands (Sec. I3.3.3P and observed Lya EWs. Magnitudes are calculated 
using two times full width half maximum apertures, including aperture corrections, centred on each identified counterpart. Errors are la, 
upper limits are 3a. Equivalent widths are calculated again from the magnitudes printed in this table. Star formation rates are from the Lya 
fluxes. Emitters marked in bold are spectroscopically confirmed, see Sec. 13.41 Candidates marked with a star have multiple counterparts. 
For these candidates, the total magnitude is given here. 



ID 


^^narrow 


Mb 


My 


Mi 




EWobs.Lya (A) 


SFRLya 

(Me/yr) 


1 


24.35 


± 





.07 


27.87 ±0.17 


27.08 ±0.07 


27.33 ±0.28 


27.07 ±0.22 


1006 ± 71 


3.94 ±0.24 


2 


25.36 


± 





17 


28.02 ±0.35 


27.29 ±0.10 


27.39 ±0.26 


> 27.52 


434 ± 82 


1.55 ±0.22 


3 


24.97 


± 





06 


> 28.08 


28.05 ±0.38 


> 27.67 


> 27.24 


912 ± 52 


2.22 ±0.11 


4 


24.19 


± 





.02 


27.63 ±0.02 


26.91 ±0.11 


27.25 ±0.30 


27.07 ±0.32 


956 ± 15 


4.55 ±0.07 


5 


25.29 


± 





03 


> 28.32 


28.45 ± 0.40 


> 28.68 


> 26.89 


881 ± 27 


1.65 ±0.04 


6 


24.03 


± 





04 


> 28.23 


27.12 ±0.15 


> 27.52 


> 27.31 


1639 ± 64 


5.28 ±0.19 


7 


25.45 


± 





21 


> 28.05 


27.47 ± 0.09 


27.34 ±0.28 


> 27.26 


436 ± 103 


1.42 ±0.24 


8 


25.23 


± 





14 


27.89 ±0.13 


27.62 ±0.21 


27.39 ± 0.37 


> 27.37 


532 ± 81 


1.74 ±0.21 


9* 


24.75 


± 





08 


26.46 ± 0.09 


25.98 ±0.06 


25.64 ±0.13 


26.03 ±0.15 


169 ± 16 


2.73 ±0.18 


10* 


25.02 


± 





12 


26.84 ± 0.09 


26.21 ±0.09 


26.24 ±0.08 


26.50 ±0.20 


179 ± 27 


2.13 ±0.22 


11 


25.28 


± 





16 


> 28.27 


> 28.36 


> 27.86 


> 27.60 


> 834 


1.66 ±0.22 


12* 


25.97 


± 





24 


> 27.76 


26.11 ±0.09 


27.43 ± 0.29 


> 27.46 


76 ±34 


0.89 ±0.18 


13 


25.63 


± 





.23 


27.24 ±0.15 


27.13 ±0.12 


27.24 ±0.28 


> 27.08 


178 ±57 


1.21 ±0.23 


14* 


23.97 


± 





05 


25.64 ± 0.05 


24.58 ± 0.03 


24.39 ± 0.04 


24.64 ± 0.06 


110 ±8 


5.58 ±0.26 


15* 


24.22 


± 





05 


26.59 ±0.09 


25.73 ±0.04 


25.25 ±0.04 


25.22 ±0.06 


296 ± 18 


4.42 ±0.21 


16 


24.92 


± 





11 


> 27.72 


26.87 ±0.12 


26.12 ±0.08 


25.46 ± 0.08 


473 ± 59 


2.33 ±0.23 


17* 


25.09 


± 





08 


26.94 ±0.16 


25.86 ±0.05 


25.62 ±0.08 


25.33 ±0.08 


138 ± 15 


1.98 ±0.14 


18 


25.37 


± 





16 


27.39 ±0.12 


26.11 ±0.04 


25.73 ± 0.04 


25.84 ±0.10 


149 ± 32 


1.54 ±0.21 


19* 


25.43 


± 





07 


27.02 ± 0.06 


26.01 ±0.03 


25.80 ± 0.08 


26.16 ±0.07 


102 ± 11 


1.45 ±0.10 


20 


25.99 


± 





15 


27.08 ±0.14 


27.04 ± 0.09 


26.23 ±0.11 


26.02 ±0.13 


90 ±23 


0.87 ±0.11 


21 


25.48 


± 





06 


> 27.88 


> 27.86 


> 27.44 


> 27.35 


> 443 


1.38 ± 0.07 


22* 


25.05 


± 





15 


26.80 ±0.12 


25.65 ±0.04 


25.65 ±0.08 


25.60 ±0.08 


115 ±26 


2.07 ±0.27 


23 


25.30 


± 





15 


> 27.78 


25.99 ±0.05 


25.47 ±0.04 


25.43 ±0.07 


174 ± 35 


1.64 ±0.21 


24 


25.46 


± 





12 


27.12 ±0.10 


26.51 ±0.06 


26.02 ±0.07 


26.04 ±0.12 


146 ± 25 


1.41 ±0.15 



42 Chapter 3 

3.4 Spectroscopy 



3.4.1 Observations and reductions 

Follow-up Multi-Object Spectroscopy (MOS) was obtained in service mode with 
FORSl/VLT UT2 over the time period December 2004 - February 2005. The total 
observing time of 5.5 hours was granted to confirm the redshift of the Lya-blob 
found in this field, as published in Nilsson et al. (2006a). In addition, we had 
the opportunity to add three of our compact emitters on the mask. The mask 
preparation was done using the FORS Instrumental Mask Simulator. Stars were 
placed on the remaining slits for calibration purposes. The MOS slitlets had a 
width of 1.4" and the combination of grism 600V and order sorting filter GG435 
was used. The grism covers the wavelength range 4650 A to 7100 A with a resolving 
power of approximately 700. The seeing varied between 0.77" - 1.2". 

The bias subtraction, fiat-fielding and wavelength calibration was performed 
using the FORSl pipeline. The individual, 2-dimensional, reduced science spectra 
were combined using a cx-clipping for rejection of cosmic ray hits. The sky was 
subsequently subtracted with MIDAS by averaging the values of all pixels on either 
side of the spectrum and expanding this value to the size of the frame. One- 
dimensional spectra were extracted by summing the column values over the spectra. 
The spectra were then flux calibrated with three stars that were on our slits, by 
measuring the stellar fluxes in the narrow-band image in 2" diameter apertures and 
comparing to the integrated flux in the narrow-band from the 1-d spectra. Because 
LEGOs are often extended in Lya (M0ller & Warren 1998; Fynbo et al. 2001) they 
are likely to have higher Lja fluxes than those listed in Table 13. 4[ 



3.4.2 Results of first spectroscopic follow-up 

The spectra of our three confirmed high-redshift Lya-emitters can be found in 
Fig. 13. 4[ Neither candidate show any other emission lines in their spectra. To 
consider if the emission line could be [Oil], we study line ratios compared to other 
lines that should be observed in such a case. In Fynbo et al. (2001), the expected 
line ratios of H/3, [Oil], [OIII] and Nell for [OII]-emitters were presented. For the 
spectroscopic sample, the 3a upper limits of the ratio log(F[o///]/F[o//]) are given in 
Table [33] if the emission lines is [Oil]. The limit will be the same for all other lines, 
as none are detected. As in Fynbo et al. (2001), we conclude that these are highly 
unlikely values for [OII]-emitters, and that the detected emission line is redshifted 
Lya. Details from the spectroscopic follow-up is presented in Table [131 

In the following sections we analyse the entire sample of confirmed LEGOs and 
LEGO candidates together. We expect the contamination of low redshift emitters 
to be small. Our previous surveys have had spectroscopic success rates of 75 - 90 % 
(Fynbo et al. 2001; Fynbo et al. 2003). Hence, we expect that more than 18 of our 
24 candidates are true Lya-emitters. 
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Figure 3.4: Spectra of confirmed LEGOs in GOODS-S. a) LEG0_G00DS-S#1 
emission line, b) LEGO-GOODS-S^^l expected position of [OIII] lines if the detected 
emission line is [Oil], c) LEG0_G00DS-S#4 emission line, d) LEGO_GOODS- 
S#4 expected position of [OIII] lines if the detected emission line is [OII],ej 
LEG0_G00DS-S#13 emission line. Our spectrum did not cover the position of 
[OIII] if the emission line is [Oil] for LEG0_G00DS-S#13. Hatched areas mark 
the positions of bright sky lines. 



Table 3.4: Names, Lya-fluxes, redshifts and line ratios for spectroscopically con- 
firmed LEGOs. The line ratios refer to the upper limit to the [OIII] line, if the 
emission line is [Oil]. 



LEGO_GOODS-S# 


FLya 




Z 


log(F[o7//]/F[o//]) 




(erg cm 










1 


3.31 X 10" 


17 


3.151 


< 


-0.62 


4 


2.79 X 10- 


17 


3.132 


< 


-0.54 


13 


2.84 X 10" 


17 


3.118 


< 


-0.55 



3.5 Basic characteristics of LEGOs 
3.5.1 SFR, surface density and sizes 

For our final sample of LEGO candidates, we calculate the star formation rate (SFR) 
as derived from Kennicutt (1983) by: 

SFR = .. Mgyr-^ (3.1) 

1.12 X 10^2 ^ ^ 

where the Ha luminosity, I/Ha, is obtained with the conversion between Lya and Ha 
luminosities of Brocklehurst (1971) of L(Lya) = 8.7x L(Ha). The SFR values of the 
LEGO candidates can be found in Table 13. 3[ The mean SFR, as derived from the 
Lya-emission for all candidates is 1.8 M0/yr. The total SFR is 43 Mq/jt, yielding 
a star formation rate density Psfr of 0.013 Mo/yr/Mpc^. This value is in very 
good agreement with other results for high redshift galaxies at this redshift of e.g. 
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Figure 3.5: Histogram plot over the size of our candidates in the narrow-band image. 
The bin size is 0.1 arcsec, the sohd hne represents the PSF of the image and the 
dotted hues the la error on the PSF. The PSF was determined from 28 objects 
with SExtractor keyword CLASS_STAR greater than 0.9 and fluxes in the range of 
our LEGO candidates. The object in the highest FWHM bin is the GOODS-S blob 
(Nilsson et al., 2006a). 



Madau et al. (1996; 0.016 Mo/yr/Mpc^), Steidel et al. (1999; 0.05 Mo/yr/Mpc^) 
and Cowie & Hu (1998; 0.01 Mg/yr/Mpc^). The results from Steidel et al. (1999) 
has been obtained from integrating the extrapolated luminosity function down to 
a luminosity of 0.1 Lj,, and the results are also corrected for dust by multiplying 
by a factor of 4.7. Hence, their results uncorrected for dust is 0.011 M0/yr/Mpc^. 
There is very good agreement between the dust uncorrected measurements of the 
SFR density at 2 ~ 3. 

We find a surface density of LEGOs at redshift z = 3.15, = 0.05 in the 
GOODS-S field of 0.53 objects per arcmin^. We can compare this with all V band 
sources in the GOODS-S field by extracting all sources with V band magnitudes 
between 25 - 28 in the available online catalog. We find 23885 such sources in 
the entire GOODS-S area, covering approximately 160 arcmin^, corresponding to a 
surface density of 149 arcmin"^. If we assume a homogeneous density of LEGOs 
between redshift 3.0 — 3.5, then the surface density of LEGOs, scaled with our 
candidate sample, will be ~ 5.3 arcmin"^. Thus, approximately 4 % of all V- 
detected sources with a magnitude of = 25 — 28 were selected as Lya-emitters in 
the redshift range 2; = 3.0 — 3.5. 

We measured the sizes of our candidates in the narrow-band images using the 
FLUX_RADIUS option in SExtractor. This gives the half width half maximum of 
each source. Object LEGO_GOODS-S# 11 was excluded because it was blended 
with another, unrelated object. The histogram of the sizes, compared to a point 
source, is presented in Fig. 13.51 Most objects appear to be barely resolved. However, 
there is a tail of larger objects extending towards the GOODS-S blob. In a future 
paper, we will present a complete morphological study of the candidate sample. 
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Figure 3.6: Left Histogram distribution of distances between candidate objects in 
the x-direction, after re-alignment by 2.1°. Histogram is binned over 50 pixels. Solid 
curve is a double gaussian fit to the data. Right Kolmogorov-Smirnov test of the 
distribution. Lines mark the simulated mean, and the 1 - 5a contours. Thick line 
represents our data. 



3.5.2 Filamentary structure 

Previous studies of Lya emitters have reported the identification of filamentary 
structures (e.g. M0ller & Fynbo 2001; Hayashino et al. 2004; Matsuda et al. 2005). 
The volume we survey here is approximately 3 times as long along the line-of-sight 
as it is wide which means that several filaments could be crossing the volume at 
different redshifts and with different position angles. If that was the case, they would 
likely blend together to wash out individual structures. Bearing this in mind we can 
still ask the question of whether the objects on our candidate list are randomly 
distributed across the field, or if they appear to be systematically aligned. 

Inspecting Fig. 13.21 we note that the candidates do in fact appear to be aligned 
in two filamentary structures approximately along the y-axis of the CCD (oriented 
N-S). To investigate whether this is a significant effect, we have calculated the 
distances between all objects and plotted a histogram of the projected distances 
between all candidates in Fig. 13. 6[ The histogram shows a definite division in 
two components, one describing the typical width of a filament, and the other the 
projected distance between the filaments. To optimise the search, we calculated 
the angle of rotation that minimised the variance of the x-coordinates around a 
mean. For the two filaments, these angles are 1.75° and 2.52° respectively. For 
the analysis, we rotate the filament by the average, 2.13°. We then generated 10'^ 
uniform random sets of coordinates, with the same number of objects, and repeated 
the same analysis of calculating distances between pairs. To establish how reliable 
our observed distribution is, we perform a Kolmogorov-Smirnov test (e.g. Peacock 
1983) on the simulated distributions. The test showed that the likelihood of the 
alignment being random is less than 2.5 x 10~^, hence a near 4 sigma detection, 
see Fig. 13.61 We fitted a double gaussian function to the histogram plot. This 
gives the typical width of the filaments as the FWHM of the first peak, and the 
distance between them as the mean distance to the second peak. We find that the 
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typical width is ~ 250 pixels, corresponding to 370 kpc at redshift 3.15, in good 
agreement with the findings of M0ller & Fynbo (2001) who find a spectroscopically 
confirmed filament with 400 kpc radius from a set of Lya-emitters at z = 3.04. 
The distance between the first and second peak is ~ 950 pixels, corresponding to 
1.4 Mpc at this redshift. To further verify this filamentary structure would require 
more spectroscopic data, that would also enable a 3D-plot of the filaments in space. 

3.6 SED fitting 

The imaging available in the GOODS-S field is extensive. The data, in 14 publicly 
available broad-bands, used here is presented in Table 13.51 



Table 3.5: Deep, multi-wavelength data available in the GOODS-S field. The fifth column refers to the 3a detection limit in the sky in a 
2x FWHM diameter aperture. The last column gives the 3a detection limit as measured in a 2" radius aperture. 



Filter/Channel 


Ac 


FWHM 


Aperture 


Aperture 


3(7 limit 


3a limit (2 


" aperture) 








radius 


correction 


(2x FWHM aperture) 
(erg • cm"^ • s"^ • Hz"^) 












(arcsec) 




(erg • cm-2 


• s-i • Hz-i) 


X-rays (Chandra) 


4.15 keV 


3.85 keV 


2.25 


— 


9.90 • 10-34 


9.90- 


10-34 


U (ESQ 2.2-m) 


3630 A 


760 A 


3.00 




1.10 • 10-30 


8.62 • 


10-31 


B (F435W, HST) 


4297 A 


1038 A 


0.12 


1.20 


7.09 • 10-32 


9.25 • 


10-30 


V (F606W, HST) 


5907 A 


2342 A 


0.12 


1.18 


4.02 • 10-32 


4.66- 


10-30 


i (F814W, HST) 


7764 A 


1528 A 


0.12 


1.25 


1.25 ■ 10-31 


1.50- 


10-29 


z' (F850LP, HST) 


9445 A 


1230 A 


0.12 


1.34 


1.88 • 10-31 


3.00 • 


10-29 


J ( VLT) 


1.25 /im 


0.6 /im 


0.60 


1.22 


1.78 • 10-30 


5.31 ■ 


10-30 


H ( VLT) 


1.65 /xm 


0.6 /xm 


0.60 


1.21 


4.11 • 10-30 


1.86 • 


10-29 


Ks {VLT) 


2.16 jjLin 


0.6 /tm 


0.60 


1.37 


4.06 • 10-30 


1.56- 


10-29 


Chi (Spitzer) 


3.58 fim. 


0.75 /tm 


1.30 


2.41 


1.14-10-31 


2.36- 


10-30 


Ch2 (Spitzer) 


4.50 ^m 


1.02 /xm 


1.80 


1.98 


2.07 • 10-32 


2.07- 


10-30 


Ch3 (Spitzer) 


5.80 /im 


1.43 /xm 


1.80 


1.60 


7.50 • 10-29 


7.50 • 


10-30 


ChA (Spitzer) 


8.00 iim 


2.91 /tm 


2.10 


1.85 


6.87 • 10-30 


6.87- 


10-30 


MLPS (Spitzer) 


24.0 ixm. 


4.70 /tm 


6.00 




1.23 • 10-28 


2.12- 


10-29 



i. 

CO 



o 
o 

Co 
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Table 3.6: Stacked magnitudes for the GOODS-S Lya-emitters. Errors in the HST 
magnitudes were set to a conservative value of 0.08. Upper limits are 3a. 



Band 


Uentr. Waveiength (Aj 


(Jbs. Magnitude 


Mag. brror 


U 


3710 


> 25.95 




B 


4297 


27.57 


0.08 


IN 


DVOO 


OA Cifi 

z4.yo 


U.Uo 


V 


5907 


26.74 


0.08 


i 


7764 


26.52 


0.08 


z 


9445 


26.56 


0.08 


J 


12500 


> 26.28 




H 


16500 


> 25.55 




Ks 


21500 


25.26 


0.29 


Chi 


35800 


> 23.06 




Ch2 


45200 


> 23.62 




Ch3 


57200 


> 24.37 




Ch4 


79000 


> 23.84 





With this data-set, we wish to perform an SED fitting, in order to constrain 
properties such as stellar mass M*, dust content Ay, metallicity and age of the LE- 
GOs. Only one of our candidates (LEG0_G00DS-S#16) is detected in bands other 
than the HST bands. This object is especially interesting as its SED is extremely 
red. It is excluded from the SED fitting, and is discussed in Section 13.6.31 For the 
rest of the sample, the LEGOs are only detected in the HST bands and hence we 
choose to stack the entire sample of 23 candidates. We can then draw conclusions on 
the general properties of this type of object. After stacking, we get a faint detection 
in the Kg band. The stacked magnitudes are given in Table 13.61 The lack of X-ray, 
MIPS 24/im and radio detections (no counterparts to any of our candidates to a 3a 
limit of 24 /iJy, Kellermann et al. in preparation) indicates that the AGN fraction 
among these objects is low. 

3.6.1 Fitting method 

We used the GALAXEV code (Bruzual & Chariot, 2003) to simulate composite 
stellar populations, in order to fit the stacked SED of the LEGOs. The fitting was 
performed according to a Monte Carlo Markov Chain method (see e.g. Gilks et 
al. 1995 for an introduction). In outline, the method works as follows; an initial 
set of parameter values is chosen according to a uniform, random and logarithmic 
distribution within the allowed parameter space. A summary of the parameter 
space is given in Table [37n Given the set of parameters, a corresponding value 
is calculated by running the GALAXEV code, creating a high-resolution spectrum 
with 6900 wavelength points from 91 A to 160 fim. To obtain the magnitudes in 
each band, we apply the transmission curves for the filters of the various observed 
wavebands; U, B, V, i, z', J, H, Kg and the four Spitzer bands, Chi — Chi, In this 
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Table 3.7: Parameter space sampled during the SED fitting. Metallicity is allowed 
to have three different values {Z/Zq = 0.005, 0.2 or 1.0). The dust components are 
the two components of the Chariot & Fall (2000) dust model used by GALAXEV. 
We fit the SED with a constant star forming model, where the star formation rate 
in solar masses per year is given by "SF-rate" . 



Parameter 


Min. value 


Max. value 


MetaUicity (Z/Zq) 


0.005 


1.0 


Dust-r 





4 


Dust-/^ 





1 


SF-rate 


0.01 


100 


Age (Gyrs) 


0.001 


1.5 



analysis, we exclude the Spitzer MIPS band as it is very difficult to stack images in 
this band because of the source confusion. At this redshift, the MIPS band is also 
contaminated by a PAH emission feature. We then normalise the output spectrum 
so that the magnitude in the model 2;'-band equals that of the observed z'-band. 
The is then calculated by comparing the magnitudes in all the other bands. We 
incorporate points with upper limits in the following way; if the predicted flux lies 
below the upper limit, no value is added to the total x^y if the predicted flux lies 
above the observed one, a is added, assuming the error on the upper limit is 0.1 
in magnitude. Hence, models with flux above the upper limit may be acceptable 
if the flux in this band is only slightly above the limit. In most cases though, the 
model will be rejected due to high x^. 

Once the has been calculated for a particular model, a new random set of 
parameters is chosen, by adding a "step vector" . The step size in each parameter is 
chosen randomly in a logarithmic interval between 1% and 100% of the total size of 
the parameter space. This step vector is equally likely to be positive or negative. The 
choice of logarithmic step sizes is a natural choice if no assumptions about the scale 
of change that affect the solution are to be made and ensures a fast convergence. 
When a step has been calculated and the new parameters have been calculated, 
the new model is accepted into the chain with a probability proportional to the 
exponential difference between the old and new x^. If the new model is accepted, 
it is added as the next step in the chain, i.e. the parameters of that particular 
model are printed in the output flle. If the model is not accepted, a copy of the old 
parameters is added as the next step. This procedure is repeated until a chain with 
30000 elements has been created. The independence of individual steps in the walk 
ensures that the resulting chain is Markovian in character, i.e. that the resulting 
chain after many iterations is a representation of the full probability distribution 
function in the chosen parameter space. The output flle can then be used to study 
the distribution in each parameter, and to determine the mean and the confldence 
levels within each parameter. It can also be used to study dependencies between 
parameters, such as e.g. in Fig. 13.71 
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3.6.2 Results from SED fitting 

The parameters we wish to fit are stellar mass, dust content, star formation rate, 
metallicity and age. Redshift is set to be the central Lja redshift of the narrow- 
band filter, i.e. z = 3.15. We use the Salpeter initial mass function (IMF) from 
0.1 to 100 M0 and the extinction law of Chariot & Fall (2000). We also in- 
corporate the effect of the Lja forest according to the model of Madau (1995). 
We use constant star formation histories. The metallicity was allowed to be ei- 
ther Z/Zq = 0.005, 0.2 or 1.0. Of the 30000 runs, 73% were with the low- 
est metallicity, 21% with the medium metallicity and 6% had solar metallicity. 
Hence the best fit models have a very low metallicity. For the rest of the anal- 
ysis, we choose to only look at the models which have the lowest metallicity, as 
these models seem to be preferred. For these models, the best fit parameters are 
M, = A.rtti X 10^ Mq, age = O-Sbtoll Gyrs, Ay = 0.26^^, and star formation 
rate SFR = O.GGlgg} Mq yr~^, where the errors are la. Contour-plots of the three 
parameters mass, age and dust are shown in Fig. 13. 7[ The degeneracies between the 
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Figure 3.7: Contour-plots of the mass, age and dust parameters in our SED fitting. 
Contours indicate 1, 2 and 3a levels. 

different parameters can be seen. In Fig. 13.81 the weighted GALAXEV spectrum of 
a subset of 10 models with ~ 1 is shown with the measured SED overplotted. 
The ?7-band and Spitzer/IRAC data are too shallow to be useful in constraining the 
data, and the narrow-band data-point is plotted only for reference. 

3.6.3 Object LEGO_GOODS-S#16 

One candidate emission hue object, LEG0_G00DS-S#16, was detected in all avail- 
able GOODS-S bands, except the U and B bands, and in X-rays. The fluxes of the 
object can be found in Table 13.81 and the thumb-nail images seen in Fig. 13. 9[ As 
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Figure 3.8: The average spectrum of 10 models with ~ 1- The spectrum is 
calculated by 1) creating the 10 spectra with the particular parameters of the models 
with GALAXEV, 2) averaging the 10 spectra together, weighted with the x^, so 
that the total flux density is F^, = [J^i^u ^ / / xl) ■ The parameters of 

the model spectrum are as indicated on the plot. They are the weighted average 
parameters of the 10 models used, weighted in the same way as the spectra. Data 
points from stacked SED. Upper limits are represented with arrows. The point well 
off the SED is the narrow-band magnitude. 
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Figure 3.9: Thumb-nail images, 18" across, of LEG0_G00DS-S#16. 



there is significant excess emission in our narrow-band filter, we expect the source 
to be either an [OII]-emitter aX z = 0.36 or a Lya-emitter aX z = 3.15. This object 
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Table 3.8: SED of object LEG0_G00DS-S#16. Upper limit is 3a. Continuum 
sources are offset from the narrow-band source by approximately 0.8 arcseconds, 
corresponding to 5.9 kpc at z = 3.15. 



Band 


Centr. Wavelength (yumj 


(Jbs. nux (/Ujy) 


Til / T \ 

r lux Error [jJoy) 


B 


0.430 


> 0.03 






U.oUd 


U.oy 


U.U4Z 


V 


0.591 


0.06 


0.008 


i 


0.776 


0.13 


0.010 


z 


0.945 


0.24 


0.018 


J 


1.25 


3.47 


0.157 


H 


1.65 


7.40 


0.420 


Ks 


2.15 


12.51 


0.565 


Chi 


3.58 


24.87 


0.188 


Ch2 


4.50 


30.86 


0.200 


Ch3 


5.80 


31.85 


0.710 


Ch4 


8.00 


21.31 


0.990 


MIPS 


24.00 


271.11 


7.001 



was first fit with the same type of stellar SED as the other sample, but with the 
redshift set to be either z = 0.36 if the narrow-band emission is [Oil] ot z = 3.15 if 
the emission is Lya. This fitting yielded no good fit, with ^ 500. 

Two types of objects could show MIR colours similar to those of this galaxy; i) 
obscured AGNs (e.g. Lacy et al. 2004; Stern et al. 2005b) and li) ULIRG/dusty 
starburst galaxies (e.g. Ivison et al. 2000; Klaas et al. 2001). The infrared colours 
can be plotted in the diagnostic colour-colour diagram of Ivison et al. (2004), see 
Fig. 13.101 In this diagram, Ivison et al. (2004) plot the colours of Arp 220, Mrk 
231 and a theoretical starburst spectrum as observed at different redshifts. The 
comparison with the colours of LEGO_GOODS-S7^16 shows the object to be more 
likely a low-redshift starburst galaxy. However, Ivison et al. (2004) do not take PAH 
emission into account. The most important PAH lines are at 3.3, 6.2, 7.7, 8.7, and 
11.2 /im (corresponding to 13.7, 25.7, 32.0, 36.1, 46.5 at z = 3.15). Especially the 
second line falls on top of the MIPS 24 /xm band. This would explain the extreme rise 
in flux in this band. In Fig. I3.14[ we see that this emission could easily add a factor of 
ten to the flux in this band, hence the 24 /im/8 fim colour of LEG0_G00DS-S#16 
can be adjusted downwards by a factor of ten. The extra line emission explains why 
a 2; ~ 3 star-burst galaxy would appear to be at lower redshifts. To understand if 
this is a low or high redshift object, we have also attempted to get a photometric 
redshift estimate for this galaxy using the HyperZ code (Bolzonella et al. 2000), 
both with and without including the Spitzer data points. The results were similar 
for both cases. When we do not include the Spitzer data, the best flt is redshift 
z = 1.7 with a of approximately 8. For the redshift z = 0.4, the flt then has a 

^ 126 and for the Lya redshift of z = 3.15, the ~ 35. When the Spitzer points 
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Figure 3.10: Diagram of Ivison et al. (2004). The x- and y-axes show the colours 
in the Spitzer bands and the sohd hues mark the locations of AGN and starburst 
galaxies (SB) depending on redshift. The redshifts are marked along the lines. The 
orange dots mark the location of a set of sub-mm galaxies presented in Ivison et al 
(2004). The sohd square marks the location of the colours of LEG0_G00DS-S#16, 
indicating a lower redshift starburst galaxy. The open square marks the colours of 
this galaxy if the MIPS 24 /im flux is decreased by a factor of 10 (see text). This 
point is indicative of a redshift 2; ~ 3 starburst galaxy. 



are included, the higher redshift is even more favoured. Hence, it seems unlikely 
that this is a lower redshift source. 

We wish to distinguish whether LEG0_G00DS-S#16 is an obscured AGN or 
a starburst galaxy. Several papers have presented infrared colours for obscured 
(and unobscured) AGN (Johansson et al. 2004; Lacy et al. 2004; Stern et al. 
2005b; Alonso-Herrero et al. 2006) and especially two papers publish selection 
criteria for obscured AGNs (Lacy et al. 2004; Stern et al. 2005b). The colours of 
LEG0_G00DS-S#16 are inconsistent with those selection criteria and we therefore 
rule out an AGN nature of this galaxy. This conclusion is further supported by the 
non-detection in X-rays. In order to study if the SED of the galaxy could be fitted 
by a starburst spectrum, we tried to fit the SED with a GRASIL (Silva et al. 1998) 
model of a starburst galaxy. GRASIL is a spectral stellar synthesis code, which 
takes into account the dust obscuration of starlight in both molecular clouds and 
the diffuse medium. Hence it is perfectly suited for the investigation of starburst 
galaxies. We could fit the photometric data of LEG0_G00DS-S#16 by a relatively 
old burst (~ 1 Gyr) at z = 3.15 with a significant amount of dust. The results are 
shown in Fig. 13.141 As can be seen in the Figure, these models reproduce the trends 
in the observed SED relatively well. This appears to be a redshift z = 3.15 dusty 
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starburst galaxy, with a region where the dust amount is smaller and Lya emission 
can escape, offset from the central parts of the galaxy. It would be of great interest 
to get sub-mm imaging of this object in order to constrain the SED better. 

3.7 Comparison to Lyman-Break Galaxies 

We wish to compare our sample of LEGOs to a sample of faint Lyman Break Galaxies 
(LBGs) in order to determine the similarities and differences of the two populations 
of high-redshift galaxies. First, we want to know if our LEGOs would be detected as 
LBGs and so we apply the LBG selection criteria for [/-band drop-outs of Wadadekar 
et al. (2006; U - B > 1.0, U - B > B - V + 1.3 and B - V < 1.2) as well as 
the criteria of Madau et al. (1996; U - B > 1.3, U - B > B - i + 1.2 and 
B — i < 1.5) to our sample. However, our [/-band data, and in the case of the 
faintest candidates also the HST data, is too shallow to get a useful measurement 
on the U — B colour. Instead, we take the best fit spectrum from the SED fitting 
(see Fig. 13. 8p and convolve it with the U (F300W), B and V filter sensitivities and 
calculate the colours. For this model spectrum, these colours become U — B = 4.51, 
B — V = 0.24 and B — i = 0.69 which well satisfy the selection criteria for U- 
band drop-outs, see Fig. 13. Ill However, many of our LEGOs are very faint and the 
stacked B magnitude is fainter than the lower limit of the selection of Madau et al. 
(1996), and about half of our sample are fainter than the V cut-off in the sample of 
Wadadekar et al. (2006), see Fig. [37[3l 

Secondly, we wish to compare the observed optical colours (restframe UV colours) 
of our LEGOs to the LBGs in order to establish if our LEGO candidates have 
the same UV continuum colours as LBGs on the red side of the Lyman break. 
In Fig. I3.13[ we plot the colours of the two samples of faint LBGs published by 
Wadadekar et al. (2006) and Madau et al. (1996) against the colours of our candi- 
dates. All samples are drawn from survey data-sets such as GOODS-S and HDF-N, 
hence there is no bias in photometry. In the plot, we see that the LEGO candidates 
are drawn from a fainter sub-sample of the high-redshift galaxy population. How- 
ever, for the brighter candidates among our sample, the LEGOs appear to have UV 
colours similar to LBG galaxies. 

3.8 Conclusion 

We have performed deep narrow-band imaging of part of the GOODS-S field. The 
image revealed a set of 24 LEGO candidates, at a redshift of 2; ~ 3.15. Of these, three 
candidates have been observed spectroscopically and are confirmed. The spatial 
distribution of the candidates appear to be in a filamentary structure, with a 4a 
confidence, however to confirm this and to plot the filament in 3D-space, we would 
need spectroscopic redshifts. We have studied the entire candidate sample in all 
bands available from X-rays to infrared in the GOODS-S data-set. From the SED 
fitting we conclude that the LEGOs on average have low metallicity {Z/Zq = 0.005), 
have stellar masses in the range of 1 — 5 x 10^ M© and low dust extinction {Ay ~ 0.3). 
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Figure 3.11: Colour-colour plot from Wadadekar et al. (2006) showing simulated 
galaxy colours, see that paper for details. The solid line marks the area (upper 
left corner) where redshift z ^ 3 LBG reside. The large star in the upper left 
corner marks the colours of our best-fit synthetic model, well within the selection 
boundaries for high-redshift LBGs. 



The candidates have ages in the range of 100 - 900 Myrs . We also find one galaxy, 
LEG0_G00DS-S#16, which is best fit by a dusty starburst galaxy a.t z = 3.15 with 
Lya-emission escaping from an area slightly offset from the central core. 

The comparison to a sample of f/-band drop-out galaxies in the GOODS-S field 
show that the colours of LEGOs are consistent with the selection criteria for U- 
band drop-outs except they are too faint to be detected from their continuum flux. 
They also have colours similar to those of LBGs at redshift z ~ 3. In agreement 
with previous results (e.g. Gawiser et al. 2006), we conclude that Lya-emitters at 
redshift z ~ 3.1 are dust- and AGN-free, star-forming galaxies with small to medium 
masses. 



56 



Chapter 3 



GOODS-S 




-4 -3 -2-10 1 

NB - V 



Figure 3.12: Colour-colour plot. Dots mark the whole sample of 2616 objects se- 
lected with SExtractor. Points with errors mark the selected emission-line objects. 
Stars (in magenta) mark the candidate sample, the triangles (in blue) the spectro- 
scopically confirmed LEGOs and the open diamond (in red) the blob (Nilsson et al., 
2006a). Dots detected in the same region of the plot as the selected sample, but 
that are not selected, consist of objects that were discarded in the visual inspection. 




Figure 3.13: Colours of our candidates (stars with error bars) compared to colours 
of the sample of faint LBGs of Wadadekar et al. (2006; diamonds) and Madau et 
al. (1996; triangles). Lines represent the best fit to the LBG data, the solid line 
the fit to the sample by Wadadekar et al. (2006) and the dotted line the fit to the 
Madau et al. (1996) sample. Left V minus i colours. Right V minus z' colours. 
Small arrows indicate upper limits for the Wadadekar et al. (2006) sample. 
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Figure 3.14: Two preliminary GRASIL model fits to the SED of LEGO_GOODS- 
S#16, both with ~ 50. Points with error bars (of the same size as the point) 
are our data, with errors being purely statistical. The point off the curve at lower 
wavelengths is the narrow-band detection. 
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hya emitters with VISTA 



4.1 Introduction 

In this chapter, I describe a project which started in February 2005. The idea was to 
make use of the large field-of-view of the camera of a future telescope called VISTA 
(Visible and Infrared Survey Telescope for Astronomy), described further in the 
next section, to find very high redshift {z ~ 9) Lya emitters through narrow-band 
imaging in the near-infrared. The cost of such filters made this a rather expensive 
project and funding for the filters very generously came from two sources; IDA - 
Instrument center for Dansk Astrofysik (DK) and the Dark Cosmology Centre (DK). 
In the spring of 2006, we extended the science case to involve all emission-lines (Ha, 
[OIII], H/5 and [Oil]) and eventually submitted a proposal for a Public Survey 
with VISTA, entitled "ELVIS - Emission Line galaxies with VISTA Survey". This 
proposal was well received by the ESO Public Survey Panel, but was asked to merge 
with a number of other surveys, with the aim to observe small, but deep, areas of 
the sky with broad-band filters. Thus, ELVIS was merged into the Ultra- VISTA 
survey. I have been involved in all steps of this project, from the original idea to 
ordering/inspecting the filters, and from proposal writing to the accepted Public 
Survey. The following sections include a short introduction to VISTA, a section on 
the design and procuration of the filters and a description of the ELVIS Survey. This 
project also highlighted the need for better estimates of number densities of very 
high redshift Lya emitters. The results of such a project are presented in Chapter O 

4.2 VISTA - Visible and Infrared Survey Tele- 
scope for Astronomy 

New advances in technology will always prompt new surveys of the night-time sky. 
The first large area sky survey in the near infrared, the Two Micron Sky Survey 
(TMSS; Neugebauer & Leighton 1969), was carried out in the 1960's and the next 
large area survey did not happen until the late 1990's with the 2MASS Survey 
(Skrutskie et al. 2006) covering the entire sky in the three main bands J, H and 
Kg. A survey called UKIDSS was started in 2005 (Lawrence et al. 2007). It has 
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Table 4.1: Technical details of VISTA. See also Emerson et al. (2004) and Dalton 
et al. (2006). 



Ml diameter 


3.95 m. 


1\ /TO 1 * J. 

iV12 diameter 


"1 (~\ A 

1.24 m. 


/ ratio at instrument 


3.25 


Mount 


Altitude- Azimuth 


Number of detectors 


16 


Type of detector 


Raytheon VIRGO HgCdTe 0.84-2.5 micron 


Number of pixels per detector 


2048 X 2048 


Pixel size 


0.34"/pixel 


Read-out noise 


20.9 e" 


Detector Quantum Efficiency 


71 % (J), 74 % (H), 75 % (K,) 


Operating temperature of camera 


~ 80 K 


Available filters 


z, Y, J, H, K„ NB1185 



several parts, both wide/shallow and small/deep fields, all in the Northern Sky. It 
is on the background of these surveys that VISTA was conceived. VISTA, originally 
planned to have both an infrared and an optical camera but later designed with 
only an infrared camera, is an almost (at the time of writing this thesis) completed 
4-m. telescope dedicated almost entirely to near-infrared surveys. For this purpose, 
it has been equipped with an infrared camera array of 16 detectors, covering a field- 
of-view of 0.6 square degrees in each single exposure, see Fig. 14.11 The detectors 
are not buttable, i.e. they cannot be placed adjacent to each other, but need to be 
separated by a large fraction of the width of the detector itself. This is because the 
detectors interfere with each other if they are too near each other. Thus, to get a 
continuous mosaic image, the camera needs to be shifted in six different positions 
(or paw-print, as each single exposure is called). A plot of the resulting exposure 
time coverage can be seen in the right panel of Fig. 14.21 

The field-of-view (FOV) of each paw-print with VISTA is unprecedented as can 
be seen in the FOV comparison in the left panel of Fig. 14. 2[ The image shows the 
FOVs of current detectors (HST/NICMOS, VLT/ISAAC and UKIRT/WFCAM) 
and a future detector (VLT/HAWK-I). As in the picture, the VISTA FOV is ap- 
proximately 3 times the FOV of WFCAM (only northern hemisphere) and 346 times 
the FOV of ISAAC. 

A schematic lay-out of the telescope and camera can be seen in Fig. 14.31 and the 
technical details of VISTA are given in Table I^TTl A point to note is that the VISTA 
detector array is stationary in the camera, which has no shutter, and the filters are 
turned into the light path with a filter wheel. The filter wheel has eight slots, of 
which one holds a blank plate, six hold the filters mentioned in Table 14.11 and the 
final slot is empty initially. The fact that there is no shutter makes positioning 
of filter sets and rotation patterns of the filter wheel crucial. When moving to 
a blue filter, or a narrow-band filter, it is undesirable to turn the wheel so that 
the detector array is exposed to the light coming through one of the redder filters. 
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Figure 4.1: VISTA IR camera array consisting of 16 detectors. The single "paw- 
print" or exposure covers a field-of-view of 0.6 sq.degrees. The array has large gaps 
between the detectors and hence needs to be moved around to fill in a full mosaic. 
From www.vista.ac.uk. 



as this will "flash" the array with light and reduce the sensitivity temporarily. 
VISTA is funded by PPARC, UK, but operated by ESO. It is placed at the "NTT" 
peak at Paranal, Chile and will be fully operational by early 2008, according to the 
plan. The majority, 75%, of VISTAs time is dedicated to Public Surveys, i.e. large 
scale surveys conducted by the science community where all data, raw and reduced, 
become public. Six surveys have been approved for the initial five year survey period. 
An overview of each survey is given in Table 14. 2[ 

4.3 The narrow-band filters 

The width of a VISTA narrow-band filter can not be narrower than AA ^ 120 A due 
to the large shifts in wavelengths over the field-of-view, see also sec. I4.3.2[ As the 
sky spectrum in the J band is full of atmospheric OH-lines, there are only a small 
number of possible wavelengths to place such a filter. These correspond to e.g. 
Zhya = 7.73, 8.22, 8.78 etc. The first redshift, and later on the last, are planned 
to be surveyed with DAzLE (The Dark Ages z Lyman-alpha Explorer, Horton et 
al. 2004). Our first idea was thus to place the filter at the intermediate redshift. 
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Figure 4.2: Left Field-of-view comparison between VISTA and other current or near 
future IR arrays as well as the size of the moon. Right Coverage map of the VISTA 
IR array after off-setting the array in six different positions to get a uniform coverage 
of the sky. The colour scheme shows effective exposure time, where dark green is 
one time the single paw-print exposure time, light green is two times, magenta is 
three times, red is four times and yellow is six times the single paw-print exposure 
time. The average, almost homogeneous, coverage is two times the single paw-print 
exposure time. Both images from www.vista.ac.uk. 




Figure 4.3: Schematic view of the VISTA IR telescope (left) and camera (right). 
From www.vista.ac.uk. 
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Table 4.2: Approved Public Surveys for the first five years of operation of VISTA. 
See also Arnaboldi et al. 2007. 



Name 


Area (deg^) 


Filters 


Science case 


ultia- V io ii\ 


n Q 
u.y 


1, J, n, t\gj iNiJiioo 


formation anH pvolntion 

very deep, small area survey 


VHS 


20000 


Y, J, H, K, 


Complete, shallow, southern 
hemisphere map 


VIDEO 


15 


z, Y, J, H, K, 


AGN surveys, galaxy clusters, 
very massive galaxies 


VVV 


520 


z, Y, J, H, K, 


Galactic bulge and plane, 
open and globular clusters, 
variable sources 


VIKING 


1500 


z, Y, J, H, K, 


Weak lensing and baryon 
acoustic oscillations 


VMC 


184 


Y, J, K, 


Survey of the Magellanic system 



However, it was discovered that this window has a significant [OII]-line in the dark 
range, and it was discarded as an option. The window corresponding to zi,ya = 8.8 
is quite wide (see Fig. 14.41) and it will only be surveyed by DAzLE at a later stage, 
hence we chose to place our filter at the wavelength corresponding to ZLyo = 8.78, 
Az = 0.10. 

4.3.1 Filter specifications 

The filter specifications were prepared in the same manner as the original VISTA 
broad-band filters specifications. The camera array will experience a slight passband 
shift over the field-of-view due to the filters operating in an f/3.25 beam, and the 
outer arrays may experience a broadening of ~ 0.5 % due to a slight field-dependent 
effect, due to the chief ray being tilted by ~ 5 degrees at the edge of the field. 
These effects will be especially pronounced with the narrow-band filter, due to the 
small passband width. We investigated how these effects would influence the filter 
specifications thoroughly prior to sending out calls for quotations, and found that 
the effect would not be severe, see sec. 14.3.21 A plot of such a study can be seen 
in Fig. 14. 6[ We also anticipated the possibility of extensively measuring the effects 
during commissioning of the telescope and following that produce corrective algo- 
rithms to be used while reducing the data. The specifications made for the filter 
were thus set to be: 

- Central wavelength should be 1185 nm ± 2 nm. 

- Peak transmission should be at least 75% (Goal: 85%). 
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Figure 4.4: Plot of the sky background OH-lines, with the intended filter curve 
over-plotted with a dashed line in arbitrary units. 



- The passband FWHM should be 10 nm ± 2 nm. 

- The transmission at 1174 nm and 1195 nm should be below 10%. 

- The transmission at 1165 nm and 1210 nm should be below 1%. 

- The average transmission between 700 nm and 1140 nm should be below 0.1%. 

- The average transmission between 1250 nm and 3000 nm should be below 
0.01%. 

The last point is especially important as "red leaks" are a severe problem with 
near-infrared narrow-band filters (Sutherland, priv. communication). The call for 
quotations was sent to four firms in October 2005, and the best bid was made by 
NDC Infrared Engineering^ The filters were subsequently ordered by the Dark 
Cosmology Centre in January 2006 and 20 filters were delivered to the Rutherford 
Appleton Laboratory in Oxford, UK, in April 2007, see Table 14.31 



http: / / www.ndcinfrarcd.com/ 
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Table 4.3: Table of delivered narrow-band filters from NDC Values given at operat- 
ing temperature. The filter marked with a star has been excluded from the selection 
due to a large gradient in the color on the edge of the filter. 



Filter 


CWL (nm) 


FWHM (nm) 


TV T * 1 m / fry \ 

Maximal T (%) 


911056R5_4 


1182 


11.1 


71 


301106R6_1 


1183 


11.6 


72 


011206R5_1-1 


1183 


12.0 


71 


021206R5_3 


1184 


11.0 


73 


011206R5_1 


1185 


11.2 


70 


260207R5_1 


1185 


12.0 


76 


281106R5_3 


1185 


10.7 


70 


230207R6_3 


1185 


11.9 


74 


UoUzU ( i\0_i 


iioo 


lis 
ii.o 


TA 

/4 


230207R5_3 


1186 


11.5 


78 


270207R5_4 


1187 


12.0 


70 


260207R5_3 


1187 


12.0 


74 


090207R5_4 


1188 


11.6 


73 


230207R5_5 


1188 


11.9 


76 


050207R6_1 


1188 


11.7 


73 


090207R5_1 


1188 


11.8 


71 


230207R5_2'^ 


1188 


11.8 


73 


230207R6_2 


1188 


11.8 


76 


260207R5_4 


1188 


11.9 


77 


311006R5_1 


1189 


10.5 


76 



4.3.2 Central wavelength and passband shift 

The optical path of VISTA is a converging beam. This means that each pixel on 
each detector sees an annulus of light coming from the primary mirror and obscured 
by the secondary mirror. But a pixel at the edge of the filter/detector array has a 
tilted view of this annulus, i.e. the centre of the annulus is shifted into the actual 
light annulus. The effect is a cosine function of the length of a vector pointing out 
a light element in the annulus divided by the refractive index (assumed to be 1.8, 
Sutherland priv. communication), integrated over the entire annulus. This causes 
the central wavelength of the filter to shift and the passband to increase. 

A schematic of the annulus can be seen in Fig. 14. 5[ The shifted central wave- 
length is given by: 

(s"-Z.. f::% "»(^)<ixdy - f::x s::%, »s(4?!)dxdy 

\ _ _\ 

/^shifted — ^ V r fi>2 p2 N 

where 



(4.1) 



/(x) = 



out 



X 



2 
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fix) = ^Jri - X2 

where Rout = 8.75 degrees, Rm = 3.85 degrees and Ac = 1185 nm (Sutherland, 
priv. communication). This is true for a central pixel at the very centre of the 
array. For any other pixel, a value of xq and yo has to be added to x and y in the 
cosine function in the top integral. The largest effect is achieved when (a;o,?/o) = 
(5.15, 5.15), which is at the corner of the field. The effect of the integrated factor is 
always to decrease the central wavelength and the effect is at its greatest a factor of 
0.5 %. The change in wavelength for ziya = 8.8 can be found in Table S3] for four 
positions. These positions correspond to e.g. the positions of detector 6 (number 1 
in table), detector 5 (number 2 in table), detector 1 (number 3 in table) and detector 
2 (number 4 in table). The shifts are symmetric in each quadrant of the array. 

Table 4.4: Central wavelength shift in the centre position and the middle point of 
four detectors according to the text. Based on an original central wavelength of the 
filter of 1185 nm. 



Shifted wavelength in four points 


Point 


Ashifted (nm) 


Centre 


1182.46 


1 


1182.34 


2 


1182.01 


3 


1181.43 


4 


1181.76 



We assume that the filter transmission function is a Gaussian function: 

1 (X-Xc)^ 

T = — = X e 5^ (4.2) 
(TV 27r 



where FWHM = 2 x V 21n2o" and Ac = central wavelength of the filter. Since 
the central wavelength shifts over each camera array, the transmission will become 
broadened and skewed towards a top-hat function. Figure 14.61 show this effect on 
the corner array in the z = 8.8 OH-line free window. The effect is a broadening in 
FWHM of the order of 10 %. 



4.3.3 Inspection of narrow-band filters 

In May, 2007, we visited Rutherford Appleton Laboratory (RAL) in order to inspect 
the filters. Some photographs from the inspection can be found in Fig. 14.81 

During the inspection, performed in the clean room at RAL, we carefully opened 
each package and checked each filter for the following issues: 

- Consistency in labelling on wrapping and filter. 

- Overall condition of filter. 
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Figure 4.5: Sketch of the annulus with inner and outer radii Rin and Rout respec- 
tively. To estimate the shift in central wavelength, we need to integrate the incident 
light on a pixel, by integrating the area of the annulus for all x and y. The centre 
of the annulus then moves with increasing distance from the centre of the camera 
array. The shift is given by Xq and ?/o- Fot an edge or corner pixel we hence need to 
integrate over the annulus, with shifted central pixels. 



0.01 5 
0.01 



0.005 



0.000 



1 .1 6x1 0*1 .1 7x1 0*1 .1 8x1 0*1 .1 9x1 0*1 .20x1 0*1 .21 x1 0^ 
Wavelength (AA) 

Figure 4.6: Plot of the sky background and the filter specifications given to the 
productions company. The dot-dashed line indicates the desired filter transmission 
curve, the dotted and dashed vertical and horizontal lines are the specifications 
made to the company (see sec. 14.3.11) . The solid line is the worst case shift at the 
corner of the camera array. As can be seen, this effect will not significantly overlap 
with the OH lines near the edge. 
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- Level of flaking. 

- Existence of "pinholes" , "comas" , scratches or other damages. 

- Discolorations. 

All filters displayed varying degrees of flaking, i.e. material from the coating at 
the edges falling off and lying loose in the package. This is very undesirable, as it 
may scratch the coating if it rubs against the surface. Two filters were mislabeled, 
or even unlabeled on the filter itself. One filter showed a large discoloration along 
the edges, indicating that the coating was incomplete around the edges. This filter 
was removed from the pool of usable filters. A few filters had small chips at the 
edges, although small enough to not cause any complications. Many filters had small 
pinholes but they were shallow enough that they can be fiat-fielded away (Dalton, 
priv. communication). Overall, all but the one with a discolored edge were accepted 
for use. 

4.3.4 Positioning of filters in the VISTA filter tray 

After studying the transmission curves of the remaining, accepted 19 filters, the 
three most blue filters were also excluded from further use, and kept as spares, as 
they overlap significantly with sky emission lines in the blue edge of the sky window, 
as seen in Fig. 14. 4[ Thus, with 16 filters left to use we attempted to "puzzle" these 
filters together so that all detectors along each column of the detector array would 
have as similar transmission curves as possible, also taking into account the shift 
that occurs of the field-of-view. The result can be seen in Fig. 14. 7[ 

4.4 ELVIS and Ultra- VISTA 

ELVIS is part of Ultra- VISTA, an ESO Public Survey that aims to do very deep 
near-infrared broad- and narrow-band imaging in one pointing in the COSMOS 
field (Scoville et al. 2006). The main science case for Ultra- VISTA is to get one 
of the deepest images to date in the near-IR bands of a large area on the sky. It 
will focus on detecting the first galaxies at very high redshift and to study faint 
objects at intermediate redshifts. ELVIS aims specifically at searching for emission- 
line galaxies at several redshifts in order to study galaxy formation and evolution, 
at what redshift re-ionisation started as well as the evolution of the cosmic star 
formation history. 

4.4.1 Science goals 

As my involvement in Ultra- VISTA is almost solely focused on the ELVIS part, I 
will here only discuss the science goals of ELVIS. There are a number of arguments 
why a large scale survey in the infrared for emission line galaxies will yield new and 
exciting insights into galaxy formation and star formation history. The objective 
is split into three parts, one looking for very high redshift Lya emitters, and the 
others to observe Ha and intermediate redshift emitters. 
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Figure 4.7: Positioning of VISTA narrow-band filters. Upper left image shows col- 
umn 1 (detector 1-4), upper right shows column 2 (detector 5-8), lower left shows 
column 3 (detector 9 - 12) and lower right shows column 4 (detector 13 - 16). Red 
lines display the OH sky background, differently coloured lines display the transmis- 
sion curves of each individual filter. Green and magenta lines represent filters to be 
placed at the uppermost and lowermost position in the tray (e.g. detector 1 and 4 in 
column 1) and blue and yellow lines show filters to be placed in the central locations 
(e.g. detector 2 and 3 in column 1). Curves have been shifted to correspond to the 
position they are to be placed in. 




Figure 4.8: Photographs from the filter inspection. To the left are the whole stack 
of filters, including two witness pieces, delivered from NDC Infrared to Rutherford 
Appleton Laboratory. The middle panel show the blocking side of the filter, in 
relation to the size of a finger. The right panel show the filter transmission side of 
the filter. Images courtesy of Wolfram Freudling. 
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LycK emitters 

The quest for detecting cosmological objects is always pushing to reach further away 
in redshift space, and further back in time. So far, no objects have been confirmed 
to be at a redshift larger than seven. Of the objects confirmed to be above redshift 
six, the vast majority are in fact Lya-emitters. To date, 22 Lya-emitters (Hu et 
al. 2002; Kodaira et al. 2003; Rhoads et al. 2004; Kurk et al. 2004; Stern et al. 
2005a; Taniguchi et al. 2005 and Kashikawa et al. 2006), 13 quasars (Fan et al. 
2006; Willott et al. 2007) and one GRB (Haislip et al. 2006) with spectroscopic 
redshift above six has been published. Another two detections of gravitationally 
lenscd Lja emitters at redshifts around z ^ 9 have been suggested in Stark et al. 
(2007). The reason for the easy detection of Lya-emitters is the relative ease with 
which they are found (narrow-band imaging) and with which the redshift can be 
confirmed (spectroscopic follow-up and detection of single emission line) . Hence, it 
is often believed that Lya-emitters will be the most successful tool in discovering 
even higher redshift sources. Some projects have already started. Three attempts 
to observe Lya-emitters at redshift 2; ~ 9 have been made (Parkes, Collins & Joseph 
1994; WiUis & Courbin 2005; Cuby et al. 2006) but the results are null-detections 
due to too small observed fields and/or to shallow flux hmits. The DAzLE project 
has been presented by Horton et al. (2004). It is designed to find Lya-cmittcrs at 
redshift z = 7.73 and z = 8.78, using a specially made instrument to be placed in 
visitor focus at the VLT. The DAzLE instrument, however, has a small field of view 
(6.83'x6.83'). As these sources will be very under-luminous due to the distance to 
them, and we thus sample only the top of the Lya luminosity function, the key to 
finding a large sample of very high redshift galaxies is large survey area and volume. 
This is the main strength of the VISTA telescope. 

There arc two main reasons why we wish to observe very high redshift Lya 
emitters; re-ionisation and galaxy formation studies. Re-ionisation is the epoch 
of the history of our Universe when its vast amounts of hydrogen gas was ionised 
after being completely neutral. This happened when the first stars, and potentially 
the first quasars, lit and galaxies started to form. But it is yet unknown exactly 
when this event took place. Recent results from studying the cosmic microwave 
background (CMB) have constrained the upper limit of the re-ionisation redshift to 
be Zj-eion 10 and observations of galaxies and quasars have set the lower limit to 
be Zreion > 6.5 (c.g. Spcrgcl et al. 2006; Malhotra & Rhoads 2004; Fan et al. 2006). 
Several authors have argued that Lya-emitting high-redshift galaxies may provide 
a tool to constrain the redshift of re-ionisation better (e.g. Miralda-Escudc 1998; 
Miralda-Escude & Rees 1998; Haiman 2002; Gnedin & Prada 2004; Haiman & Cen 
2005). These authors show that absorption in the Gunn- Peterson trough (i.e. the 
absorption on the short wavelength side of the Lya line due to the intergalactic 
medium; Gunn & Peterson 1965) may extend to the red side of the Lya emission 
line and cause damping wings. Haiman (2002) and Haiman & Cen (2005) suggest 
that this damping will change the shape of the luminosity function of Lya-emitters 
before and after re-ionisation. McQuinn et al. (2007) suggest that the clustering of 
very high redshift Lya emitters is an independent and powerful tool to diagnose the 
level of re-ionisation. With a large sample of Lya-emitters at redshift z — 8.8, it 
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will be possible to analyse the Lya emission line, make a comparison of luminosity 
functions of Lya-emitters at different redshifts, and determine the level of clustering, 
thus revealing when re-ionisation happened, or at a minimum placing constraints 
on the time of re-ionisation. 

When it comes to galaxy formation and evolution, possible observables include 
number density of Lya-emitters, SFR and psFR and clustering effects. An inter- 
esting question to answer is how many Lya-emitters existed at a specific time? 
Taniguchi et al. (2005) study the number densities of several Lya-surveys between 
z ~ 3.4 to 6.6 and find no evolution with comoving number density. Another key 
question related to galaxy formation and evolution is that of the shape of the star 
formation history. Observations at lower redshifts indicate that Lya-emitters are 
moderately star-forming, dust free galaxies with little or no AGN content. It is thus 
of interest to investigate if the SFR was different in this class of objects at very 
high redshift, what percentage of the total star formation happened in this type of 
galaxy at very high redshift and what the star formation density was at this epoch. 

Emission-galaxies at intermediate redshift 

At lower redshift, ELVIS will also be sensitive to emission hne galaxies with much 
higher surface density; Ha ai z — 0.80, [OIII] z — 1.36, H/3 at z — 1.43, and 
[OH] at z = 2.17. These redshifts span most of the time when the cosmic star 
formation occurred, associated with the formation of the stars and galaxies we see 
in the Universe today. The [OIII]-, H/3-, and [OII]-lines can all be used as tracers of 
star formation, and [OIII] and [Oil] are also frequently used as tracers of AGN and 
Seyfert galaxies. The line profiles and ratios of these emission fines give information 
on the kinematics of the emitting gas, and of the properties of the AGN (Heckman et 
al. 1981; Boroson 2005; Gu et al. 2006). The forbidden oxygen lines are metallicity 
dependent, but also affected by AGN. Nevertheless, in particular [Oil] is still a 
good tracer of star formation and hence we will have an interesting handle on the 
star formation density at z — 2.2, which is complementary to broad-band surveys 
targeting similar redshifts (e.g. Adelberger et al. 2004). 

Ha emitters 

Ha is a one of the best, direct tracers of the instantaneous star formation rate 
and it is particularly useful as it is relatively unaffected by metallicity and dust 
extinction. Several authors have attempted to obtain the Ha luminosity function 
at redshifts close to z = 1, which is believed to be the peak of star formation in the 
Universe, but sample sizes are still small (e.g. Tresse et al. 2002; Doherty et al. 
2006; Hopkins et al. 2000; Yan et al. 1999). Evidence exists for strong evolution 
in the SFR from redshift zero to one, by perhaps an order of magnitude (Hopkins 
2004, and references therein). However, many surveys are limited to the bright end 
of the luminosity function due to the small sample sizes and it is therefore not clear 
whether this evolution in global SFR is a property of galaxies with well above average 
SFRs, or whether it also extends to average SFR galaxies. The Ha luminosity of 
normal galaxies is well studied in the local universe. Comparing the full range of 
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Table 4.5: Expected flux depths from the Ultra- VISTA survey. 





Deep Survey 


ELVIS 


Ultra Deep Survey 


Filter 


Y J H K, 


NB1185 


Y J H K, 


Total exp. time (h) 


48 48 48 48 


180 


320 320 320 320 


5(7 depth, AB 


25.7 25.5 25.1 24.5 


24.1 


26.7 26.6 26.1 25.6 



Ha luminosities at higher redshifts to those of local galaxies could therefore provide 
crucial insight into the evolution of the SFR. Another question recently posed is that 
of the apparent "downsizing" of star forming galaxies (Cowie et al. 1996; Heavens et 
al. 2004; Juneau et al. 2005), claiming that more massive galaxies form earlier than 
less massive galaxies. To date, no survey has simultaneously reached the sensitivity 
and the area necessary to fully address the questions of global SFR and "down- 
sizing" . The sensitivity of ELVIS will reach Ha luminosities more than three orders 
of magnitude below L*(Ha). This will enable unprecedented constraints to the Ha 
luminosity function ai z = 0.8, in particular the value of the faint-end slope which is 
still a matter of debate, hence giving a direct observational evidence as to whether 
the "down-sizing" scenario is true or merely an observational bias. 

4.4.2 Survey plan 

Ultra- VISTA has three parts; a "deep" survey, a "very deep" survey and ELVIS. 
All parts will observe the COSMOS field (R.A. = 10''00™28.6^ Dec = +02°12'21" 
(J2000)). The VISTA camera array consists of 16 detectors, placed in a square 
pattern with large gaps in between, see Fig. 14. 1[ The gaps can be filled in by 
shifting the array in a six step pattern. A full mosaic will then be 1.6 sq. degrees. 
The broad-band components of Ultra- VISTA will observe a full mosaic and four 
strips respectively. The "deep" survey will observe a full mosaic but the "very 
deep" survey will only shift the detectors along the y-axis, creating four strips. 

For ELVIS, it was debated for some time which configuration to choose; the 
continuous mosaic or the four strips. The total integration time was assumed to be 
fixed and so the full mosaic would be more shallow and the configuration with four 
strips would go deeper. Based on the expected number counts and clustering results 
(see Chapter [5]) it was decided that ELVIS would also observe in four strips but to 
a deeper flux limit. The expected depths reached by the different Ultra- VISTA 
surveys can be seen in Table 14.51 

The survey is planned to start early 2008, with a first release to the public about 
six months later, followed by yearly releases until the survey is completed. If the 
telescope operates at peak efficiency and Ultra- VISTA is awarded all the time in its 
RA range, the program will be completed in five years. 

4.4.3 Expectations 

The expectations, i.e. the expected number of objects detected, in ELVIS is divided 
into two parts; Lya-emitters, and lower redshift emitters. 
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Lya emitters 

The attempt to predict the number of Lya emitters detected with ELVIS sparked 
the project presented in Chapter [5l which presents results from two theoretical 
models of very high redshift Lya emitters; one semi-analytic model based on ACDM 
and one phenomenological model. The semi-analytic model GALFORM has been 
presented in Le Delliou et al. (2005; 2006) and follows a hierarchical evolution 
of galaxies. It has been shown to well reproduce the Lya luminosity functions 
at lower redshifts. The phenomenological model has been presented in Thommes 
Sz Meisenheimer (2005) and assumes that Lya luminosity is proportional to star 
formation rate, which in turn is proportional to the baryonic mass of the galaxy. 
The results of comparing the two models and extrapolating results at higher redshifts 
are described in Chapter [5l and for ELVIS the models predict the detection of 3 - 
13 Lya emitters. 

Lower redshift emitters 

To estimate the number of Ha-emitters that may be found, we use the Ha luminosity 
function of Tresse et al. (2002) with Hq = 73, Qm = 0.3, fl\ = 0.7 and without 
reddening correction (Doherty et al. 2006). Using this luminosity function in a field 
of view of 0.9 deg^, with a line flux limit of 3.7 x 10~^^ erg cm~^ s~^ and over a 
redshift range z = 0.802 — 0.820 we calculate that of the order 3500 Ha-emitters 
would be detectable. The expected density of intermediate redshift galaxies will 
probably be smaller than the density of Ha emitters because of the higher distance 
moduli and (on average) lower equivalent widths, even though the higher redshift 
increases the observed equivalent widths. We estimate that it will be possible to 
detect about an order of magnitude fewer such objects compared to Ha emitters, 
hence of the order of 350 intermediate redshift galaxies. 
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Predicting results from very high 
redshift Lya surveys 

This paper has been accepted for pubhcation in Astronomy & Astrophysics on 31 
August, 2007. The authors are Nilsson, K.K., Orsi, A., Lacey, C.G., Baugh, CM., 
& Thommes, E. 

5.1 Abstract 

Context Many current and future surveys aim to detect the highest redshift {z > 7) 
sources through their Lyman-Lya (Lya) emission, using the narrow-band imaging 
method. However, to date the surveys have only yielded non-detections and upper 
limits as no survey has reached the necessary combination of depth and area to 
detect these very young star forming galaxies. 

Aims We aim to calculate model luminosity functions and mock surveys of Lyo; 
emitters ai z >1 based on a variety of approaches calibrated and tested on obser- 
vational data at lower redshifts. 

Methods We calculate model luminosity functions at different redshifts based on 
three different approaches: a semi-analytical model based on CDM, a simple phe- 
nomenological model, and an extrapolation of observed Schechter functions at lower 
redshifts. The results of the first two models are compared with observations made 
at redshifts z ~ 5.7 and z ~ 6.5, and they are then extrapolated to higher redshift. 
Results We present model luminosity functions for redshifts between z = 7 — 12.5 
and give specific number predictions for future planned or possible narrow-band 
surveys for Lya emitters. We also investigate what constraints future observations 
will be able to place on the Lya luminosity function at very high redshift. 
Conclusion It should be possible to observe z = 7 — 10 Lya emitters with present 
or near-future instruments if enough observing time is allocated. In particular, 
large area surveys such as ELVIS (Emission Line galaxies with VISTA Survey) will 
be useful in collecting a large sample. However, to get a large enough sample to 
constrain well the z > 10 Lya luminosity function, instruments further in the future, 
such as an ELT, will be necessary. 

7K 
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5.2 Introduction 

One of the most promising ways of detecting very high redshift {z > 5), star- 
forming galaxies is via narrow-band imaging surveys targeting Lyman-a (Lya). In 
particular, redshifts z ~ 5.7 and 6.5 have been extensively surveyed by several 
groups (e.g. Ajiki et al. 2003; Hu et al. 2004; Shimasaku et al. 2005; Ouchi et 
al. 2005, 2007; Malhotra et al. 2005; Taniguchi et al. 2005; Tapken et al. 2006; 
Kashikawa et al. 2006). The current redshift record for a spectroscopically confirmed 
Lya emitter (LEGO - Lya Emitting Galaxy-building Object; see M0ller & Fynbo 
2001) is z = 6.96 (lye et al. 2006) although Stark et al. (2007) have suggested 
the discovery of two LEGOs at z = 8.99 and 9.32. The reason why narrow-band 
surveys are restricted to a discrete number of narrow redshift windows is the night 
sky OH emission fines. According to the OH fine atlas of Rousselot et al. (2000), 
at Lya redshifts ZLya ^ 7 (A > 9800 A) there are only a few possible wavelengths 
where a narrow-band filter can fit in between the OH sky lines. These correspond 
to ZLya ~ 7.7, 8.2, 8.8, 9.4 and 10.1 — 10.5. Several future surveys will target these 
windows in the sky aiming to detect very high redshift galaxies. Three narrow-band 
surveys for Lya at redshift z ~ 8.8 have already been completed (Parkes, Collins 
& Joseph 1994; WiUis & Gourbin 2005; Guby et al. 2007) but have only yielded 
upper limits. Future surveys planned for these redshifts include DaZle (Dark Ages 
Z Lyman-a Explorer, Horton et al. 2004) and ELVIS (Emission-Line galaxies with 
VISTA Survey, Nilsson et al. 2006b). Observations of very high redshift LEGOs 
have been proposed as an excellent probe of reionisation, through its effects on the 
Lya emission line profile (e.g. Miralda-Escude 1998; Miralda-Escude & Rees 1998; 
Haiman 2002; Gnedin & Prada 2004), the luminosity function (e.g. Haiman & Gen 
2005; Dijkstra, Wyithe & Haiman 2007b) and the clustering of sources (McQuinn 
et al. 2007). 

We here focus on Lya emission from star-forming galaxies, where the Lya pho- 
tons are emitted from gas which is photo-ionised by massive young stars. During re- 
cent years, theoretical work on Lya emitting galaxies has made significant progress. 
There are three main aspects to these studies: i) predicting the numbers of star- 
forming galaxies as a function of star formation rate and redshift, ii) calculating 
the fraction of the Lya photons which escape from galaxies into the IGM and Hi) 
calculating the factor by which the Lya fiux is attenuated by scattering in the IGM 
on its way to the observer. Accurate treatments of ii) and Hi) are complicated 
because Lya photons are resonantly scattered by hydrogen atoms, with the con- 
sequences that absorption of Lya by dust in galaxies is hugely amplified, thereby 
reducing the escape fraction, and that even a small neutral fraction in the IGM can 
be effective at scattering Lya photons out of the line-of-sight, thus attenuating the 
flux. Because of these compfications, most theoretical papers have chosen to con- 
centrate on only one aspect, adopting simplified treatments of the other two aspects. 
Haiman & Spaans (1999) made predictions of the number counts of Lya emitting 
galaxies by combining the Press-Schechter formalism with a treatment of the inho- 
mogeneous dust distribution inside galaxies. Barton et al. (2004) and Furlanetto et 
al. (2005) calculated the numbers of Lya emitters in cosmological hydrodynamical 
simulations of galaxy formation, but did not directly calculate the radiative trans- 
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fer of Lya photons. Radiative transfer calculations of the escape of Lya photons 
from galaxies include those of Zheng & Miralda-Escude (2002), Ahn (2004) and Ver- 
hamme, Schaerer & Maselli (2006) for idealised geometries, and Tasitsiomi (2006) 
and Laursen & Sommer-Larsen (2007) for galaxies in cosmological hydrodynamical 
simulations. The transmission of Lya through the IGM has been investigated by 
Miralda-Escude (1998), Haiman (2002), Santos (2004) and Dijkstra, Lidz & Wyithe 
(2007a), among others. Several authors (e.g. Haiman, Spaans & Quataert 2000; 
Fardal et al. 2001; Furlanetto et al. 2005) have studied the effect of cold accretion 
to describe the nature of so-called Lya blobs (Steidel et al. 2000; Matsuda et al. 
2004; Nilsson et al. 2006a), see also sec. K7\ 

Two models in particular, dissimilar in their physical assumptions, have been 
shown to be successful in reproducing the observed number counts and luminosity 
functions of Lya emitting galaxies at high redshifts: firstly, the phenomenological 
model of Thommes & Meisenheimer (2005) which assumes that Lya emitters are 
associated with the formation phase of galaxy spheroids, and secondly the semi- 
analytical model GALFORM (Cole et al. 2000, Baugh et al. 2005), which follows the 
growth of structures in a hierarchical, ACDM scenario. The GALFORM predictions 
for Lya emitters are described in detail Le Delliou et al. (2005, 2006) and Orsi et al. 
(in prep.), who show that the model is successful in reproducing both the luminosity 
functions of Lya emitting galaxies in the range 3 < 2; < 6 and also their clustering 
properties. 

In this paper we aim to provide model predictions to help guide the design of 
future planned or possible narrow-band surveys for very high redshift Lya emit- 
ters. We make predictions based on three approaches: the semi-analytical and phe- 
nomenological models already mentioned, and an extrapolation from observations 
at lower redshift. In section [5?3] we describe the different models used to make the 
predictions, and in section [57il we present the predicted number counts and compar- 
isons with observed luminosity functions at lower redshifts. In section 15.51 we make 
number predictions for some specific future surveys. A brief discussion regarding 
what can be learned from these future surveys is found in section 15.61 We give our 
conclusions in section 15.71 

Throughout this paper, we assume a cosmology with Hq = 70 km s~^ Mpc"^, 
flrn = 0.3 and = 0.7, apart from the mock surveys discussed in section which 
use GALFORM models matched to the cosmology of the Millenium Run (Springel 
et al. 2005), (which has Hq = 73 km s"^ Mpc'^, = 0.25 and = 0.75). 

5.3 Models 

We use three different approaches to predict the numbers of high redshift {z > 7) 
Lya emitters. The models are based on very disparate assumptions. The first model 
is the semi-analytical model GALFORM (Le Delliou et al. 2005, 2006), the second 
is the phenomenological model of Thommes & Meisenheimer (2005), and the third 
model is based on directly extrapolating from observational data at lower redshifts. 
Both the semi-analytical and phenomenological models assume that the fraction 
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of Lya photons escaping from galaxies is constant, and that the IGM is transparent 
to Lja. The simple expectation is that before reionisation, the IGM will be highly 
opaque to Lya, and after reionisation it will be mostly transparent. However, var- 
ious effects can modify this simple behaviour; e.g. Santos (2004) finds that the 
transmitted fraction could be significant even before reionisation, while Dijkstra, 
Lidz & Wyithe (2007a) argue that attenuation could be important even after most 
of the IGM has been reionised. The WMAP 3-year data on the polarisation of the 
microwave background imply that reionisation occurred in the range z ^ 8 — 15 
(Spergel et al. 2007), i.e. the IGM may be mostly transparent to Lyct at the red- 
shifts of most interest in this paper. In any case, what is important for predicting 
fluxes of Lya emitters is the product of the escape fraction from galaxies with the 
attenuation by the IGM. The two effects are in this respect degenerate. 

5.3.1 Semi-analytical model 

The semi-analytical model GALFORM (Cole et al. 2000; Baugh et al. 2005), 
which is based on ACDM, has been shown to be successful in reproducing a range 
of galaxy properties at both high and low redshift, including Lya emitters in the 
range z = 3 - 6 (Le Delliou et al. 2005; 2006). A full description of GALFORM 
is given in these earlier papers, so we only give a brief summary here. GALFORM 
calculates the build-up of dark halos by merging, and the assembly of the baryonic 
mass of galaxies through both gas cooling in halos and galaxy mergers. It includes 
prescriptions for two modes of star formation - quiescent star formation in disks, and 
starbursts triggered by galaxy mergers - and also for feedback from supernovae and 
photo-ionisation. Finally, GALFORM includes chemical evolution of the gas and 
stars, and detailed stellar population synthesis to compute the stellar continuum 
luminosity from each galaxy consistent with its star formation history, IMF and 
metallicity (see Cole et al. 2000 for more details). The unextincted Lya luminosity 
of each model galaxy is then computed from the ionising luminosity of its stellar 
continuum, assuming that all ionising photons are absorbed by neutral gas in the 
galaxy, with case B recombination. 

The semi-analytical approach then allows us to obtain the properties of the Lya 
emission of galaxies and their abundances as a function of redshift, calculating the 
star formation histories for the entire galaxy population, following a hierarchical 
evolution of the galaxy host haloes. In addition, when incorporated into an N-body 
simulation, we also obtain spatial clustering information. This model has been 
incorporated into the largest N-body simulation to date, the Millennium Simulation 
(Springel et al. 2005), to predict clustering properties of Lya galaxies. These results 
will be presented in a forthcoming paper (Orsi et al., in prep.). 

The version of GALFORM which we use here is the one described in Baugh et 
al. (2005) and Le Delliou et al. (2006), with the same values for parameters. The 
parameters in the model were chosen in order to match a range of properties of 
present-day galaxies, as well as the numbers of Lyman Break and sub-mm galaxies 
at 2; ~ 2 — 3. We assume a Kennicutt IMF for quiescent star formation, but a top- 
heavy IMF for starbursts, in order to reproduce the numbers of sub-mm galaxies. 
The only parameter which has been adjusted to match observations of Lya emitters 
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is the Lya escape fraction, which is taken to have a constant value fesc = 0.02, 
regardless of galaxy dust properties. Le Delliou et al. (2006) show that the simple 
choice of a constant escape fraction fesc = 0.02 predicts luminosity functions of 
Lya emitters in remarkably good agreement with observational data at 3 < z < 6. 
Le Delliou et al. (2006) also compared the predicted Lya equivalent widths with 
observational data at 3 < z < 5, including some model galaxies with rest-frame 
equivalent widths of several 100 A, and found broad consistency. For this reason, 
we use the same value fesc = 0.02 for making most of our predictions at z > 7. 
However, since the value of the escape fraction at 2; > 7 is a priori uncertain in the 
models (e.g. it might increase with redshift if high redshift galaxies are less dusty) 
we also present some predictions for other values of fesc- 

Reionisation of the IGM affects predictions for the numbers of Lya emitters in 
deep surveys in two ways: i) feedback from photo-ionisation inhibits galaxy for- 
mation in low-mass halos and ii) reionisation changes the opacity of the IGM to 
Lya photons travelling to us from a distant galaxy, as discussed above. GALFORM 
models the first effect in a simple way, approximating reionisation as being instan- 
taneous at redshift Zreion (see Le Delliou et al. 2006 for more details). We assume 
Zreion = 10, in line with the WMAP 3-year results (Spergel et al. 2007). As was 
shown in Le Delliou et al. (2006; see their Fig. 8), as far as the feedback effect 
is concerned, varying Zreion over the range 7 < Zreion ^ 10 does not have much 
effect on the bright end of the Lya luminosity function most relevant to current and 
planned surveys. For example, varying Zreion between 7 and 10 changes the predicted 
luminosity function at Liya > 10^^'^ erg s~^ by less than 10% for 2; ~ 7 — 10. 



5.3.2 Phenomenological model 

The phenomenological model of Thommes & Meisenheimer (2005; TM05 hereafter) 
assumes that the Lya emitters seen at high redshift are galaxy spheroids seen during 
their formation phase. We summarise the main features here, and refer the reader 
to TM05 for more details. The model is normalised to give the observed mass 
function of spheroids at 2; = 0, which is combined with a phenomenological function 
that gives the distribution of spheroid formation events in mass and redshift. Each 
galaxy is assumed to be visible as a Lya emitter during an initial starburst phase of 
fixed duration (and Gaussian in time), during which the peak SFR is proportional to 
the baryonic mass and inversely proportional to the halo collapse time. The effects 
of the IMF and the escape fraction on the Lya luminosity of a galaxy are combined 
into a single constant factor (i.e. the escape fraction is effectively assumed to be 
constant). With these assumptions, the luminosity function of Lya emitters can be 
computed as a function of redshift. The free parameters in the model were chosen 
by TM05 to match the observed number counts of Lya emitters at 3.5 < 2; < 5.7 
(analogously to the choice of fesc in the GALFORM model). This model does not 
include any effects from reionisation. 
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Table 5.1: Parameters of the fitted Schechter function in previously published pa- 
pers. References are 1) van Breukelen et al. (2005), 2) Gronwall et al. (2007), 3) 
Ouchi et al. (2007), 4) Dawson et al. (2007), 5) Malhotra & Rhoads (2004), 6) 
Shimasaku et al. (2006), and 7) Kashikawa et al. (2006). References 3 — 6 fit for 
three faint end slopes {a = —1.0, —1.5 and —2.0), but here we only reproduce the 
results for fits with a = —1.5 as we fix the slope in our calculations. Malhotra & 
Rhoads (2004) do not give error bars on the fits. Dawson et al. (2007) fix the slope 
to a = —1.6. 
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5.3.3 Observational extrapolation 

Our third approach is to assume that the Lya luminosity function is a Schechter 
function at all redshifts, following 

(j){L)dL = 0*(L/L*)"exp(-L/L*)rfL/L'^ (5.1) 

and to derive the Schechter parameters a, 0* and L* at high redshifts by extrap- 
olating from the observed values at lower redshifts. For our extrapolation, we use 
fits to observations at redshift 2; ~ 3 (van Breukelen et al. 2005; Gronwall et al. 
2007; Ouchi et al. 2007), z = 3.7 (Ouchi et al. 2007), z = 4.5 (Dawson et al. 2007), 
z ^ 5.7 (Malhotra & Rhoads 2004; Shimasaku et al. 2006; Ouchi et al. 2007) and 
z ^ 6.5 (Malhotra & Rhoads 2004; Kashikawa et al. 2006), as found in Table O 
We make linear fits to log(f)* and logL* vs z, and extrapolate to higher redshift. 
For simplicity, we assume a fixed faint end slope of a = —1.5. We do not make any 
corrections for any possible effects of reionisation or IGM opacity. The extrapolated 
values are given in Table 15. 2[ 

5.4 Luminosity functions 

The possible Lya redshifts between z = 7 and 2; = 10 where a narrow-band filter can 
be placed are ZLya = 7.7, 8.2, 8.8, and 9.4. Redshifts beyond 10 are unreachable with 
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Table 5.2: Extrapolated parameters of the observed Schechter function at higher 
redshifts. The faint end slope is fixed to a = —1.5. 



Redshift 


log 0*Mpc"^ 


log L*ergs/s 


7.7 


-3.73 ±0.50 


42.88 ±0.24 


8.2 


-3.80 ±0.50 


42.89 ±0.24 


8.8 


-3.88 ±0.50 


42.91 ±0.24 


9.4 


-3.96 ±0.50 


42.92 ± 0.24 


12.5 


-4.38 ±0.50 


42.99 ± 0.24 



ground-based instruments of the near-future. However, one possibility for z > 10 
surveys may be the James Webb Space Telescope (JWST, see section 15.5.31) and so 
we also make predictions for redshift ZLya = 12.5. 

First, we compare the Lya luminosity functions predicted by the semi-analytical 
(GALFORM) and phenomenological (TM05) models with current observational 
data at 2; ~ 6. This comparison is shown in Fig. 15.21 where we compare the models 
with the cumulative luminosity functions measured in several published surveys at 
z = 5.7 and z = 6.5. We can see that both models match the observational data 
reasonably well, once one takes account of the observational uncertainties. The error 
bars on the observational data points, omitted in the plot in order to not confuse the 
points, are large, at the bright end of the luminosity function due to small number 
statistics, and at the faint end due to incompleteness in the samples. The shallow 
slopes at the faint ends of the Taniguchi et al. (2005) and Ouchi et al. (2007) 
luminosity functions may be due to spectroscopic incompleteness. Both models fit 
the observations well. Hence we conclude that both of these models can be used to 
extrapolate to higher redshifts. 

We now have three methods of extrapolating to higher redshifts, when the direct 
extrapolation of the Schechter function from lower redshifts is included. In Fig. 15.31 
we plot the predicted luminosity functions ai z = 7.7, 8.8 and 12.5 computed by 
these three methods. For other redshifts, the curves may be interpolated. For 
GALFORM, we show predictions for the standard value of the escape fraction f^sc = 
0.02 in the left panel, and for a larger value fesc = 0.2 in the right panel. This 
illustrates the sensitivity of the predictions to the assumed value of fesc at high 
redshift. The predictions from the other two models are plotted identically in both 
panels, since they do not explicitly include the escape fraction as a parameter. 
GALFORM predictions for the numbers of Lya emitters at z > 7 were also given in 
Le Delliou et al (2006). We can see that the predictions from the different methods 
are fairly similar at z = 7.7, but gradually diverge from each other with increasing 
redshift. For the highest redshift, z = 12.5, the TM05 model fails in producing a 
prediction due to numerical problems. We note that making predictions for z = 12.5 
is challenging, for several reasons. Even though only ~ 200 Myrs separate the ages of 
the Universe between redshift 8.8 and 12.5, the Universe went through an important 
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transition at this time as reionisation occurred (Spergel et al. 2007). However, we 
do not know exactly how and when this happened. Also, during this epoch the 
structure in the dark matter (and hence also in galaxies) was building up very 
rapidly. This underlines the interest of obtaining observational constraints at these 
redshifts. 

The hatched regions in Fig. 15.31 show the region of the luminosity function dia- 
gram that has been observationally excluded ai z = 8.8 by Willis & Courbin (2005) 
and Cuby et al. (2006). The former survey was deeper but in a smaller area, whereas 
the latter was more shallow over a larger area, hence the two-step appearance of the 
hatched area. From the plot, it is obvious that their non-detections are perfectly 
consistent with our theoretical models, although the GALFORM model with the 
non-standard escape fraction /e^c = 0.2 is marginally excluded. 

5.5 Future surveys 

In this section, we discuss more specific predictions for several planned and possible 
future surveys. For all calculations, we assume a simple selection on the flux of 
the Lya emission line, with no additional selection on the equivalent width (i.e. we 
include all galaxies with EWiya > 0). We also assume no absorption by the neutral 
hydrogen in the IGM which would reduce the measured fluxes and for GALFORM 
predictions we assume an escape fraction of fesc = 0.02. The predictions from the 
GALFORM and TM05 models for these future surveys as well as some published 
surveys are summarised in Table 15. 3[ 

5.5.1 DaZle — Dark ages z Lyman-Lya Explorer 

DaZle is a visitor mode instrument placed on the VLT UTS (Horton et al. 2004). 
The instrument is designed to use a narrow-band differential imaging technique, 
i.e. observing the same field with two very narrow filters with slightly offset central 
wavelength. Objects with LyLya in one of the filters can then be selected from 
the differential image of both filters. The field-of-view of DaZle is 6.83' x 6.83' and 
it is expected to reach a flux level of 2 x 10^^^ erg s^^ cm^^ (5cr) in 10 hours of 
integration in one filter. This corresponds to a luminosity limit at redshift z = 7.7 of 
log (LLya) = 42.13 erg s~^. The two initial filters are centred on ziya = 7.68 and 7.73 
(with widths Az = 0.006 and 0.025 respectively) and at this redshift, the surveyed 
volume becomes 1340 Mpc^ per pointing per filter pair. Thus, from Fig. 15.31 we 
can conclude that DaZle will discover ~ 0.16 — 0.45 candidates ai z = 7.7 with one 
pointing and filter pair. 

5.5.2 ELVIS - Emission Line galaxies with VISTA Survey 

ELVIS0 is part of Ultra- VISTA, a future ESO Public Survey with VISTA| (Visible 
and Infrared Survey Telescope for Astronomy). Ultra- VISTA is planned to do very 



iwww.astro.ku.dk/~kim/ELVIS.html 
^www. vista. ac.uk 
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deep near-infrared broad- and narrow-band imaging in the COSMOS field. It will 
observe four strips with a total area of 0.9 deg^. The narrow-band filter is focused on 
the ZLya = 8.8 sky background window with central wavelength Ac = 1185 nm, and 
redshift width Az = 0.1. The flux limit of the narrow-band images is expected to 
reach 3.7 x 10"^*^ erg s~^ cm^^ {5a) after the full survey has been completed. Ultra- 
VISTA will run from early 2008 for about 5 years and all the data will be public. 
ELVIS is presented further in Nilsson et al. (2006b). ELVIS will survey several 
different emission-lines (e.g. Ha at redshift z = 0.8, [OIII] at redshift z = 1.4 and 
[Oil] at redshift z = 2.2) as well as the LyLya line. 

When the survey is complete, the final mosaic will reach a LyLya luminosity 
limit of log{LLya) = 42.53 erg s~^. The volume surveyed will be 5.41 x 10^ Mpc^. 
From Fig. 15.31 we see that ELVIS should be expected to detect 3 — 20 LEGOs at 



A possibility even further into the future is to use the James Webb Space Telescopqj 
(JWST). JWST is scheduled for launch in 2013 and will have excellent capabilities 
within the near- and mid-infrared regions of the spectrum. Two of the instruments 
aboard JWST could be used for narrow-band surveys; NIRCam, the near-infrared 
camera, and TFI, the tunable filter imager (for a review on JWST see Gardner 
et al. 2006). NIRCam will have 31 filters, of which nine are narrow-band filters. 
The filter with shortest wavelength has central wavelength Ac = 1.644 /xm (F164N; 
ZLya = 12.5, Az = 0.135). TFI will have tunable filters with variable central 
wavelength, however it is only sensitive at wavelengths larger than A ~ 1.6 fim. 
NIRCam is expected to reach a flux limit of ~ 1 x 10^^^ erg s~^ cm~^ {5a) in 
10000 s of exposure time. Hence, a flux limit of ~ 5 x 10^^^ erg s~^ cm~^ {5a, 
log {LLya{z = 12.5)) = 42.00 erg s~^) could be reached in 10 hours, assuming that the 
sensitivity is proportional to the square root of the exposure time. TFI is expected to 
be able to reach a flux limit almost a factor of two deeper in the same time, however 
it has a field-of-view of only half of the NIRCam (which is 2 x 2.16' x 2.16'). In 
one NIRCam pointing at redshift z = 12.5, approximately 1640 Mpc^ are surveyed. 
Again, from Fig. 15.31 we can estimate that we will detect < 0.1 galaxies per 10-hour 
pointing with NIRCam. However, the number of detections depends strongly on 
the escape fraction which is unknown at such high redshifts, and thus the number 
of detected galaxies can be larger. 

5.6 Constraints on the early Universe 

Of the surveys at these redshifts that have been presented in previous articles (Hor- 
ton et al. 2004; Wilhs & Courbin 2005; Cuby et al. 2006; Nilsson et al. 2006b), 
or are conceivable (JWST, see section 15.5.31) only ELVIS will detect a large enough 
sample to start to measure the luminosity functions and the extent of reionisation 



z = 8.8 



5.5.3 JWST 




www .jwst.nasa.gov 
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at these redshifts and to study the fraction of PopIII stars in the population. We 
here discuss these issues with respect to ELVIS. 

From the semi-analytical modelling, we can make mock observations of the 
ELVIS survey. The procedure to produce these catalogues is explained in detail 
in Orsi et al. (in prep.), but the outline of the process is that galaxies from GAL- 
FORM are placed in matching dark matter haloes in the Millenium N-body sim- 
ulation (Springel et al. 2005), which is a cubical volume in a CDM universe of 
comoving size 500 Mpc/h, thus creating a mock universe with simulated galaxies 
which includes all the effects of clustering. We can then make mock observations of 
this simulated Universe, including the same limits on flux, redshift, sky area etc. as 
for any real survey, and from these observations produce mock galaxy catalogues. 
From the mock catalogues, we can in turn make mock luminosity functions of Lya 
emitters at redshift z = 8.8. In Fig. 15. II we plot the "observed" luminosity functions 
in the 112 mock catalogues taken from different regions of the Millenium simulation 
volume. Note that for making these mock catalogues, GALFORM was run with 
the same cosmological parameters as in the Millenium simulation itself, which are 
slightly different from the "concordance" values assumed elsewhere in this paper, 
as described in the Introduction (this is why the mean luminosity function for the 
whole simulation volume which is plotted in Fig. 15.11 is slightly different from the 
GALFORM prediction for z = 8.8 plotted in Fig. 15. 3p . We used escape fraction 
fesc = 0.02. The figure shows that the spread in number density between the dif- 
ferent mock catalogues is large, almost a factor of ten in number density in each 
luminosity bin. This is a consequence both of the small numbers of galaxies in the 
mock surveys and of galaxy clustering, which causes "cosmic variance" between dif- 
ferent sample volumes. The prediction from GALFORM is therefore that it will be 
difficult to accurately measure the luminosity function of Lya emitters at z = 8.8 
even using the sample from the large area ELVIS survey. In particular, there will 
be no useful constraint on the faint-end slope a. This is simply a consequence of 
the flux limit of narrow-band surveys, i.e. even if we use the median values of the 
luminosity function from all the mocks, then Schechter functions with slopes in the 
range —1 to —2 all give acceptable fits, as illustrated in Fig. 15. 1[ However, if all 
the data are combined in one luminosity bin, it should be possible to measure 0* 
assuming values for a and L*. The possibility of including data points from several 
surveys at different luminosities (e.g. also lensing surveys that probe the faint end 
of the luminosity function) would also significantly improve the results. 

Two suggested methods of constraining reionisation from observations of Lya 
emitters, without requiring spectroscopy, are to measure the clustering of Lya- 
sources and to compare the Lya and UV continuum luminosity functions at these 
redshifts (Kashikawa et al. 2006; Dijkstra, Wyithe & Haiman 2007b; McQuinn et 
al. 2007). McQuinn et al. (2007) show that large HII bubbles may exist during 
reionisation, and that these will enhance the observed clustering of Lya emitters 
in proportion to the fraction of neutral hydrogen in the Universe. A sample of 
~ 50 emitters will be enough to constrain the level of reionisation using this effect 
(McQuinn, priv. communication), almost within reach of the ELVIS survey. A 
future, extended version of ELVIS would be able to place very tight constraints on 
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41.0 41.5 42.0 42.5 43.0 43.5 
log(LLya[erg s"']) 

Figure 5.1: Luminosity functions at z = 8.8 for a set of mock ELVIS surveys com- 
puted using GALFORM. The 112 mock surveys are identical apart from being taken 
from different regions in the MiUennium simulation volume. The open circles show 
number counts in each mock catalogue, in four luminosity bins. The black dot with 
error bars shows the median of the mocks in each bin, with the error bars showing 
the 10-90% range. The thin lines are best fit Schechter functions to the median 
points with different assumed faint end slopes. The thick sohd line shows the "true" 
luminosity function, as measured from model galaxies in the total Millenium simu- 
lation volume. 



reionisation. In Kashikawa et al. (2006) and Dijkstra, Wyithe & Haiman (2007b) the 
use of the combination of the UV and Lya LFs to constrain the IGM transmission is 
explored. Lya emission will be much more susceptible to IGM absorption than the 
continuum emission and thus the ratio between the two LFs will give information 
on the level of IGM ionisation. However, with increasing redshift for Lya, the 
continuum emission will be increasingly difficult to observe, and it is unclear if this 
method will be feasible for surveys such as ELVIS. 

It is possible that galaxies at z = 8.8 still have a significant population of pri- 
mordial PopIII stars. A test for the fraction of primordial stars is the amount of 
Hell 1640 A emission (Schaerer 2003; Tumhnson, SchuU & Venkatesan 2003). De- 
pending on models, these authors predict that the Hell 1640 A emission line should 
have a flux between 1 — 10 % of the flux in the Lya line. For ELVIS z = 8.8 Lya 
emitters, the Hell 1640 A line is redshifted to 1.61 yum. Due to the many OH sky 
emission lines in this region of the spectrum, it would be desirable to try to observe 
the Hell 1640 A hue from a space-based observatory such as JWST. According to 
the JWST homepage, NIRSpec will achieve a sensitivity in the medium resolution 
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Table 5.3: Number of predicted/observed objects per observed field in several 
present and future surveys from two theoretical models. Data from Subaru XMM 
Deep Field (SXDS) are from Ouchi et al. (2005), Sliimasaku et al. (2006) and 
Kasliikawa et al. (2006). GALFORM predictions are made assuming an escape 
fraction of fesc = 0.02. 



Name 


z 


Area 


Luminosity limit 


GALFORM 


TM05 


Observed 






(arcmin^) 


(5(7, erg s~^) 






number 


SXDS-0 


5.7 


8100 


1042.40 


443 


339 


515 


SXDS-S 


5.7 


775 


1042.40 


112 


86 


83 


SXDS-K 


6.5 


918 


1042.27 


108 


57 


58 


DaZle 


7.7 


47 


1042.13 


0.45 


0.16 




ELVIS 


8.8 


3240 


1042.50 


20 


2.8 




Cuby06 


8.8 


31 


1043.10 


0.0003 


0.0 





W&C05 


8.8 


6.3 


1042.25 


0.015 


0.002 





JWST 


12.5 


9.3 


1042.00 


0.018 






mode on 


an emission line at 


1.6 /im of ~ 7 X 


10^19 gj,g g-1 


cm~^ 


(lOcr) for an 



exposure time of 10^ s (30 hours). Thus, if the Lya emission line has a flux of 
~ 5 X 10~^^ erg s~^ cm~^, the Hell 1640 A will be marginally detected with JWST 
in 30 hours of integration, depending on the ratio of Hell 1640 A to Lya flux. The 
NIRSpec sensitivity increases at longer wavelengths, but the increasing luminosity 
distance to galaxy candidates with Hell 1640 A emission at longer wavelengths will 
most likely counteract this effect. 

5.7 Discussion 

We summarise our predictions for number counts of Lya emitters in narrow-band 
surveys in Fig. 15.41 We also summarise the numbers of detected objects for specific 
current and future surveys in Table 15.31 A few comments can be made on the 
differences in predictions between the two models. Firstly, as can be seen in Fig. 15.41 
and also Fig. 15.31 the luminosity functions have steeper faint-end slopes in the 
GALFORM models than in the TM05 models. Secondly, the GALFORM and TM05 
models predict similar amounts of evolution at a given flux over the range z = Q — 9 
where they can be compared. 

Several factors enter into the error bars of our predictions. One problem is the 
uncertainties in, and disagreement between, the observed lower redshift luminosity 
functions which are used to calibrate the theoretical models. There are many caveats 
in producing Lya luminosity functions, of which the selection function is the most 
difficult to correct for. The problem arises from that the filter transmission curve is 
not box-shaped, but rather gaussian. Thus, only brighter objects will be observed 
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log(LLya [erg s"']) 



Figure 5.2: Plot of luminosity functions at redshifts z = 5.7 and 6.5. Red points 
and lines are at redshift z = 5.7, black points/lines at redshift z = 6.5. Points are 
observations by Ajiki et al. (2003; redshift 5.7, squares), Hu et al. (2004; redshift 
5.7, pluses), Taniguchi et al. (2005; redshift 6.5, crosses), Shimasaku et al. (2006; 
redshift 5.7, diamonds), Kashikawa et al. (2006; redshift 6.5, stars), Malhotra & 
Rhoads (2004; redshift 5.7 and 6.5, triangles) and Ouchi et al. (2007; redshift 
5.7, filled squares). Solid lines show the GALFORM model (with escape fraction 
fesc = 0.02), dot-dashed lines the TM05 model. Note that the Taniguchi et al. 
(2005) and the Ouchi et al. (2007) samples are the spectroscopic samples only. 



at the wings of the filter, and these will be observed to have smaller than intrinsic 
luminosities. Secondly, the equivalent width (EW) limit that the survey is complete 
to depends on the depth of the broad-band images used for the selection. Thirdly, 
if the sample is a photometric sample, it is possible that there are lower redshift 
interlopers, where the emission line is e.g. [Oil], in the sample. Finally, the samples 
are still so small that we have to deal with small number statistics. All of these 
problems cause the observed luminosity function at lower redshifts to be uncertain. 

Both theoretical models (semi-analytical and phenomenological) have uncertain- 
ties resulting from how they model the galaxy formation process, and also from the 
assumption that the fraction of Lya photons escaping from galaxies is constant and 
does not change with redshift. In addition, neither model includes attenuation of 
the Lya flux due to neutral hydrogen in the IGM. This attenuation would be ex- 
pected to be strong at z > Zreion, when the IGM is neutral, and weaker at 2; < Zreion, 
when most of the IGM is ionised. The degree of attenuation depends on a number 
of different effects, as analysed in Santos (2004), and discussed in Le Delliou et al. 
(2006), and is currently very uncertain. Nonetheless, this attenuation is expected 
to produce observable effects on the evolution of the Lya luminosity function, if 
reionisation occurs within the redshift range covered by future observations, and so 
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Figure 5.3: Predicted Lya luminosity functions at z > 7. Red lines are extrapo- 
lations from observed luminosity functions at lower redshift, green lines are TM05 
models and black lines are GALFORM models. Different linestyles show different 
redsliifts z = 7.7, 8.8 and 12.5. No prediction is shown for the TM05 model at 
z = 12.5. Hatched area marks observational upper limits from Willis & Courbin 
(2005) and Cuby et al. (2006), both at redshift z = 8.8. In the left panel, the GAL- 
FORM predictions are shown for escape fraction f^sc = 0.02 (our standard value), 
while in the right panel, they are shown for /g^c = 0.2. The predictions from the 
other two methods are identical in both panels. 




6 8 10 12 14 



Figure 5.4: Summary of predictions. The plot shows the number of Lya emitting 
galaxies expected per square degree per redshift interval Az = 0.1 as a function of 
redshift and observed flux limit. The predictions of GALFORM are shown in black 
and of the TM05 model in red. The different line styles are for different flux limits. 
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estimating the reionisation redshift and the neutral fraction after reionisation are 
included in the science goals of these surveys. 

It is apparent that the key to acquiring a large sample of Ly«-emitting galaxies at 
redshifts greater than 7 is both depth and area. In a recent paper, Stark, Loeb & Ellis 
(2007) suggest that one of the most efficient means of finding very high redshift Lya 
emitters is through spectroscopic surveys focused on gravitational lensing clusters. 
Lensing surveys could easily reach down to a luminosity limit of 10^°'^ erg s~^ in a 
few tens of hours. However, the surveyed volumes are very small, of the order of a 
hundred Mpc'^. For a lensed survey, the area in the source plane is reduced by the 
same factor that the flux is amplified, so in principle one gains in the total number of 
objects detected relative to an unlensed survey if the luminosity function is steeper 
than A^(> L) oc L'^. In the GALFORM and TM05 models, the asymptotic faint- 
end slope is shallower than this, but at higher luminosities, the slope can be steeper. 
For example, GALFORM predicts that at z = 10, the average slope in the luminosity 
range 10^^-10^^ erg s~^ is close to N{> L) oc (see Fig. 8 in Le Delliou et al. 
2006), so that a lensing amphfication of 10 results in 10 times more objects being 
detected, with intrinsic luminosities 10 times lower, compared to an unlensed survey 
with the same area and flux limit. Therefore lensing and narrow-band surveys are 
complementary to each other as they probe different parts of the luminosity function. 
With either type of survey, reaching a signiflcant sample of redshift z ~ 7 — 8 should 
be possible in the next few years with telescopes/instruments in use or soon available. 

An interesting type of object found recently in narrow-band surveys are the 
Lya blobs, large nebulae with diameters up to 150 kpc and Lya luminosities up to 
]^q44 gj,g g-i -^jI;]^ QTj. without counterpart galaxies (e.g. Steidel et al. 2000; Matsuda 
et al. 2004; Nilsson et al. 2006a). Several mechanisms have been proposed to explain 
this phenomenon, including starburst galaxies and superwinds, AGN activity or cold 
accretion. It is interesting to consider if such objects would be detected in any of 
these surveys, assuming they exist at these redshifts. A typical Lya blob will have 
a luminosity of ~ 10^^ erg s^^ and a radius of, say, 25 kpc. This will result in a 
surface brightness of ~ 5 x 10^^ erg s~^ kpc^^. Thus, a narrow-band survey will have 
to reach a flux limit, as measured in a 2" radius aperture of ~ 1.3 x 10"^^ erg s~^ at 
redshift z — 8.8, corresponding to logL = 42.11. (An aperture radius of around 2" 
is expected to be roughly optimal for signal-to-noisc.) For lower or higher redshifts, 
this limit is higher or lower respectively. Thus, ELVIS will not be able to detect 
Lya blobs unless they arc brighter and/or more compact at higher redshift than a 
typical blob at lower redshift. DaZle and JWST could in principle detect this type 
of object, but only if they are very abundant in the very high redshift Universe, due 
to the small survey volumes of these instruments. It is of course highly uncertain 
what properties such Lya blobs would have at 2; ~ 7 — 9, or their space density, but 
it appears unlikely that the future surveys presented here would detect any such 
objects. 

To flnd compact Lya emitters at redshifts 2; > 10 in signiflcant numbers we will 
probably have to await instruments even further in the future. If a future 40-m 
ELT (Extremely Large Telescope) was equipped with a wide-fleld NIR imager and 
a narrow-band fllter of similar width to ELVIS, it could reach a luminosity limit of 
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L ~ 10^^'^ erg at redshift z = 10.1 (where a suitably large atmospheric window 
exists) in approximately 20 hours. Using the GALFORM model for z = 10, the 
number density should be N(>L)^ 4 x 10~^ Mpc~^ at this luminosity limit. Thus, 
to get a sample of ten Lya emitters would require imaging an area on the sky 
of approximately 16 square arcminutes, assuming a narrow-band filter with redshift 
range 10.05 < ZLya < 10.15. This could be achieved with one pointing if the detector 
has a field-of-view of 6 arcmin on a side, as suggested by the ESO ELT Working 
Group0. It should of course be noted that these are very tentative numbers, but 
they display the possibilities of far future instruments. 



'*http://www.eso.org/projects/e-elt/Publications/ELTJNSWG_FINAL_REPORT.pdf 
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Selection methods for Lya 
emitters 

This paper is in preparation. The authors will be Nilsson, K.K., MoUer, O., Hayes, 
M., M0ller, P., & Fynbo, J.P.U. 

6.1 Abstract 

Conteoct Narrow-band surveys for Lya emitters is an increasingly popular method 
to find high rcdshift galaxies. However, different groups have presented different 
methods to find Lya galaxy candidates. We here present a method to determine 
the best selection criteria, combining narrow- and broad-band observations. 
Aims To find the optimal selection sub-space in colour-colour space for Lya emitters, 
i.e. to determine the best colour selection criteria for this type of object. 
Methods We simulate a galaxy population by constructing a large number of galaxy 
spectra and convolve with filter profiles to get the "observed" fiuxes. Some spectra 
have Lya emission superposed, which is then observed through a narrow-band filter 
profile. We study the distribution in A'^ dimensional colour space and attempt to 
find the angle in which the scatter in the normal galaxy population is the least and 
the emission-line candidates are most easily and most accurately selected. 
Results Wc present the three most favourable combinations of galaxy colours to be 
used to select i) Lya emitters from field galaxies and ii) Lya emitters from [Oil] 
emitters. We also present the optimal selection criteria to use with these colours. 
Conclusions All of the best combinations of colours for selecting Lya emitters from 
field galaxies include sampling the UV slope of the galaxy with two broad-band 
filters, one on either side of the Lya line. It is concluded that selection with 
two narrow-band/broad-band colours is superior to selection with only one narrow- 
band/broad-band colour. 

6.2 Introduction 

Narrow-band imaging surveys have become increasingly successful in discovering 
high and very high redshift galaxies. There are now in total several hundreds of 

Q1 
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Table 6.1: Parameter space for field galaxy population. 



Min. value Max. value Number of steps 



Redshift (z) 
Dust (E(B-V)) 
Age (Gyr) 



3 50 

1 20 

0.001 10 50 (log) 



spectroscopically confirmed Lya emitters at redshifts z ~ 3 (e.g. Steidel et al. 
2000; Fynbo et al. 2003; Matsuda et al. 2005; Venemans et al. 2007; Nilsson 
et al. 2007), z ~ 4.5 (Finkelstein et al. 2007), z ~ 5.7 (Malhotra et al. 2005; 
Shimasaku et al. 2006; Tapken et al. 2006) and z ~ 6.5 (Taniguchi et al. 2005; 
Kashikawa et al. 2006). Most candidates are selected using one narrow-band filter 
and one broad-band filter, where objects that are bright in the narrow-band filter 
but comparatively faint or non-detected in the broad-band filter are selected. In 
some cases, the selection is based on determining the equivalent width (EW) of 
the potential emission line. However, the calculation of the EW of an emission- 
line observed through a narrow-band /broad-band filter set-up has been shown to be 
uncertain (Hayes & Ostlin 2006). In 1993, M0ller & Warren (1993) suggested for the 
first time that two colours should be used for detecting emission-line galaxies, i.e. 
that any object should be observed with a narrow-band filter and two broad-band 
filters. This method was subsequently explored in Fynbo et al. (2003) and Nilsson 
et al. (2007). However, present day and future multi- wavelength surveys such as e.g. 
the GOODS (Giavalisco et al. 2004) and COSMOS (Scoville et al. 2006) surveys 
warrant a review on the selection criteria for narrow-band surveys in fields where 
public, multi-wavelength data exists. 

This Letter is organised as follows; in section 16.31 we present the models we use 
to create colour plots for "normal" galaxies and Lya emitters as well as the method 
by which we find the most efficient selection colours. In section 16.41 we present our 
results and in sections 16.51 we discuss the results and draw conclusions. 

Throughout this paper, we assume a cosmology with Hq = 72 km s~^ Mpc^^ 
(Freedman et al. 2001), flm = 0.3 and = 0.7. Magnitudes are given in the AB 
system. 



6.3.1 Creating a mock sample 

To study which selection criteria are most successful in determining Lya emitter 
candidates, we create model spectra using the Starburst99 models (Leitherer et al. 
1999; Vazquez & Leitherer 2005). For the "normal" galaxy population, we create 
50'000 spectra, evenly covering the parameter space given in Table 16. 1[ The age 
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Table 6.2: Redshifts used for Lya emitters. The redshifts are selected to represent 
a wide range of redshifts that have or will be surveyed. Representative references 
are given. 



Redshift (zLya) 


Reference (e.g.) 


2.4 


Francis et al. 2004 


3.1 


Nilsson et al. 2007 


4.5 


LALA; Finkelstein et al. 2007 


5.7 


Shimasaku et al. 2006 


6.5 


Kashikawa et al. 2006 


7.7 


DaZle; Horton et al. 2004 


8.8 


ELVIS; Nilsson et al. 2006b 



of each model galaxy is constrained to be less than the age of the Universe at the 
redshift it is at. The mass of a galaxy is a simple multiplication factor and as we here 
only work with colours, i.e. flux ratios, the mass is irrelevant except in the cases 
when the mass is so small that the flux in a certain waveband is lower than the 
flux detection limit. We discuss this further in sec. 16. 5[ All spectra are calculated 
using single stellar populations (SSP), i.e. all the stars are created instantaneously. 
This will cause the spectrum to be blue to the extreme at a young age and red to 
the extreme when old, thus widening the distribution of the fleld galaxy population 
with respect to other more extended SFHs. Particularly at late times colours will 
be redder, with fewer hot stars present in the population. 

For the Lya emitters, we create galaxy spectra with the same range of ages and 
dust amounts as for the fleld galaxies and redshift them to the redshift in question. 
We then add Lya by calculating the expected Ha EW from the age of the starburst, 
convert this to a Lya EW using case B recombination (Brocklehurst 1971) and add 
such a line to the model spectra. The Lya EWs range between — 240 A (Chariot 
& Fall 1993). As we are also interested in distinguishing between Lya emitters and 
interloper galaxies such as [OII]-emitters, we also add an [OII]-emission line to a 
sample of galaxy spectra with EWs ranging between — 100 A, as [Oil] is the most 
common interloper in Lya emission surveys. We create 50'000 Lya emitters and 
[Oil] emitters respectively. We wish to decide the optimal selection criteria for a 
variety of Lya redshifts, see Table 16. 2[ that have been or will be surveyed. 

For each galaxy spectrum, we convolve the spectrum with a set of standard 
fllter transmission curves in order to calculate magnitudes and colours. The fllters 
used correspond, for broad-bands, to the fllter set used in the GOODS Surveys, see 
Table 16.31 and Fig. 16. H and are very similar to standard Johnson- Cousins fllters. 
The narrow-band fllters are deflned to be gaussian shaped, with central wavelength 
centred on the particular Lya redshift desired and with full-width half-maximum 
(FWHM) of 1% of the central wavelength. 
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Table 6.3: Filters used in this paper. These filters are used in the convolution with 
the synthetic spectra in order to calculate fluxes and colours in typical observed 
bands. 



Filter name 


Central Wavelength 


FWHM 


U {ESO 2.2-m) 


3630 A 


760 A 


B (F435W, HST) 


4297 A 


1038 A 


V (F606W, HST) 


5907 A 


2342 A 


z (F814W, HST) 


7764 A 


1528 A 


z' (F850LP, HST) 


9445 A 


1230 A 


J ( VLT) 


1.25 /im 


0.6 //m 


H ( VLT) 


1.65 //m 


0.6 //m 


Ks (VLT) 


2.16 jim 


0.6 jim 


Chi (Spitzer) 


3.58 /xm 


0.75 /im 


Ch2 (Spitzer) 


4.50 fim. 


1.02 iim 


ChS (Spitzer) 


5.80 //m 


I A?, 


ChA {Spitzer) 


8.00 nm 


2.91 //m 




Wavelength (/xm) 



Figure 6.1: Transmission of selection filters from the GOODS-S data-set. 



6.3.2 Optimal distinction of Lyo; emitters 

In past work (e.g. Fynbo et al. 2003; Nilsson et al. 2007) a distinction between 
LycK emitting candidates and non-emitters was performed by making a colour cut 
in two narrow- minus broad-band colours. The reasoning was that, due to the 
emission in the narrow-band (NB) filter, the difference between Lya emitters and 
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non-emitters should be maximal. However, this was based on a purely qualitative 
argument, rather than a careful analysis. In this Letter we want to perform a more 
quantitative study using the evolutionary models developed. The data at hand 
are three populations of galaxies with, in total, 78 colour combinations created 
from mock magnitudes in the bands presented in Table 16.31 combined with a NB 
filter. Thus, we wish to discriminate between three populations of galaxies in a 78 
dimensional space. The goal is to find two colours in which the populations can 
be plotted where the largest separation between the three populations occur. We 
choose two colours as this allows a significant number of combinations - more than 
one colour - but does not require a large amount of additional data for the selection. 

For each combination of colour, the catalogues of galaxies is smoothed using a 
gaussian filter onto a 100 x 100 grid with a step size of 0.1 in magnitude difference. 
We then optimise the selection of Lya emitters by calculating a "merit" number M 
which is the ratio of distance between the mean colour of each population over the 
area of the overlap between the populations, see Eq. 16. 1[ 



In this equation, C are the colours and A is the area covered by the population 
in colour space. The best selection colours will have the highest merit number 
M. To calculate the optimum, we need to know the distribution of samples of both 
populations. This can in principle be obtained observationally. However, in the case 
here, for our evolutionary models, we do not have that information. The sample 
of normal galaxies will for example outline a region in the parameter space where 
outliers are extreme types of galaxies. Within these boundaries, we do not know 
the actual number distribution of galaxies in the colour space. Rather we assume 
that within the region covered by the boundaries, the distribution is uniform. 



For each redshift, all combinations in two dimensions that include the narrow-band 
magnitude in both combinations (in order to minimise the number of observed bands 
needed) are ranked according to their merit number, M, as described in the previous 
section. We also run the code with only [OII]-emitters and Lya emitters with the aim 
to find an optimal combination of colours for separating the two types of emitters 
photometrically. For the emitters, we make a number of cuts in the properties of 
the spectra. We set the maximum age of a Lya emitter to be 40 Myrs and the 
maximum age of an [OH] emitter to be 100 Myrs. The maximum E(B-V) allowed 
is 0.3 for both Lya and [OH] emitters. Also, for both populations of emitters, the 
minimum EW of the emission line is taken to be 10 A. 

An illustration of the best combinations for redshift z = 2A can be seen in 
Fig. 16. 2[ The plot shows the top combinations of colours to select between Lya 
emitters and field galaxies or [OH] emitters. As can be seen, the overlap is very small 
with the field galaxy population. The overlap between the two emitter populations 
is larger. The best three selection combinations for our selected redshifts can be 



M 



< Ciya > — < CpG > 
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Figure 6.2: Illustration of selection method for redshift z = 2.4. The left panel shows 
the optimal selection colours to distinguish between Lja emitters and field galaxies. 
The right panels shows the best combination to distinguish between Lja emitters 
and [Oil] emitters. Black points are field galaxies, red points are Lya emitters and 
green points are [OII]-emitters. 



found in Table 16.41 In this table, the merit number has been normalised to the most 
favourable value for each redshift. The merit number for the distinction between Lja 
emitters and [OII]-emitters has been calculated the same way as for the distinction 
with field galaxies. For each best solution, we also want to find the actual selection 
criteria to apply. Again, as the actual distribution of galaxies within colour space is 
unknown, it is impossible to decide different confidence levels at which the criteria 
can be applied. Instead, the only possibility is to calculate the 100% confidence 
limit, i.e. where there is no contamination of the two populations. For the best 
solutions in Table 16.41 these criteria are given in Table 16. 5[ 



6.5 Discussion 

The motivation of this Letter is to study how to best select Lya-emitter candidates 
photometrically. The necessity of such a review is two-fold; firstly, several surveys 
have or are being conducted that offer a large set of public, multi-wavelength data in 
fields or varying size on the sky. These data-sets can be used to make very accurate 
photometric redshifts and should hence provide excellent opportunities to select 
Lya emitting candidates from narrow-band surveys with great accuracy. Secondly, 
narrow-band surveys will try to reach fainter fluxes and larger area/samples at all 
redshifts, and also to much higher redshifts {ziya ^7). In both cases spectroscopic 
follow-up to confirm the candidates can be challenging or even beyond the capacity 
of spectrographs. Without spectroscopic confirmation, it is necessary to make very 
accurate photometric selections with as few interlopers as possible. 

Hence, we have here reviewed the most efficient colour combinations that can be 
made with a large multi-wavelength data-set for selecting Lya emitters from field 
galaxies, as well as distinguishing between Lya emitters and [OII]-emitters. For the 
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Table 6.4: Best selection methods for Lya emitters at different redshifts. The 
columns give the first, second and third best selection colours. NB denotes the 
narrow-band filter at the specified redshift. The numbers in the parentheses are the 
relative goodness of the selection technique normalised to the best solution. The 
upper part of the table give the solutions for selecting between Lya emitters and 
field galaxies and the lower part between Lya emitters and [OII]-emitters. 



z Best solution Second solution Third solution 

Lya emitters vs. Field galaxies 

2.4 NB-U vs. NB-B (1.00) NB-B vs. NB-V (0.85) NB-B vs. NB-i (0.67) 
3.1 NB-B vs. NB-V (1.00) NB-V vs. NB-i (0.86) NB-B vs. NB-i (0.83) 

4.5 NB-B vs. NB-J (1.00) NB-B vs. NB-K^ (1.00) NB-B vs. NB-H (1.00) 
5.7 NB-U vs. NB-H (1.00) NB-B vs. NB-H (0.76) NB-B vs. NB-Chl (0.76) 
6.5 NB-U vs. NB-H (1.00) NB-V vs. NB-H (0.83) NB-V vs. NB-Ch2 (0.83) 

7.7 NB-U vs. NB-K^ (1.00) NB-B vs. NB-H (0.93) NB-B vs. NB-K, (0.93) 

8.8 NB-B vs. NB-K, (1.00) NB-V vs. NB-K, (0.55) NB-i vs. NB-C/i3 (0.37) 

Lya emitters vs. [011]-cmitters 

2.4 NB-U vs. i-z' (1.00) NB-U vs. NB-B (0.98) NB-B vs. i-z'i (0.93) 
3.1 NB-C/il vs. U-B (1.00) NB-K, vs. U-B (0.57) NB-C/i2 vs. U-B (0.57) 

4.5 NB-K, vs. U-V (1.00) NB-J vs. B-V (0.51) NB-K, vs. B-V (0.51) 
5.7 NB-U vs. Y-Ch3 (1.00) NB-U vs. Y-Ch2 (1.00) NB-U vs. Y-Chl (1.00) 
6.5 NB-U vs. V-Chl (1.00) NB-U vs. V-Ch2 (1.00) NB-U vs. i-ChA (1.00) 

7.7 NB-U vs. V-Chl (1.00) NB-U vs. V-Ch2 (1.00) NB-U vs. B-Chl (0.94) 

8.8 NB-U vs. U-C/i2 (1.00) NB-B vs. \J-Ch3 (0.96) NB-U vs. B-Ch2 (0.96) 



first part of the results, it is interesting to note that all of the best selection criteria 
for Lya emitters/field galaxies involve one combination with a broad-band filter on 
the blue side of the narrow-band filter and one combination with a broad-band filter 
on the red side of the narrow-band filter. For the best combination for all redshifts, 
the filter on the red side of the Lya line is also always located on the blue side of 
the Balmer break in the spectrum. This result can be understood intuitively, as two 
broad-band magnitudes on either side of the Lya line, but both blue-ward of the 
Balmer break, will constrain the UV slope of the galaxy and thus one easily gets 
a measure of the equivalent width of the emission line when this value is combined 
with the narrow-band flux. 

For redshift z — 3.1 the best selection colour is identical to what has been 
used in other publications (Fynbo et al. 2003; Nilsson et al. 2007). Fynbo et al. 
(1999; 2002) have also used a two-colour approach in trying to find Lya emitters at 
redshifts z ^ 2. They applied the colours NB-U and NB-i with successful results. 
This criteria is similar to the third best solution that we find for redshift z — 2A. At 
other redshifts, the method of selecting emitters using two filters has not previously 
been used. However, it is obvious from our results that a selection with two colours 
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Table 6.5: Selection criteria. The selection criteria here are the "100% confidence 
levels" so that the populations are separate with the combinations of colours corre- 
sponding to the best selection colours as presented in Table 16. 4[ 



Kedsnirt 


beiection 


criteria (Lya/FG) 


2A 


iNr>-r)< - 


-O.z n (iNB-Bj— O.ox (iNB-U j< 0.25 


6.1 


i\lr>-V< - 


-0.7 U (iNB- Vj +0.45 X (IN B-B J < —0.7 


4.5 


i\B-J< — 


1 O 1 1 /' ATT") T\ 1 (\ AO-,/ A.TT) T~) \ ^ 1 O 

1.2 U (i\lB-J j +0.48 X (IN B-B J < —1.2 


0. 1 


iNJD-rl< - 


-i.Z U (^iNiD-rl J + U.oUX (^iNrS-U J< — L.Z 


6.5 


NB-H< - 


-1.2 U (NB-H)+0.30x(NB-U)< -1.2 


7.7 


NB-K, < 


-1.2 U (NB-K,)+0.22x(NB-U)< -1.2 


8.8 


NB-K, < 


-1.0 U (NB-k/)+0.23x(NB-B)< -1.0 




Selection 


criteria (LyQ;/[OII]) 


2.4 


NB-U< - 


-1.6 


3.1 


U-B> 0.8 


U (U-B)-0.13x(NB-C/il)> 0.8 


4.5 


U-V> 1.2 


: U (U-V)-0.33x(NB-K,)> 1.0 


5.7 


NB-U< - 


-2.5 U (V-C/i3)+2.71x(NB-U)< -5.9 


6.5 


NB-U< - 


-3.0 U (V-C/il)+1.60x(NB-U)< -3.2 


7.7 


NB-U< - 


-4.0 U (U-C/il)+2.15x(NB-U)< -4.0 


8.8 


NB-U< - 


-4.0 U (U-C/i2)+2.15x(NB-U)< -4.6 



including narrow-band/broad-band combinations is better than a selection with only 
one narrow-band/broad-band colour as is most commonly used. This is obvious since 
our method finds selection criteria for all redshifts (in Table 16. 5p where both colours 
are included. If one colour would be the best selection combination, the criteria 
would have been a cut parallel to either the x— or y—axis. This is not the case. 
This is also an intuitive result. More information will always yield a more accurate 
result. In future work it would be interesting to further quantify how much better 
results a two colour selection yields compared to a one colour selection. 

As for the distinction between Lya emitters and [Oil] emitters, for several red- 
shifts the best selection requires the narrow-band minus U band colour. This is 
presumably because Lya emitters with redshifts above z ~ 2 should at all times 
have extremely faint emission in U, as this filter will sample the far-UV, whereas 
[OII]-emitters are at lower redshift and will thus have continuum emission in the 
observed U filter. 

The method used here of maximising the merit number M has the advantage 
that it is very simple to calculate and to understand. It also involves very few bands, 
in our case only three for the selection between Lya emitters and field galaxies, and 
so it can be applied to almost any survey. The disadvantage is that for large multi- 
wavelength surveys such as GOODS or COSMOS, it does not make use of all the 
data available. Better solutions may be found if more bands are allowed. The 
method is not very rigorous mathematically, nor is it a simple task to define the 
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merit number. A more careful analysis would take a more complex approach, using 
a principal component analysis such as e.g. Fisher Linear Discriminants (Fisher 
1936). 

A caveat with the input to the method is the unknown mass distribution. In 
order to facilitate the calculations made here, all galaxies were assumed to have a 
mass of = IO^^Mq. The catalogues were also cut at a magnitude of magAB = 30 
which introduces another arbitrary cut in mass properties. Thus, we do not know 
the actual distribution of galaxies in colour space, and we will miss extreme cases 
of colours in low or high mass galaxies. The next step would thus also be to try 
to incorporate a mass function of galaxies into the model. Yet another caveat is 
the lack of understanding of the properties of Lya emitters and [Oil] emitters. The 
true age and dust content distribution of these galaxies are not known. Previous 
SED fitting of Lya emitters have shown them to be in general young, low mass and 
dust free (Gawiser et al. 2006; Lai et al. 2007; Finkelstein et al. 2007; Nilsson et 
al. 2007). However, this may be a selection bias arising from the selection criteria 
which often include high EW and blue colours. Finally, in future work it would be 
desirable to try to incorporate the effects of IGM opacity in the spectra. This effect 
will be most important at very high redshifts, but also very uncertain as we have 
no conclusive information about how and when re-ionisation happened. 
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Chapter 7 
Conclusions 



The technique of finding high redshift galaxies through their Lya emission and 
narrow-band imaging, and hence the class of high redshift galaxies called Lya emit- 
ters, is a relatively young branch of observational cosmology. Even though it was 
suggested already by Partridge & Peebles in 1967, the first successful observations 
did not occur until the beginning of the 1990's. Since the year 2000, many Lya emit- 
ters have been published though, reaching in total more than 500 spectroscopically 
confirmed galaxies between redshift 2 < z < 7, including the highest redshift galaxy 
recorded to date at z = 6.96 (lye et al. 2006). The relative ease with which this type 
of galaxy can be found, i.e. with moderate exposure times on medium sized to large 
telescopes, promises to make this technique one of the best selection techniques for 
future high and very high redshift galaxy surveys. This in turn will enable a better 
understanding of the young Universe and galaxy formation and evolution through 
the study of these galaxies, and an insight into the so-called "Dark Ages" , before 
the Universe was re-ionised. 

However, although the sample of Lya emitters is steadily growing, relatively little 
is known about the nature of these galaxies, or their relation to other classes of high 
redshift galaxies. In the first part of this thesis, we attempt to study a sample of 
redshift z = 3.15 Lya emitters found through narrow-band imaging in GOODS-S, 
a field with a large amount of public, multi- wavelength data (Chapter [2] and [3]). 
The main objective to get this narrow-band image was to use the multi-wavelength 
data to try to understand the nature (e.g. the masses, ages, star formation rates, 
redshift distribution and spatial distribution) of these galaxies, and to compare their 
properties to those of other classes of high redshift galaxies. In part, this goal has 
been achieved. The following main results have been presented in this thesis. 

• SED fitting analysis: properties of Lya emitters at 2^ = 3.15 

In Chapter [3], the multi-wavelength data was used to perform an SED fitting of 
the spectra of the Lya emitters, thus allowing a study of their properties, including 
stellar mass, metallicity and dust content. We found that the stacked SED of these 
Lya-emitters, which covered the spectra of the galaxies from the restframe far-UV to 
restframe near-infrared, was best fitted by a very low metallicity {Z = 0.005 x Zq), 
had low dust extinction {Ay = 0.26'tQ \\), medium ages (0.851q;42 Gyrs) and medium 
stellar masses (M.,, = x 10^ M©). These results are in reasonably good 



101 



102 



Chapter 7 




5.0x10° 1.0x10^ 1.5x10 



2.0x10^ 2.5x10^ 



Stellar mass (Mg) 



Figure 7.1: Mass function of GOODS-S LEGOs. 



agreement compared to results from other studies (Gawiser et al. 2006; Lai et al. 
2007; Finkelstein et al. 2007), see also Table fT73l In one respect our results differ 
from the others; our age is higher. However, the age of an SED is the parameter 
which is the hardest to constrain, as can also be seen in Fig. 13. 7[ Thus, relatively 
little weight should be put on these differences, as well as the actual ages themselves. 

To expand on the work presented in the paper in Chapter [3l we here evaluate the 
mass function of our Lya emitters. The properties cited above are for the stacked 
SED, hence a sort of average of all the properties of all candidates. It is interesting 
to consider what the distribution of masses for these objects look like. To do this, we 
take the mass derived from the SED fitting, multiply it with the number of objects 
in the stack (this number is 23) and calculate individual masses by assuming that 
the stellar mass of the galaxy is proportional to the V band luminosity. The reason 
for choosing the V band as a mass estimator is because it is in this band that we 
have the largest signal-to-noise detections of the candidates. This is not the best 
estimator, as the observed V band corresponds to the restframe UV and the galaxy 
light in this regime can be heavily absorbed by dust, i.e. a more massive galaxy 
with a large amount of dust would appear smaller in this calculation. However, it is 
the only band in which we can do this analysis, as it is the only band with adequate 
signal-to-noise and most detections. The resulting distribution can be seen in the 
histogram plot in Fig. 17.11 It appears that a large fraction, more than 20%, of the 
total mass of the sample is concentrated to one object (LEG0_G00DS-S#14), and 
more than half the sample has masses below M,, < 3 x 10* Mq. With these results 
in mind, it will be of great interest in the future to try to improve the statistics on 
this kind of Lya emitter mass function, to see if this type of massive Lya emitters 
are very rare occurrences or if they are more abundant. Whatever the result, it can 
be said with certainty that the majority of Lya emitters have small masses. 

We can compare the results of our SED fitting with similar results from LEG 
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surveys. Three such surveys have been pubhshed in Papovich, Dickinson & Ferguson 
(2001), Shapley et al. (2001) and Verma et al. (2007). In all papers, the authors 
use samples of spectroscopically confirmed LBGs at redshift 2; ~ 3 or 5 respectively. 
They use Bruzual & Chariot (2003) models, with a Calzetti extinction law and 
metallicity fixed at solar levels. They fit for star formation history, dust, age and 
stellar mass. The median results for the sample of 74 LBGs presented in Shapley et 
al. (2001) are E{B -V) = 0.155, age = 320 Myrs and ~ 2 x 10^° Mq. Papovich, 
Dickinson & Ferguson (2001) find similar results, but with a larger dust content of up 
to Av ~ 2. Verma et al. (2007) find their LBGs to be less massive (M^ ^ 2 x 10^ Mq) 
and younger, but with larger star formation rates. To compare the dust results, one 
has to take into account that the ratio between Ay and E{B — V) is R ~ 3.1. 
Thus, the dust extinction in the LBGs of Shapley et al. (2001) is approximately 
Ay 0.48, in good agreement with the results of Verma et al. (2007). This 
is similar, or a bit higher than the value found in Chapter [31 The ages of the 
LBGs are slightly lower than those of the Lya emitters, although this is a difficult 
quantity to constrain, as discussed above. Finally, a large difference also lies in the 
stellar masses. LBGs at 2; ~ 3 are much more massive than Lya emitters at the 
same redshift. Another difference is the star formation rate. Papovich, Dickinson 
& Ferguson (2001), Shapley et al. (2001) and Verma et al. (2007) all find star 
formation rates in the range of a few ten to several hundred solar masses per year, 
whereas we find star formation rates less than ten solar masses per year. Thus, there 
is a clear difference in the properties of LBGs and Lya emitters at this redshift; LBGs 
are more massive, more star forming and possibly more dusty than Lya emitters 
although the two types of galaxies have similar ages. 

• LEGO_GOODS-S#16: A Lya emitting dusty starburst galaxy? 

During the analysis of the multi- wavelength data of the candidates in the GOODS-S 
field, we discovered one candidate with a very unusual SED profile, LEGO_GOODS- 
S#16. This candidate was undetected in the bluest band, the HST B band, and 
marginally detected in the V band with the greatest signal-to-noise. But as we 
followed the candidate in the redder bands, its flux increased dramatically, see also 
Fig. 13.91 Although with a small offset between Lya and counterpart centroid, it 
appears as if this is a very dusty starburst galaxy with Lya emission emerging 
from one side. We are hoping to follow this detection up with future mm/sub- 
mm observations, that could hopefully give more insight into the nature of this 
anomalous object. 

• Discovery of filamentary structure: Two parallel filaments at z = 3.15 

A further result from that paper was the discovery of a filamentary structure, con- 
sisting of two, seemingly parallel filaments stretching beyond the length of the image. 
Filaments have been observed before (M0ller & Fynbo 2003; Matsuda et al. 2005), 
but never two parallel filaments. It would be of great interest to get spectroscopy 
of the candidates in order to confirm the structure, and plot it in three dimensions. 
Ultimately, further imaging at the ends of the filament would be interesting, to map 
out the structure completely. 



104 



Chapter 7 



• First detection of cold accretion: a Lya blob in GOODS-S 

In Chapter [2] we present the discovery of a Lya blob. This blob has an angular 
diameter of up to ~ 100 kpc and a Lya line luminosity of almost 10^^ erg s~^. After 
studying the complete set of deep multi-wavelength data available, we saw that this 
blob, unlike all other published blobs, had no obvious optical/infrared counter-part. 
It was realised, after taking all possible explanations into account, that this was 
the first observed Lya blob that could only be described by cold accretion in a 
satisfactory way. 

• ELVIS: a promising future survey for very high redshift Lya emitters 

A large part of the time during my PhD has been spent on the project which 
evolved into "ELVIS". The initial phase of the project included writing numerous 
applications for funding of the filters and observing time. I was also involved in the 
designing of the filters and the initial contacts with the filter producing companies. 
At later stages, I was involved in writing a Public Survey proposal and in the testing 
of the filters. After comparing the survey parameters (i.e. area, flux limit expected, 
filter specifications) of different near-infrared surveys for very high redshift Lya 
emitters, including ISAAC surveys, DaZle, ELVIS, JWST etc., it is clear that ELVIS 
will be one of the best opportunities to observe the very high redshift Universe. The 
combination of large area (0.9 deg^) and relatively deep flux limit for a near-infrared 
narrow-band image [Fu^^^^ = 3.7 x 10~^^ erg s~^ cm~^) should enable ELVIS to 
detect a few times ten Lya emitters at 2; = 8.8. These detections can be used 
to study galaxy formation and will probably give one of the best constraints on re- 
ionisation at this redshift from observed high redshift galaxies. Moreover, the survey 
will produce a vast catalogue of Ha emitters, as well as other emission line objects 
at intermediate redshift. This catalogue will be a treasure for anyone interested 
in studying the star formation and metallicity history of the Universe and galaxy 
evolution in general. Again, this sample will be gathered in the COSMOS field, 
a field with a vast amount of public, multi-wavelength data, including the very 
deep near-infrared data collected by Ultra- VISTA. This will enable further multi- 
wavelength studies of, in particular, the lower redshift galaxy sample. 

• Predictions for very high redshift Lya surveys 

In Chapter [5] we studied very high redshift Lya luminosity functions. This was 
a follow-up work from the preparatory work done for ELVIS. In this chapter, we 
extrapolate two published theoretical models for Lja emitters at lower redshifts 
to higher redshifts. We also try to extrapolate the parameters of observed Lya 
luminosity function fits to Schechter functions. The result is a plot of 7 < z < 12.5 
luminosity functions which can be used for groups planning future surveys for high 
redshift Lya. Using these luminosity functions, we made specific number predictions 
for several present and future narrow-band surveys. The conclusion of the project 
was that detections of small samples of z > 7 galaxies are to be expected in the next 
five years, but to get a good handle on the very high redshift Universe, yet more 
future observatories will be necessary. 
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• How to best select Lya emitters 

In the final chapter of this thesis, Chapter [6l we explore the colour space spanned 
by a population of mock galaxies, including "normal" galaxies, Lya emitters and 
[OII]-emitters. The aim of this work in progress is to determine how to best select 
Lya emitters from a general population of observed galaxies and how to distinguish 
them from interloper emitters such as [OII]-emitters. The motivation was to see if 
it is possible to select Lya emitters in a confident way so that spectroscopy is un- 
necessary, using the strength of present and future public multi-wavelength surveys 
such as GOODS and COSMOS. In our method, we calculated a "merit number" 
which quantifies how good any combination of colours is in separating two popula- 
tions of galaxies. This merit number is related to the projected distance between 
the two populations and to the area of the overlap between the two projections. 
As a result, we presented the three best selection combinations to distinguish Lya 
emitters from field galaxies or [Oil] emitters. The conclusion of the project is a 
set of selection colours and selection criteria that may be used by the community 
both to find Lya emitters with observed data, as well as when preparing observ- 
ing proposals. Two interesting results from the project were the conclusions that 
for selection, two colour combinations including the narrow-band magnitude offer 
better results than one colour combination and that these two combinations should 
include a broad-band filter on either side of the Lya line, but both blue- ward of the 
Balmer break. This is work that will be refined in the future. 

Summary 

To summarise, in the first part of the thesis, the emphasis was on using the available 
multi-wavelength data in GOODS-S to study the nature of Lya emitters at z ~ 3. 
This data turned out to be crucial in several aspects of my work. It was necessary, 
for instance, to determine that the Lya blob that was detected was indeed the first 
one to be best explained by cold accretion. This would have been impossible if the 
broad-band coverage had been smaller. The X-ray and mid-IR data ruled out an 
AGN nature of the blob and the optical/near-IR data excluded the possibility of a 
star-forming galaxy. Further, this data was also crucial in performing an accurate 
SED fitting for the candidates. It is true that many of the bands only included 
upper limits, but in many cases these limits constrained the SED. Especially the 
infrared data was useful in the SED fitting, in order to constrain the masses of 
the objects. It was also using the multi-wavelength data that the special SED of 
LEG0_G00DS-S#16 was identified as a potential Lya emitting dusty starburst 
galaxy. If this is confirmed, LEG0_G00DS-S#16 will be the first dusty starburst 
galaxy detected through its Lya emission. The increase in public multi-wavelength 
data also prompted the work presented in the final chapter, where we tried to find 
new selection criteria for Lya emitters at several different redshifts for surveys where 
a large complimentary data-set is available. 

In the second part of this thesis, a project to detect very high redshift {z = 8.8) 
Lya emitters was presented. This is a project which has not yielded any observa- 
tional results yet, although it prompted the work presented in Chapter [5l where 



106 



Chapter 7 



new Lya luminosity functions are derived and discussed in the context of future 
narrow-band surveys. ELVIS is a project which I am hoping to be part of for the 
coming years, and which will produce very interesting results within the next five 
years. 

The topic of Lja emitters is a field of research which is growing intensively at the 
moment, with increasing interest from the science community. It is a very promising 
method to easily detect high redshift, low luminosity, low mass, dust- and AGN-free 
star forming galaxies. It is also one of the most, if not the most, promising methods 
in detecting the highest redshift galaxies in the Universe. In the next chapter I write 
about what future ideas, plans and hopes I have regarding this subject. 



Chapter 8 

Future ideas, plans and hopes 



The subject of exploring the Universe using Lya as a tool has only begun, but it is 
increasing in popularity and it is quite imaginable that many new and interesting 
revelations about galaxy formation will come from this field in the future. A few 
aspects that I see as being crucial in understanding these galaxies, and which I hope 
to be able to work on are mentioned here. 

Understanding the stellar population in Lya emitters and their role in 
the scheme of high redshift galctxies. 

Today, the sample of confirmed or candidate Lya emitters is growing rapidly. The 
number of spectroscopically confirmed LEGOs at redshifts between z — 1.8 — 9 
approaches 600 and probably an equal amount of candidates have been presented in 
the literature. However, the exact nature of these galaxies is not yet known to any 
great detail. A few papers have attempted SED fitting of the continuum of LEGOs 
at redshifts ~ 3 (Gawiser et al. 2006; Nilsson et al. 2007), ~ 4.5 (Finkelstein et 
al. 2007) and ~ 5.7 (Lai et al. 2007). Unfortunately, the samples on which these 
analysis has been done are still very small and larger samples arc needed before 
statistically robust results will come. However, the future is promising. Several 
large public multi-wavelength surveys are under way. A project that will make a 
significant improvement in the field is the COSMOS survey at redshift z = 2.3 that 
we have started with the WFI imager in Chile. This study will greatly increase the 
number of LEGOs with sufficient multi-wavelength coverage for a large statistical 
study of the continuum SED. 

Another burning issue is that of the scheme of high redshift galaxies. There are 
today many different methods of finding high redshift galaxies, and just as many 
"classes" of objects; LEGOs, LBGs, DLAs, sub-mm galaxies, GRB host galaxies, 
DRGs, EROs etc. etc. How are these classes of galaxies related? It has been 
suggested that LEGOs are the predecessors of LBGs (Stark, Loeb & Ellis 2007) or 
that they are part of the same population (Giavalisco 2002), but neither suggestion 
has been proved. And where does the population of Lya emitting and dusty sub- 
mm galaxies fit in? During which phases of galaxy formation/evolution do galaxies 
emit Lya? These are very important questions that I hope will be, if not explained, 
than at least better understood in the next decade. 
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Understanding how Lya emitters trace Large Scale Structure and using 
it to derive information about the geometry of the Universe. 

LEGOs appear to be very good tracers of the Large Scale Structure (LSS) of the 
Universe. They have been observed in candidate proto-clusters (Venemans et al. 
2007 and references therein) and in filamentary structures (M0ller & Fynbo 2001; 
Matsuda et al. 2005; Nilsson et al. 2007). Are LEGOs more prone to follow LSS 
than other classes of galaxies? If so, why? If it is true that LEGOs are tracers 
of LSS, then they could be used as an independent test on cosmological models, 
as described by Weidinger et al. (2002). Further, large scale surveys such as the 
COSMOS survey and also specifically designed instruments such as the VIRUS 
instrument for the HETDE}4j survey may be able to make use of this property of 
LEGOs to draw conclusions about the cosmological parameters that describe our 
Universe. 

Understanding the evolution of Lya emitters with redshift. Are they the 
same population at different redshifts? 

As mentioned earlier, the sample of Lya emitting galaxies now extend between 
redshifts from 1.8 to almost 9. This covers almost 3 billion years in time, or almost 
25% of the age of the Universe. Are the galaxies that emit Lya the same across 
this period of time, or are they different types of objects? Are the mechanisms that 
create Lya the same? Yet another question is how LEGOs evolve with time, if they 
are only short-lived "fiares", or if it is one of many stages in galaxy formation. It 
would also be interesting to compare the properties of high redshift LEGOs with 
Lya emitting galaxies in the local Universe. Can we detect Lya emitters at redshift 
zero? Again, larger samples and better knowledge of the nature of these objects at 
different redshifts will hopefully add understanding within this subject. 

Exploring the star formation history of the Universe. 

Lya is a direct tracer of the star formation rate of a galaxy, albeit potentially 
affected by dust. Hence, using observations of Lya at different redshifts will put 
constraints on the star formation rate and star formation density history of the 
Universe. Together with observations from other emission-line galaxies, which are 
sometimes by-products of narrow-band surveys for Lya, and from other means of 
estimating star formation rates in galaxies, we can expect to get more and more 
accurate measurements of the star formation history of the Universe. 

Detecting the first stars, understanding the very first stages of galaxy 
formation and how and when re-ionisation occurred. 

LEGOs are truly one of the best tools in detecting the first stars and the first 
galaxies. As we move to higher and higher redshifts, the continuum of galaxies will 
become more and more faint and reach a level that is almost impossible to observe 
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with the instruments of today and to some extent of tomorrow as welL However, 
it has been suggested that the equivalent width of very high redshift galaxy Lja 
emission can be very high (Schaerer 2003; Dijkstra & Wyithe 2007) due to the 
Population III stars at those redshifts. It will then be possible to observe Lya 
emission to very high redshifts, even with todays instruments. ELVIS is one such 
example. With ELVIS we will hopefully discover a large sample of redshift z — 8.8 
galaxies and observe the first stages of galaxy formation. 

Re-ionisation is also a popular topic of study at the moment. When did it 
happen? How? Two sources of information about the re-ionisation are conceivable. 
One is mapping the Hll-emission (e.g. Ihev et al. 2002; Furlanetto, Zaldarriaga 
& Hernquist 2004; McQuinn et al. 2006) in the radio regime and the other one is 
Lya observations. With ELVIS we will hopefully be able to constrain the level of 
ionised gas at redshift z = 8.8, providing key information to help improve models 
of re-ionisation. 



In short, the future for Lya studies appears to be bright, with many interesting 
results to look forward to. I hope I can be part of them. 
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A Lyman-o; blob in the GOODS South field: evidence for cold accretion onto a 

dark matter halo 

by K.K. Nilsson, J.P.U Fynbo, P. M0ller, J. Sommer-Larsen & C. Ledoux 
Astronomy & Astrophysics, 452, L23-L26 (2006) 

In this paper, the observations, image reduction and source extraction were done by 
J. Fynbo (JF) and P. M0ller (PM). The spectroscopic observations were prepared 
by K. Nilsson (KN). The spectroscopic data was delivered bias subtracted, fiat- 
fielded and wavelength calibrated. KN did the spectrum extraction and stacking, 
sky subtraction, extraction of 1-D spectrum and fiux calibration. KN also did most 
of the analysis of the multi-wavelength data and the imaging/spectroscopic data 
available and the photometric redshifts. The calculations in the beginning of the 
Discussion were made by Jesper Sommer-Larsen (JSL), who also did the simulation 
of cold accretion described in the Discussion. The text in the paper was mainly 
written by KN with support from JF and PM, the plots were all made by KN. 
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A multi- wavelength study of z = 3.15 Lyman-a emitters in the GOODS South 

Field 

by K.K. Nilsson, P. M0ller, O. Moller, J.P.U. Fynbo, M.J Michalowski, D. Watson, 
C. Ledoux, P. Rosati, K. Pedersen & L.F. Grove 
Astronomy & Astrophysics, 471, 71-82 (2007) 

The results of this paper was based on the same narrow-band imaging as in Paper I, 
hence the original narrow-band image data reduction was performed by Johan Fynbo 
(JF) and Palle M0ller (PM). The source extraction was performed by Kim Nilsson 
(KN), as well as counterpart galaxy identification. The spectra were the same as 
in Paper I, i.e. basic data reductions were supplied from the ESO pipeline but KN 
performed stacking, sky subtraction, 1-D spectra extraction and flux calibration. All 
photometry in all bands was made by KN. The discovery of the filamentary structure 
and the statistical test was performed by KN. For the SED fitting, KN stacked the 
candidates and extracted photometry of the SED, and also gave input to the SED 
fitting code. However, the code itself was written by Ole Moller (OM). KN did most 
of the analysis on LEG0_G00DS-S#16, except for the GRASIL fits which were 
made by Michal Michalowski (MM). The comparison to LBGs was performed by 
KN. KN wrote most of the text, excluding small parts of the introduction, counter- 
part selection (section 2.2 and 2.3) and SED fitting method section (section 5.1). 
KN made all plots, excluding Fig. 8 and 12. 
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The two models (the semi-analytical and the phenomenological model) used in this 
paper have been presented in earlier publications and are not new to this work. 
The calculations made from the semi-analytical model were made by Alvaro Orsi 
(AO), Cedric Lacey (CL) and Carlton Baugh (CB) and the calculations from the 
Thommes & Meisenheimer model were made by Eduard Thommes (ET). Kim Nils- 
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the Schechter functions. The calculation for Cosmic Variance for ELVIS were made 
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which was written by AO. KN made all plots except Fig. 3 and 4. 
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Selection methods for Lyman-a emitters 
by K.K. Nilsson, O. MoUer, M. Hayes, P. M0ller & J.P.U. Fynbo 

in prep. 

In this paper, the mock sample of galaxies was created by Matthew Hayes (MH) 
and the code to calculate the merit number of each combination of selection colours 
was written by Ole Moller, all with input from Kim Nilsson (KN). KN wrote the 
text and made the plots. The selection criteria in Table 16.51 were also retrieved by 
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